| library(Design) | |||||||
| reri <- function(datsam) | |||||||
| { | |||||||
| fitlr <- glm(outcome ~ determinant A * determinant B, family=binomial, data=datsam) | |||||||
| reri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1 | |||||||
| } | |||||||
| summary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975)) | |||||||


雷达卡





京公网安备 11010802022788号







