#kast med terning x=sample(c(1:6),100,replace=T) x[1:10] mean(x) #terning med 3 sider middelvaerdi og varians #15 kast x=sample(c(1:3),15,replace=T) 2 1 2 3 3 2 1 2 2 3 1 1 1 2 1 x=sample(c(1:3),10000000,replace=T) x mean(x)#empiriske gennemsnit 2.3 var(x)#s^2 0.68 #kast med 5 terninger x=sample(c(1:6),5,replace=T) x sum(x) #min(x) #c(1:6) #c(1,2,3,4,5,6) totalsum=rep(0,1000)#just a vector of 0's #totalsum #minoejne=rep(0,1000) for (i in 1:1000){ x=sample(c(1:6),5,replace=T) totalsum[i]=sum(x) #minoejne[i]=min(x) } totalsum[1:10] mean(totalsum) 17.51+2*sd(totalsum)/sqrt(100000) binare=as.numeric(totalsum>=15) mean(binare) mean(binare)+2*sd(binare)/sqrt(100000) sum(totalsum>=15) 773/1000 mean(totalsum>=15) 78103/100000 mean(totalsum) 17.445 s=sd(totalsum) 3.70 s s/sqrt(10) sum(totalsum>=15)/1000 z=minoejne>2 I=as.numeric(z) sum(I)/1000 mean(I) sd(I) 0.3386 sqrt(0.132*(1-0.132)) sd(I)/sqrt(100000) totalsum[1:10] mean(totalsum)#eksakt svar er 5*3.5=17.5 var(totalsum) 14.3975 sd(totalsum) 3.79 17.5016-2*3.806271/sqrt(100000) 17.5016+2*3.806271/sqrt(100000) z=totalsum>14 I=as.numeric(totalsum>14) mean(I) 0.7788 0.793 var(I) 0.1643 mean(I)-2*sd(I)/sqrt(100000) mean(I)+2*sd(I)/sqrt(100000)