在
数据挖掘与R语言海藻预测一章中有如下代码:cv.rpart <- function(form,train,test,...) {
m <- rpartXse(form,train,...)
p <- predict(m,test)
mse <- mean((p-resp(form,test))^2)
c(nmse=mse/mean((mean(resp(form,train))-resp(form,test))^2))
}
cv.lm <- function(form,train,test,...) {
m <- lm(form,train,...)
p <- predict(m,test)
p <- ifelse(p < 0,0,p)
mse <- mean((p-resp(form,test))^2)
c(nmse=mse/mean((mean(resp(form,train))-resp(form,test))^2))
}
res <- experimentalComparison(
c(dataset(a1 ~ .,clean.algae[,1:12],'a1')),
c(variants('cv.lm'),
variants('cv.rpart',se=c(0,0.5,1))),
cvSettings(3,10,1234))
请教下experimentalComparison()函数是怎么调用上面function自编函数的?
function自编函数中form,train,test分别由哪个数据传输入的?还有为什么要用resp()函数,而不直接用test值,不就是用预测值p和test值做比较么?
哪位前辈帮忙解答下,万分感谢。