風俗嬢ですが再起して頑張り医学部に合格しました [転載禁止]©2ch.net

**卵の名無しさん** · 2015/04/06(月) 10:00:54.93

風俗板にて

風俗辞めた後に医学部に合格したと嘘をついて袋叩きにされてる様子
http://aoi.bbspink.com/test/read.cgi/club/1424999821/724
760 ：名無しさん＠ピンキー：2015/04/03(金) 17:34:23.58 ID:hEZXryQ10
>>724
アホか、お前
＞医者になった嬢も居る、
そんなの絶対いないわ、ボケ。
お前は医者と看護師の区別もつかんのか(大笑)

**卵の名無しさん** · 2017/05/04(木) 18:33:14.57

//2x2.stan
data{
int<lower=0> N[2]; //行毎の標本和（(介入+対照)）
int n[2,2];// 行：介入の有無、列：結果の有無
}

parameters{
simplex[2] p[2];//ex. 総和sum(p[1]=1)となる確率のベクトル
}

model{
for(i in 1:2){//行
for(j in 1:2){//列
n[i,j]~ binomial(N[i],p[i][j]);//ex. n[1,2] ~ Binom(N[1],p[1][2])
}//行ごと(=介入群、対照群ごと）に二項分布
}
/*
n[1,1] ~ binomial(N[1],p[1][1]);
n[1,2] ~ binomial(N[1],p[1][2]);
n[2,1] ~ binomial(N[2],p[2][1]);
n[2,2] ~ binomial(N[2],p[2][2]);
*/
}

**卵の名無しさん** · 2017/05/04(木) 18:33:21.26

generated quantities{
real d;
real delta_over;
real p11;
real p10;
real p01;
real p00;
real RR;
real OR;
p11 = p[1][1];//介入あり効果あり、(暴露あり病気あり)
p10 = p[1][2];//介入あり効果なし、(暴露あり病気なし)
p01 = p[2][1];//介入なし効果あり、(暴露なし病気あり)
p00 = p[2][2];//介入なし効果なし、(暴露なし病気なし)
d = p11 - p01;//暴露あり病気あり割合－暴露なし病気あり割合
delta_over = step(d);// ifelse(d>0,1,0)に相当
RR = p11/p01;
OR = (p11/p10)/(p01/p00);
}

**卵の名無しさん** · 2017/05/04(木) 18:33:49.02

library("rstan")
rstan_options(auto_write = TRUE)
options(mc.cores = parallel::detectCores())

options(scipen = 10)

# UC Berkeley gender bias
# Men 8442 44%
# Women 4321 35%

n = matrix(round(c(8442*.44,8442*.56,4321*.35,4321*.65)),ncol=2,byrow=TRUE) ; n
N=apply(n,1,sum) ; N
data <- list(n=n,N=N)
model2x2=stan_model('2x2.stan')
fit2x2 <- sampling(model2x2,data=data,seed=1234)

print(fit2x2,pars=c('d','RR','OR'),digits_summary=dig,probs=c(.025,.975))
stan_dens(fit2x2,pars=c('d','RR','OR'),separate_chains = TRUE)

Epi::twoby2(n)
fisher.test(n)
prop.test(n[,1],c(8442,4321),correct=FALSE)