底辺私立医大を卒業した医者って頭悪いよね? Part32 (867レス)
前次1-
抽出解除 必死チェッカー(本家) (べ) レス栞 あぼーん

394: 09/14(土)07:57 ID:4aTGKsaG(1/5) AAS
# 信頼区間からβ分布の形状係数を算出する
ci2ab=\(l,u,cl=0.95){ # ci : [l,u], cl : confidence level
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),f1)
opt=optim(opt$par,f1)
par=opt$par
ls=HDI(qbeta,shape1=par[1],shape2=par[2])
省4
395: 09/14(土)08:24 ID:4aTGKsaG(2/5) AAS
改良

# 信頼区間からβ分布の形状係数を算出する
ci2ab=\(l,u,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){
省12
396
(2): 09/14(土)08:48 ID:4aTGKsaG(3/5) AAS
臨床応用問題

尿素呼気試験 感度90-100%、特異度80-99%
外部リンク[html]:www.crc-group.co.jp

あるクリニックで
通算のP-CAB利用の除菌成功率は1次除菌が88.5%(502/567)であった。
判定は尿素呼気試験 感度90-100%、特異度80-99%
外部リンク[html]:www.crc-group.co.jp
で行った。

このクリニックで一次除菌成功と判定されたときにピロリ菌が持続感染している確率の
平均値、中央値、95%信頼区間を求めよ。
省2
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){
省15
400: 09/14(土)09:43 ID:4aTGKsaG(5/5) AAS
LearnBayesのbeta.select関数(分位数から形状母数を算出)だと
エントリーに0や1があると計算してくれない。
95信頼区間を2.5%-97.5%の分位数として算出させようとするとエラーが返ってくる。
> LearnBayes::beta.select(list(p=0.025,x=0.90),list(p=0.975,x=1.00))
Error in approx(prob2[ind], logK[ind], p2) :
need at least two non-NA values to interpolate
前次1-
スレ情報 赤レス抽出 画像レス抽出 歴の未読スレ AAサムネイル

ぬこの手 ぬこTOP 1.478s*