密度症例対照研究のシミュレーションまとめ 2 各パラメータにおける分散についてグラフにしてみた。◎e1についてe1が0.5に近い値をとっているときのほうが分散の値が小さいことは確かなようである。◎p/qについてp/qの値が大きいほうが分散の値は大きくなった…

累積症例対照研究のシミュレーション 2 関数のバグはとれた。 mat<-function(N,e,p,q){ a<-N*10000 b<-e*10000 r<-runif(a,0,1) r[1:b][r[1:b]>=1-p]<-2 r[(b+1):a][r[(b+1):a]>=1-q]<-2 s<-rep(c(1,0),c(b,a-b)) return(rbind(r,s)) } box<-function(N,e,p,…