source
## From: A two class data example Venables and Ripley pp. 194--199
library(MASS); data("birthwt"); attach(birthwt)
race <- (factor(race, labels = c("white", "black", "other")))
ptd <- factor(ptl > 0)
ftv <- factor(ftv)
levels(ftv)[-(1:2)] <- "2+"
bwt <- data.frame(low = factor(low), age, lwt, race, smoke = (smoke > 0)
, ptd, ht = (ht > 0), ui = (ui > 0), ftv)
detach(birthwt)
birthwt.glm <- glm(low ~ ., family=binomial(link=logit), data=bwt)
birthwt.step2 <- stepAIC(birthwt.glm, ~ .^2
+ I(scale(age)^2) + I(scale(lwt)^2), trace = F )
ht01.twoclass <- data.frame(observed = bwt$low
, predicted = predict(birthwt.step2
, type = "response"))
write.table(ht01.twoclass, file = "ht01.twoclass.txt")