# Use the ovarian cancer data
data(Xdata, package="CGEN")
# Fit using a stratification variable
fit <- snp.logistic(Xdata, "case.control", "BRCA.status",
main.vars=c("oral.years", "n.children"),
int.vars=c("oral.years", "n.children"),
strata.var="ethnic.group")
# Compute the effects
effects <- snp.effects(fit, "oral.years", var.levels=c(0, 2, 3))
printEffects(effects)
Run the code above in your browser using DataLab