laborki6.txt

----------
ZADANIE 1
----------
(a)

sluby=read.csv2("sluby.csv");
boxplot(sluby,names=c("Żona","Mšż"),col=c("red","blue"));

(b)

i<-25; #badana liczba
alfa<-0.05; #poziom istotnosci
N<-length(sluby$ZONA); #liczebnosc probki
q<-0.5; #mediana
a<-ceiling(N*q-qnorm(1-alfa/2)*sqrt(N*q*(1-q)));
b<-ceiling(N*q+qnorm(1-alfa/2)*sqrt(N*q*(1-q))); #przedzial ufnosci [a,b]

test<-function(){
if (sort(sluby$ZONA)[a]<=i && i<=sort(sluby$ZONA)[b])
cat("Liczba",i,"wpada do przedziału ufnosci dla zon\n");
if (sort(sluby$MAZ)[a]<=i && i<=sort(sluby$MAZ)[b])
cat("Liczba",i,"wpada do przedziału ufnosci dla mezow\n");
}
test();

----------
ZADANIE 2 
----------

wzrost=read.csv2("wzrost.csv");
hist(wzrost[,1],prob=T,15,col="yellow",main="Histogram wzrostu",xlab="Wzrost");
lines(density(wzrost[,1]),col="red",lwd=3);
x<-seq(130,220,by=0.1);
y<-dnorm(x,mean(wzrost[,1]),sd(wzrost[,1]));
lines(x,y,col="blue",lwd=3);

----------
ZADANIE 3
----------

(a)

k<-15;
dl<-(max(dane)-min(dane))/(k-1);
podzial<-dl*0:k-dl/2+min(dane);
table(cut(dane,podzial));
hist(dane,podzial,col="violet",main="Histogram wzrostu",xlab="Wzrost");

(b)

k<-6;
d<-signif(length(dane)/k,6);
podzial2<-sort(dane)[1+d*0:k];
hist(dane,podzial2,col="yellow",main="Histogram wzrostu",xlab="Wzrost");(density(wzrost[,1]),col="red",lwd=3);
x<-seq(130,220,by=0.1);
y<-dnorm(x,mean(wzrost[,1]),sd(wzrost[,1]));
lines(x,y,col="blue",lwd=3);

----------
ZADANIE 3
----------

(a)

k<-15;
dl<-(max(dane)-min(dane))/(k-