# NOT RUN {
data(regularity)
# ---- predicting regularity with a logistic regression model
library(rms)
regularity.dd = datadist(regularity)
options(datadist = 'regularity.dd')
regularity.lrm = lrm(Regularity ~ WrittenFrequency +
rcs(FamilySize, 3) + NcountStem + InflectionalEntropy +
Auxiliary + Valency + NVratio + WrittenSpokenRatio,
data = regularity, x = TRUE, y = TRUE)
anova(regularity.lrm)
# ---- model validation
validate(regularity.lrm, bw = TRUE, B = 200)
pentrace(regularity.lrm, seq(0, 0.8, by = 0.05))
regularity.lrm.pen = update(regularity.lrm, penalty = 0.6)
regularity.lrm.pen
# ---- a plot of the partial effects
plot(Predict(regularity.lrm.pen))
# predicting regularity with a support vector machine
library(e1071)
regularity$AuxNum = as.numeric(regularity$Auxiliary)
regularity.svm = svm(regularity[, -c(1,8,10)], regularity$Regularity, cross=10)
summary(regularity.svm)
# }
Run the code above in your browser using DataLab