# 5. přednáška library(RcmdrPlugin.TeachingDemos) pnorm((140.5-136.1)/6.4)-pnorm((134.5-136.1)/6.4) pnorm(140.5,136.1,6.4)-pnorm(134.5,136.1,6.4) # běžně používané kritické hodnoty (alpha = c(0.5,0.25,0.10,0.05,0.01)) qnorm(1-alpha/2) qt(1-alpha/2,100) qt(1-alpha/2,20) qt(1-alpha/2,5) # aproximace binomického normálním x = 0:60 # celá čísla od 0 do 60 xx = seq(0,60,length=201) # 201 hodnot mezi 0 a 60 plot(x,dbinom(x,60,1/6),col=2,pch=1,xlab="k",ylab="",main="") lines(xx,dnorm(xx,10,sqrt(50/6)),col=2) points(x,dbinom(x,60,3/6),col=3,pch=2) lines(xx,dnorm(xx,30,sqrt(90/6)),col=3) points(x,dbinom(x,60,4/6),col=4,pch=3) lines(xx,dnorm(xx,40,sqrt(80/6)),col=4) legend("topright",legend=c("bi(60,1/6)","bi(60,3/6)","bi(60,4/6)"), lty=1,col=2:4,pch=1:3) # data(simulaceTisic) # opar <- par(mfrow = c(2,2), oma = c(1.1, 0, 1.1, 0), las = 1) hist(x,xlim=range(x),main="populace",xlab="",col="yellow") hist(x1,xlim=range(x),main="n=1",xlab="",col="yellow") hist(x10,xlim=range(x),main="n=10",xlab="",col="yellow") hist(x100,xlim=range(x),main="n=100",xlab="",col="yellow") par(opar) # opar <- par(mfrow = c(2,2), oma = c(1.1, 0, 1.1, 0), las = 1) hist(x,xlim=range(x),main="populace",xlab="",col="yellow") hist(x1,main="n=1",xlab="",col="yellow") hist(x10,main="n=10",xlab="",col="yellow") hist(x100,main="n=100",xlab="",col="yellow") par(par)