底辺私立医大を卒業した医者って頭悪いよね? Part32 (902レス)
底辺私立医大を卒業した医者って頭悪いよね? Part32 http://egg.5ch.net/test/read.cgi/hosp/1723605276/
上
下
前
次
1-
新
通常表示
512バイト分割
レス栞
397: 卵の名無しさん [sage] 2024/09/14(土) 09:04:20.00 ID:4aTGKsaG 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) } http://egg.5ch.net/test/read.cgi/hosp/1723605276/397
上
下
前
次
1-
新
書
関
写
板
覧
索
設
栞
歴
あと 505 レスあります
スレ情報
赤レス抽出
画像レス抽出
歴の未読スレ
AAサムネイル
Google検索
Wikipedia
ぬこの手
ぬこTOP
0.006s