底辺私立医大を卒業した医者って頭悪いよね？ Part32 (950ﾚｽ)
底辺私立医大を卒業した医者って頭悪いよね？ Part32 http://egg.5ch.net/test/read.cgi/hosp/1723605276/
上
下
前
次
1-
新
通常表示
512ﾊﾞｲﾄ分割
ﾚｽ栞
387: 卵の名無しさん [sage] 2024/09/12(木) 05:55:42.06 ID:Gblj4tVi 裏口容疑者には算出できない問題 こういう実地臨床に役立つ計算は楽しい。 " 正しく薬を服用すれば1回目の除菌療法の成功率は約68〜92％といわれています。 https://www.pylori-story.jp/sterilization/sterilization/point/ 尿素呼気試験の感度は90-100% 特異度80-99% https://i.imgur.com/KY052as.pnghttps://www.crc-group.co.jp/crc/q_and_a/67.html 上記数値はHighest Probability Densityでの95%信頼区間であると想定するとき 尿素呼気試験で陰性であった患者にピロリ菌感染が持続している確率とその95%信頼区間を求めよ。 " http://egg.5ch.net/test/read.cgi/hosp/1723605276/387
388: 卵の名無しさん [sage] 2024/09/12(木) 09:39:39.91 ID:Gblj4tVi # オミクロン株流行下の抗原定性検査の感度と特異度が判明 # https://www.aist.go.jp/aist_j/press_release/pr2023/pr20230130/pr20230130.html st2pv=\(prior,sl=0.53,su=0.73,tl=0.995,tu=1.00,cred=0.95){ # prior 検査前確率 感度[sl,su] 特異度[tl,lu] cred:信用レベル mv2ab<-function(m,v){# mean, variance to a,b a=(-m^3+m^2-m*v)/v b=(m^3-2*m^2+m*v+m-v)/v c(a=a,b=b) } absen=mv2ab((sl+su)/2,((su-sl)/2/qnorm(1/2+cred/2))^2) abspc=mv2ab((tl+tu)/2,((tu-tl)/2/qnorm(1/2+cred/2))^2) k=1e6 s=rbeta(k,absen[1],absen[2]) t=rbeta(k,abspc[1],abspc[2]) # 検査陽性の事後確率 postp=\(p,s,t) p*s/ (1-t+p*(s+t-1)) # p:事前確率 s:感度 t:特異度 # 検査陰性の事後確率 postn=\(p,s,t) p*(s-1)/(-t+p*(s+t-1)) # p:事前確率 s:感度 t:特異度 ppv=postp(prior,s,t) npv=postn(prior,s,t) ppv95=HDInterval::hdi(ppv,cred) npv95=HDInterval::hdi(npv,cred) cat('PPV\n') print(signif(ppv95[1:2],3)) cat('\nNPV\n') print(signif(npv95[1:2],3)) invisible(list(ppv,npv,ppv95,npv95)) } st2pv(0.5) http://egg.5ch.net/test/read.cgi/hosp/1723605276/388
389: 卵の名無しさん [sage] 2024/09/12(木) 10:15:18.07 ID:b25APUh4 >>387 ほとんどの数学は役に立つよ 役立たせ方を知らない、考えようとしてないだけ http://egg.5ch.net/test/read.cgi/hosp/1723605276/389
390: 卵の名無しさん [sage] 2024/09/12(木) 12:13:51.85 ID:3v9jdxfx 湘南鎌倉総合病院には元旭川医科大学勤務で国立大の弘前大学卒の元AV女優浅丘りなこと田中茉里子（旧姓古郡）がいるもんね https://www.skgh.jp/department/hbps/ https://www.suruga-ya.jp/product/detail/131150111 http://www.oma-aozora.jp/ikizama/kisamai/welcome_02.html http://egg.5ch.net/test/read.cgi/hosp/1723605276/390
391: 警備員[Lv.4][新芽] [] 2024/09/12(木) 19:29:36.15 ID:AgPCF0df >>387 脳内医者はダンマリ決め込むしかない問題>>327 http://egg.5ch.net/test/read.cgi/hosp/1723605276/391
392: 卵の名無しさん [sage] 2024/09/13(金) 14:31:38.52 ID:FlUFi/rO library(LearnBayes) # 累積密度関数の逆関数InvCDF(qbeta,qgamma等)と省略しない引数から # Highest Density Intervalの信頼区間を求める 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) } ab=beta.select(list(p=0.025,x=0.80),list(p=0.975,x=0.99)) ab curve(dbeta(x,ab[1],ab[2])) HDI(qbeta,shape1=ab[1],shape2=ab[2]) f=\(ab){ (pbeta(0.99,ab[1],ab[2]) - pbeta(0.80,ab[1],ab[2]) - 0.95)^2 } par=optim(par=ab,f)$par ; par HDI(qbeta,shape1=par[1],shape2=par[2]) http://egg.5ch.net/test/read.cgi/hosp/1723605276/392
393: 卵の名無しさん [sage] 2024/09/13(金) 16:58:13.01 ID:Ro8w2WW/ 湘南鎌倉総合病院には元旭川医科大学勤務で国立大の弘前大学卒の元AV女優浅丘りなこと田中茉里子（旧姓古郡）がいるもんね https://www.skgh.jp/department/hbps/ https://www.suruga-ya.jp/product/detail/131150111 http://www.oma-aozora.jp/ikizama/kisamai/welcome_02.html http://egg.5ch.net/test/read.cgi/hosp/1723605276/393
394: 卵の名無しさん [sage] 2024/09/14(土) 07:57:34.52 ID:4aTGKsaG # 信頼区間からβ分布の形状係数を算出する 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]) cat('lower =',ls[1],' upper =',ls[2],'\n') cat('α =',round(par[1],2),' β =',round(par[2],2),'\n') invisible(par) } http://egg.5ch.net/test/read.cgi/hosp/1723605276/394
395: 卵の名無しさん [sage] 2024/09/14(土) 08:24:14.38 ID:4aTGKsaG 改良 # 信頼区間からβ分布の形状係数を算出する 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){ 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]) 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') invisible(par) } ci2ab(0.80,0.99,0.95) http://egg.5ch.net/test/read.cgi/hosp/1723605276/395
396: 卵の名無しさん [sage] 2024/09/14(土) 08:48:01.20 ID:4aTGKsaG 臨床応用問題 尿素呼気試験 感度90-100%、特異度80-99% https://www.crc-group.co.jp/crc/q_and_a/67.html あるクリニックで 通算のP-CAB利用の除菌成功率は1次除菌が88.5％（502/567）であった。 判定は尿素呼気試験 感度90-100%、特異度80-99% https://www.crc-group.co.jp/crc/q_and_a/67.html で行った。 このクリニックで一次除菌成功と判定されたときにピロリ菌が持続感染している確率の 平均値、中央値、95%信頼区間を求めよ。 計算に必要な条件は適宜設定してよい。 例 尿素呼気試験の値は95%信頼区間の値であると設定して算出。 http://egg.5ch.net/test/read.cgi/hosp/1723605276/396
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
398: 卵の名無しさん [sage] 2024/09/14(土) 09:20:27.93 ID:YRPV5tb8 >>396 感度と特異度のpdfを描出 https://i.imgur.com/piuX3Rs.png https://i.imgur.com/hSYCE7O.png 区間の中心を平均値や中央値に設定しての算出より誤差が少なかった。 http://egg.5ch.net/test/read.cgi/hosp/1723605276/398
399: 卵の名無しさん [sage] 2024/09/14(土) 09:36:05.16 ID:6IJRlZll 湘南鎌倉総合病院には元旭川医科大学勤務で国立大の弘前大学卒の元AV女優浅丘りなこと田中茉里子（旧姓古郡）がいるもんね https://www.skgh.jp/department/hbps/ https://www.suruga-ya.jp/product/detail/131150111 http://www.oma-aozora.jp/ikizama/kisamai/welcome_02.html http://egg.5ch.net/test/read.cgi/hosp/1723605276/399
400: 卵の名無しさん [sage] 2024/09/14(土) 09:43:15.44 ID:4aTGKsaG 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 http://egg.5ch.net/test/read.cgi/hosp/1723605276/400
401: 警備員[Lv.4][新芽] [] 2024/09/14(土) 09:57:01.23 ID:9hxLr8jT 尿瓶ジジイ都合の悪いことは全てダンマリで板違いのチンパン数学()を垂れ流すことしかもうできない模様 http://egg.5ch.net/test/read.cgi/hosp/1723605276/401
402: 卵の名無しさん [] 2024/09/14(土) 10:42:48.87 ID:kTVOq2TW 底辺私立医大卒の医者は、なぜ国立医大に行かなかったのだろう⁈ 可哀想なのは患者だけで無く、底辺私立医大卒の家族たち!! 日々の何気ない生活の中で、時には知り合いから『お宅のご主人の大学は??』『おまえのお父さんの大学どこ??』という質問をうける事もあるだろう。 そんな時、阿保医者の家族達は『知らないんだょね』と誤魔化すのが精一杯、、知らない筈はないが… 全ての患者も、優秀な医者に診て欲しいと思っているから、底辺私立医大卒は寝る間を惜しんでもっと勉強、再教育すべきと思う ※国も医者のカテゴリ分けを検討すべきだな、、 国立医大卒はクラスA、 有名私立医大卒はクラスB、 その他底辺私立医大卒はクラスC、 って感じか？？？ http://egg.5ch.net/test/read.cgi/hosp/1723605276/402
403: 卵の名無しさん [sage] 2024/09/14(土) 11:03:08.76 ID:6Ky8Mu7c >>402 出身大学と医者のレベルは必ずしも一致しない 論文読めば分かるだろ http://egg.5ch.net/test/read.cgi/hosp/1723605276/403
404: 卵の名無しさん [sage] 2024/09/14(土) 11:35:02.46 ID:6IJRlZll 湘南鎌倉総合病院には元旭川医科大学勤務で国立大の弘前大学卒の元AV女優浅丘りなこと田中茉里子（旧姓古郡）がいるもんね https://www.skgh.jp/department/hbps/ https://www.suruga-ya.jp/product/detail/131150111 http://www.oma-aozora.jp/ikizama/kisamai/welcome_02.html http://egg.5ch.net/test/read.cgi/hosp/1723605276/404
405: 卵の名無しさん [sage] 2024/09/14(土) 12:04:35.95 ID:v9WwgBOq >>396 この計算から ピロリ菌再発は何人にひとりくらいと患者に説明すればよいか? http://egg.5ch.net/test/read.cgi/hosp/1723605276/405
406: 卵の名無しさん [sage] 2024/09/14(土) 12:06:13.71 ID:v9WwgBOq >>405 Rでの試算 # 信頼区間からβ分布の形状係数を算出する 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) } abs=ci2ab(0.90,1.00) abt=ci2ab(0.80,0.99) k=1e6 s=rbeta(k,abs[1],abs[2]) t=rbeta(k,abt[1],abt[2]) # 検査陽性の事後確率 postp=\(p,s,t) p*s/ (1-t+p*(s+t-1)) # p:事前確率 s:感度 t:特異度 # 検査陰性の事後確率 postn=\(p,s,t) p*(s-1)/(-t+p*(s+t-1)) # p:事前確率 s:感度 t:特異度 res=postn(1-502/567,s,t) 1/mean(res) 1/median(res) 普段、俺が説明している数字になったのが( ・∀・)ｲｲ!! http://egg.5ch.net/test/read.cgi/hosp/1723605276/406
上
下
前
次
1-
新
書
関
写
板
覧
索
設
栞
歴
あと 544 ﾚｽあります
ｽﾚ情報
赤ﾚｽ抽出
画像ﾚｽ抽出
歴の未読ｽﾚ
AAｻﾑﾈｲﾙ
Google検索
Wikipedia
ぬこの手
ぬこTOP
0.008s