TOP カテ一覧 スレ一覧 100〜終まで 2ch元 削除依頼
★難関 糖尿病内科専門医対策スレッド２
放射線科総合(１例目)
医学界の在日朝鮮人
臨床統計もおもしろいですよ、その3
患者と付き合ったことのある医者っているの？
慈恵医大について語ろう
「続 」安定剤 睡眠薬の恐怖
盗撮医師 多田素久の盗撮（千葉大学卒）
医学部医学科ランキング
【誤診】千葉県のヤブ医者スレ【手遅れ】
- 855 ：
- これで、
# 2x2のクロス表からstanを使ってMCMCサンプリングしてFisherテストしてp.valueの配列からその分布を表示させる
library(rstan)
options(mc.cores = parallel::detectCores())
rstan_options(auto_write = TRUE)
contingency2stan2pf <- function(mat,alpha=0.05,print=T){
N=apply(mat,1,sum)
data <- list(mat=mat,N=N)
stanModel=readRDS('contingency.rds')
fit=sampling(stanModel,data=data, chain=1, iter=25000,warmup=5000)
pars=c('y1','n1_y1','y0','n0_y0','d','RR','OR','NNT')
if(print) print(fit, pars=pars,digits=4,prob=c(0.025,0.5,0.975))
ms=rstan::extract(fit)
p.sim=NULL
for(i in 1:length(p11)){
MAT=round(with(ms,matrix(c(y1[i],n1_y1[i],y0[i],n0_y0[i]),ncol=2,byrow=TRUE)))
p.sim=append(p.sim,fisher.test(MAT)$p.value)
}
BEST::plotPost(p.sim,compVal = 0.05, xlab='p.value')
cat('Pr(p.value < ',alpha,') = ', mean(p.sim<alpha),'\n')
cat('Original p.value = ',fisher.test(mat)$p.value,'\n')
print(summary(p.sim))
invisible(p.sim)
}
サンプル数を増やしてｐを小さくして信用させるという手法を見破れるかな。
広大2年「神経解剖学」追試１２０人「全員不合格」
【愛仁会】尼崎だいもつ病院
美人でH〜Iカップの女医女子医大生がいない理由
２ちゃんねる医道審議会
療養病床勤務医・地域包括ケア勤務医・老健勤務医が集うスレ☆Part.6
☆★★芸能界はフランクリン界
☆★ 札幌ひばりが丘病院 医師 薬剤師 看護師の免許取消、逮捕 情報集約スレ ★★
NTT東日本関東病院
【過去】医者板で話題になった女性たち【現在】
【コロナ収束第二波】マジでギリギリな開業医【患者増減】9
--------------------
[Ampere]NVIDIA GeForce RTX30XX総合 Part3
■長澤まさみ と 野田洋次郎の噂 ﾜﾝｵｸtakaの噂★30
非処女だけど婚活相手とは結婚するまでセックスはしたくない女
（●＾o＾●）うまいもん（●＾o＾●）
TBS★近藤夏子★新人
【ポケ森】キャンプ場おひろめスレ 〜 転載禁止 part7
モテないけど結婚相談所に入った
一般人には見分けがつかないアニメキャラ貼ってけwwww
【ペコちゃん】岡崎朋美選手専用スレ【万年３位】
自慢の腕時計コレクションをうpするスレ
ISFnet Japan｜(株)アイエスエフネット
LAMEコマンドラインオプションを語れ！その42
【PSP/VITA】DJMAX Portable 60曲目
今日の朝ご飯を書き込むスレ 5
【天体観測Diner】アフィリア・スターズ Part.03
色っぽい、色気がある人の特徴
【NMB48】大段舞依応援スレ★9【まいち】
明石商 ｖｓ．履正社 （令和元年夏近畿頂上決戦）
「トロフィー不要論者＝ゲームにやりこみなど不要と主張する人」だけど、それは違うと思うなぁ。
【妖怪ウォッチ】キャラクターアンチスレ
TOP カテ一覧 スレ一覧 100〜終まで 2ch元 削除依頼
