有人可以解释为什么predict
无效speedglm
?
工作代码:
library(speedglm)
mtcars2<-mtcars
mtcars2$gear<-as.factor(mtcars2$gear)
mtcars_train<-mtcars2[1:10,]
mtcars_test<-mtcars2[11:nrow(mtcars2),]
model<-speedglm(formula = cyl ~ gear,data = mtcars_train,family=poisson(link="log"))
pred<-predict(object = model, newdata = mtcars_test)
pred
我只得到NA。
是否可以修改glm
对象的类,使其使用常规predict
函数?我的意思是,speedglm
对象与常规glm
对象的形式相同,对吗?
class(model)
#[1] "speedglm" "speedlm"