[過去ﾛｸﾞ] 当直医のスレ Part24 [無断転載禁止]©2ch.net (1002ﾚｽ)
上下前次1-新

---

このｽﾚｯﾄﾞは過去ﾛｸﾞ倉庫に格納されています｡
次ｽﾚ検索 歴削→次ｽﾚ 栞削→次ｽﾚ 過去ﾛｸﾞﾒﾆｭｰ

---

722: 2017/10/08(日)20:04 ID:jiwdTQ4W(11/21) AA×


[240|320|480|600|100%|JPG|べ|ﾚｽ栞|ﾚｽ消]
722: [sage] 2017/10/08(日) 20:04:54.22 ID:jiwdTQ4W ## くじの数を増やせば一様分布になるか？ Takarakuji <- function(k=10^2,.Na=1000,.Nb=1000,.pa=0.1,.pb=0.01, .n=100,.Wa=100,.Wb=1000,Histogram=TRUE,Print=TRUE,...){ Na=.Na ; Nb=.Nb pa=.pa ; pb=.pb Wa=.Wa ; Wb=.Wb n=.n aa=0:n A0=c(rep(1,Na*pa),rep(0,Na-Na*pa)) B0=c(rep(1,Nb*pb),rep(0,Nb-Nb*pb)) Get <- function(a)sum(sample(A0,a))*Wa + sum(sample(B0,n-a))*Wb MAX=replicate(k,which.max(sapply(aa,Get))-1) tbl_MAX=table(MAX) MODE=tbl_MAX[which.max(tbl_MAX)] mode=as.numeric(names(MODE)) if(Histogram) { hist(MAX,freq=FALSE,main=paste('Na=',Na,'Nb=',Nb),...) lines(density(MAX),lwd=1) } if(Print){ print(round(c(mean=mean(MAX),sd=sd(MAX),mode=mode),2)) } invisible(c(mean=mean(MAX),sd=sd(MAX),mode=mode)) } par(mfrow=c(2,2)) Takarakuji(.Na=10^3,.Nb=10^3,col='lightblue') Takarakuji(.Na=10^4,.Nb=10^4,col='wheat') Takarakuji(.Na=10^5,.Nb=10^5,col='maroon') Takarakuji(.Na=10^6,.Nb=10^6,col='gray') # どれも一様分布とは程遠い http://egg.5ch.net/test/read.cgi/hosp/1506113587/722
くじの数を増やせば一様分布になるか？ どれも一様分布とは程遠い

---

上下前次1-新書関写板覧索設栞歴

---

あと 280 ﾚｽあります
ｽﾚ情報 赤ﾚｽ抽出 画像ﾚｽ抽出 歴の未読ｽﾚ AAｻﾑﾈｲﾙ

---

ぬこの手 ぬこTOP 0.037s