本人在用R软件进行逻辑回归的时候,发现R软件跑出的结果与文献上跑出的结果不一致,求解~~~~左图是文献的结果,右图是自己跑出的结果。右图为啥会多234这种数字呢?????代码如下:setwd("i:/xw/1222")
w=read.table("12221.txt",head=T)
x1=ncol(w)
for(i in 1:x1)w[,i]=factor(w[,i])
a=glm(收视终端~.,w,family="binomial")
z=(predict(a,w,type="response")>0.5)
u=rep("NO",956);u[!z]="AB" #把预测结果转换成原先的值(AB或NO)
zz=table(w[,11],u) #2乘2矩阵,对角线外的数目为分错的数目
(sum(zz)-sum(diag(zz)))/sum(zz) #计算错误率
zz;summary(a)