X <- matrix(c(4,2,99996,99998),nr=2) fisher.test(X, workspace=10000000) require(epitools) oddsratio(c(4,99996,2,99998),method="fisher") detach(package:epitools) require(vcd) OR <- oddsratio(X,log=F) ORL <- summary(oddsratio(X)) ORCI <- confint(OR) M <- c("オッズ比の点推定量"," (p=",") 95%信頼区間 = [ "," , "," ]\n") cat(M[1],OR,M[2],ORL[1,4],M[3],ORCI[1],M[4],ORCI[2],M[5])