data(GlaucomaMVF)
response <- function (data) {
attach(data)
res <- ifelse((!is.na(clv) & !is.na(lora) & clv >= 5.1 & lora >=
49.23372) | (!is.na(clv) & !is.na(lora) & !is.na(cs) &
clv < 5.1 & lora >= 58.55409 & cs < 1.405) | (is.na(clv) &
!is.na(lora) & !is.na(cs) & lora >= 58.55409 & cs < 1.405) |
(!is.na(clv) & is.na(lora) & cs < 1.405), 0, 1)
detach(data)
factor (res, labels = c("glaucoma", "normal"))
}
mypredict.inclass <- function(object, newdata){
res <- predict.inclass(object = object, cFUN = response, newdata = newdata)
return(res)
}
errorest(clv+lora+cs~., data = GlaucomaMVF, model=inclass,
predict=mypredict.inclass, iclass="Class", estimator="cv",
pFUN = rpart)
Run the code above in your browser using DataLab