TOP カテ一覧 スレ一覧 100〜終まで 2ch元 削除依頼
数学の洋書vol.2
現代数学の系譜 工学物理雑談 古典ガロア理論も読む64
素人には 8÷2(2+2) を16と答える馬鹿が居るらしい
ポエムはここに書いてね 2
くだらねぇ問題はここへ書け
不等式への招待 第9章
無職だから最近数学の勉強をしている
ゲーデルの不完全性定理を論破する方法
現代数学はインチキのデパート
【大学院へ】 30過ぎて、数学の道へ その2
17 :
PCR検査の感度を0.7、特異度を0.9とする。
広島県で第一号の感染発見例は
県の検査で1回陰性、病院の検査で2回陽性、症状軽快した現時点で陰性(何回やったか報道がないので1回陰性とする)であるという。
ここで問題:
検査前の感染確率の分布が一様分布であると仮定して、
現在患者が感染している確率とその95%CIを計算してみた。
"

https://i.imgur.com/4CzTEWA.jpg

sn=0.7 # sensitivity
sp=0.9 # specificity
plus=2 # how many positive result?
minus=2 # how many negative result?
n=1e7 # how large the simulation
p0=runif(n,0,1)
oz0=p0/(1-p0) # prob -> odds
pLR=sn/(1-sp) # TP/FP
nLR=(1-sn)/sp # FN/TN
oz1=oz0*pLR^plus*nLR^minus # Bayesian formula
p1=oz1/(1+oz1) # odds -> prob
BEST::plotPost(p1,showMode =T) # show mode instead of mean
BEST::plotPost(p1,showMode =F)
HDInterval::hdi(p1) # Highest Density Interval
quantile(p1,c(.025,0.5,.975)) # 95%CI by quantile
summary(p1) # mean, median
MAP <- function(x) {
dens <- density(x)
mode_i <- which.max(dens$y)
mode_x <- dens$x[mode_i]
mode_y <- dens$y[mode_i]
c(x=mode_x, y=mode_y)
}
MAP(p1)['x'] # show mode

数学博士中退後就職できません。ご意見ください。
点の長さは0である。しかし点の集まりである線には長さが存在する。
雑談はここにかけ!【54】
数学の本 第86巻
プラスマイナスゼロ
人工知能・機械学習のスレッド@数学板
数学って結局暗記ゲーなの?
現代数学の系譜 古典ガロア理論を読む36
■■■■■■■■■■■■■   人工太陽
【未解決問題】奇数の完全数が存在しないことの証明3
--------------------
旧世界二大悪 ナチスドイツ帝国と大日本帝国
Google検索エンジンの劣化について語るPart1
芸能人格付けチェック! PART9
今日買ったゲームソフトを報告するスレ Part17
メトロンブログを語るスレ▼▼2
【DAXEL】結城友奈は勇者である12玉目【総合】
【マイクロ】MicroATX PC 総合 27台目【M-ATX】
月ノ美兎応援スレ
ミ☆魔法の天使 クリィミーマミ 33☆彡
☆★☆★☆神なんて実在しないよね(*´・ω・)(・ω・`*)ネ- 339
【速報】PS5リークが来たぞ
オヴェルスの翼
[test] 書き込みテスト 専用スレッド [てすと]
ラッキーセブンディスタンス陸上競技大会
【CLUB】クラブパナソニック 10pt【Panasonic】
ヤマザキ 豆大福風パンが大好きな人が集まるスレ!
太田上田 Part10
新規開業クリニック(動物病院版) 9件目
■自作・オリジナルPCケース Vol52■
買ったソフトを報告して語ろう 3
TOP カテ一覧 スレ一覧 100〜終まで 2ch元 削除依頼