> or<-function(p,q,r,s){data<-matrix(c(p,q,p+q,r,s,r+s,p+r,q+s,p+q+r+s), nr=3) + colnames(data)<-c("Before","After","Total") + rownames(data)<-c("sexual","other","Total") + print(data) + OR<-q*r/(p*s) + RRL<-OR*exp(-qnorm(0.975)*sqrt(1/p+1/q+1/r+1/s)) + RRU<-OR*exp(qnorm(0.975)*sqrt(1/p+1/q+1/r+1/s)) + cat("OR:",OR," [95%CI:",RRL," - ",RRU," ]\n")} > or(107,33557,255,45702)ただ、WaiWaiへの地域別アクセス数のデータがあれば、地域別被害件数の増加と比較して、統計学的な有罪判決が出せるかもしれませんね。
Before After Total
sexual 107 255 362
other 33557 45702 79259
Total 33664 45957 79621
OR: 1.749864 [95%CI: 1.395625 - 2.194017 ]
それにしても、「毎」も「売」もハングルでは、"매"と読まれるようで、何たる皮肉!
>毎日新聞は「売日新聞」に紙名を変更するとよい。
0 件のコメント:
コメントを投稿