分からない問題はここに書いてね457

分からない問題はここ ..

5:ﾄみた。同様に確からしければ1/90=0.01111=1.111111%の近似値が返ってくるはず。 # 同様に確からしくないシミュレーション rm(list=ls()) r12=c(1,1,3,2,2,3) r13=c(1,2,1,3,2,3) r23=c(1,2,2,3,1,3) r56=c(2,1,1,2,3,3) sim <- function(x,n=6){ # 部屋の割当がxと等しいかT/Fを返す r=numeric(n) # １～6人の部屋番号（１～３）の配列 for(i in 1:n){ j=which(c(sum(r==1)<n/3 , sum(r==2)<n/3 , sum(r==3)<n/3 ))　#　定員に達していない部屋から r[i]=as.numeric(sample(as.character(j),1)) # ランダムに割り当てる } all(x==r) } mean(replicate(1e5,sim(r12)))*100 # 百分率%表示 mean(replicate(1e5,sim(r13)))*100 mean(replicate(1e5,sim(r23)))*100 mean(replicate(1e5,sim(r56)))*100 実行結果 > mean(replicate(1e5,sim(r12)))*100 # 百分率%表示 [1] 1.384 > mean(replicate(1e5,sim(r13)))*100 [1] 0.949 > mean(replicate(1e5,sim(r23)))*100 [1] 0.915 > mean(replicate(1e5,sim(r56)))*100 [1] 1.874

次ページ

続きを表示

1を表示