# NOT RUN {
data("birthwt",package="MASS")
birthwt$Category.s <- ifelse(birthwt$bwt <= 2500, 3,
ifelse(birthwt$bwt <= 3000, 2,
ifelse(birthwt$bwt <= 3500, 1, 0)))
T.abs.mc <- rpartScore(Category.s ~ age + lwt + race + smoke +
ptl + ht + ui + ftv, data = birthwt)
plotcp(T.abs.mc)
T.abs.mc.pruned<-prune(T.abs.mc,cp=0.02)
plot(T.abs.mc.pruned)
text(T.abs.mc.pruned)
T.quad.mr <- rpartScore(Category.s ~ age + lwt + race + smoke + ptl + ht +
ui + ftv, split = "quad", prune = "mr", data = birthwt)
# }
Run the code above in your browser using DataLab