底辺私立医大を卒業した医者って頭悪いよね？ Part32 (897ﾚｽ)
上下前次1-新

---

397: 09/14(土)09:04 ID:4aTGKsaG(4/5) AAS
pdfを描出するように改良

# 信頼区間からβ分布の形状係数を算出する
ci2ab=\(l,u,verbose=FALSE,cl=0.95){ # ci : [l,u], cl : confidence level
HDI=\(InvCDF=qbeta,cred=0.95,...){
opt=optimize(\(p) InvCDF(p+cred,...) - InvCDF(p,...),c(0,1-cred))
lwr=InvCDF(opt$min,...)
upr=lwr+opt$obj
c(lwr,upr)
}
f=\(ab){
LU=HDI(qbeta,cred=cl,shape1=ab[1],shape2=ab[2])
(LU[1]-l)^2 + (LU[2]-u)^2
}
opt=optim(runif(2,1,100),f)
opt=optim(opt$par,f)
par=opt$par
ls=HDI(qbeta,cred=cl,shape1=par[1],shape2=par[2])
if(verbose){
cat('mean =',round(par[1]/sum(par),3),'lower =',round(ls[1],3),' upper =',round(ls[2],3),'\n')
cat('α =',round(par[1],2),' β =',round(par[2],2),'\n')
curve(dbeta(x,par[1],par[2]),type='h',col=2,n=250,bty='l',ann=FALSE,axes=FALSE)
axis(1)
}
invisible(par)
}

---

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

---

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

---

ぬこの手 ぬこTOP 0.008s