[過去ﾛｸﾞ] 臨床統計もおもしろいですよ、その１ [無断転載禁止]©2ch.net (747ﾚｽ)

臨床統計もおもしろいですよ、その１ [無断転載禁止]©2ch.net http://egg.5ch.net/test/read.cgi/hosp/1493809494/

上下前次1-新
通常表示 512ﾊﾞｲﾄ分割 ﾚｽ栞
抽出解除 必死ﾁｪｯｶｰ(本家) (べ) ﾚｽ栞 あぼーん

---

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

---

269: 卵の名無しさん [sage] 2017/11/23(木) 15:22:14 ID:Ue5tZuwc par(mfrow=c(2,2)) theta=0.5 NN=1000 N=1:NN flip=numeric() for(i in N){ flip=append(flip,rbinom(1,1,theta)) } z=cumsum(flip) z_N=z/N plot(z_N,type='l',ylim=c(0,1)) pv=numeric() for(i in N){ pv[i]=binom.test(z[i],N[i],theta)$p.value } plot(pv,type='l') abline(h=0.05,col='blue',ylim=c(0,1)) bf=numeric() for(i in N){ bf[i]=beta(z[i]+1,N[i]-z[i]+1)/beta(1,1)/(theta^z[i]*(1-theta)^(N[i]-z[i])) } plot(log(bf),type='l') abline(h=log(3),lty=3) abline(h=log(1/3),lty=3) http://egg.5ch.net/test/read.cgi/hosp/1493809494/269

270: 卵の名無しさん [sage] 2017/11/23(木) 15:22:43 ID:Ue5tZuwc hdi=NULL for(i in N){ y=flip[1:i] s=rep(1,i) data=data.frame(y,s) Ntotal=i Nsubj=1 dataList=list(y=y,s=s,Ntotal=Ntotal,Nsubj=Nsubj) js=genMCMC2(data) hdi=rbind(hdi,HDIofMCMC(as.matrix(js))) } saveRDS(hdi,'hdi_sequential') hdi=readRDS('hdi_sequential') plot(hdi[,1],type='l',ylim=c(0,1),main='95% HDI') lines(hdi[,2]) abline(h=0.5,col=4) plot(apply(hdi,1,diff),type='l',main='HDI width') http://egg.5ch.net/test/read.cgi/hosp/1493809494/270

---

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

---

ｽﾚ情報 赤ﾚｽ抽出 画像ﾚｽ抽出 歴の未読ｽﾚ AAｻﾑﾈｲﾙ

---

Google検索

Wikipedia

---

ぬこの手 ぬこTOP 1.544s*