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

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

上下前次1-新
通常表示 512ﾊﾞｲﾄ分割 ﾚｽ栞

---

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

---

255: 卵の名無しさん [sage] 2017/11/15(水) 10:17:36 ID:4qVYF7j+ y= c(1,1,1,1,1,1,1,1,1) N=length(y) omega1=0.75 omega2=0.25 kappa1=12 kappa2=12 modelString=' model { for (i in 1:N){ y[i] ~ dbern( theta ) } theta <- equals(m,1)*theta1 + equals(m,2)*theta2 theta1 ~ dbeta( omega1*(kappa1-2)+1 , (1-omega1)*(kappa1-2)+1 ) theta2 ~ dbeta( omega2*(kappa2-2)+1 , (1-omega2)*(kappa2-2)+1 ) m ~ dcat( mPriorProb[] ) mPriorProb[1] <- p mPriorProb[2] <- 1-p p ~ dunif(0,1) } ' dataList=list(y=y,N=N, omega1=omega1,omega2=omega2, kappa1=kappa1, kappa2=kappa2) writeLines( modelString , con="TEMPmodel.txt" ) jagsModel = jags.model( file="TEMPmodel.txt" , data=dataList) update( jagsModel) codaSamples = coda.samples( jagsModel , variable=c('m','theta','theta1','theta2'), n.iter=50000) summary(codaSamples) js=as.matrix(codaSamples) tapply(js[,'theta'],js[,'m'],length) sum(js[,'m']==1)/nrow(js) http://egg.5ch.net/test/read.cgi/hosp/1493809494/255

---

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

---

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

---

Google検索

Wikipedia

---

ぬこの手 ぬこTOP 0.173s*