吴老师的《复杂数据统计方法——基于R的应用》中§3.2的数据,使用bagging分类原代码
> set.seed(1044)
> a=bagging(NSP~.,w)
> z0=table(w[,23],predict(a,w)$class)
错误于predict(a, w)$class : $ operator is invalid for atomic vectors
此外: 共有25个警告 (用warnings()来显示)
> z0;(E0=(sum(z0)-sum(diag(z0)))/sum(z0))
错误: 找不到对象'z0'
> barplot(a$importance)
错误于barplot.default(a$importance) : 'height'要么是向量,要么是矩阵
> set.seed(1004)
> E0=rep(0,5);E1=E0
> for(i in 1:5){m=mm[[i]]
+ n0=2126-length(m);n1=length(m)
+ a=bagging(NSP~.,w[-m])
+ E0[i]=sum(w[-m,23]!=predict(a,w[-m,])$class)/n0
+ E1[i]=sum(w[m,23]!=predict(a,w[m,])$class)/n1}
错误于`[.default`(w, -m) : 只有负下标里才能有零
> mean(E0);mean(E1)
[1] 0
[1] 0
数据已经删改过了,还要在完全不懂的情况下要给别人讲QAQ,可肿么办。。。