if (FALSE) {
##load the data set
data(data_pltr)
## set the parameters
args.rpart <- list(minbucket = 40, cp = 0)
family <- "binomial"
Y.name <- "Y"
X.names <- "G1"
G.names <- paste("G", 2:15, sep="")
## build a maximal tree
fit_pltr <- pltr.glm(data_pltr, Y.name, X.names, G.names, args.rpart = args.rpart,
family = family,iterMax = 5, iterMin = 3)
## Coerce a tree into a glm model using the confonding factor
fit_glm <- tree2glm(fit_pltr$tree, data_pltr, Y.name, X.names,
family = family)
summary(fit_glm)
}
Run the code above in your browser using DataLab