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)) |