## Not run:
# set.seed(13)
# library(envlpaster)
# library(aster2)
# data(generateddata)
# m.null <- aster(resp ~ 0 + varb, fam = fam, pred = pred,
# varvar = varb, idvar = id, root = root, data = redata)
# m1 <- aster(resp ~ 0 + varb + mass + timing,
# fam = fam, pred = pred, varvar = varb, idvar = id,
# root = root, data = redata)
# m2 <- aster(resp ~ 0 + varb + mass + timing +
# I(mass^2) + I(timing^2) + I(mass*timing),
# fam = fam, pred = pred, varvar = varb, idvar = id,
# root = root, data = redata)
# anova.table <- anova(m.null,m1,m2); anova.table
# beta <- m1$coef
# a <- grepl( "offsp", names(beta))
# a <- a + grepl( "surviv", names(beta))
# b <- which(a == 1)
# target <- c(1:length(beta))[-b]
# nnode <- ncol(m1$x)
# data.aster <- asterdata(data, vars, pred, rep(0,nnode),
# fam, families = list("bernoulli", "poisson",
# fam.zero.truncated.poisson()))
# selection(parm = beta, index = target, model = m1,
# data = data.aster, alpha = 0.05, type = "canonical",
# method = "eigen")
# ## End(Not run)
Run the code above in your browser using DataLab