我想使用BMA
包来实现样本预测。
这是我的代码:
out.3m<-bic.glm(x3,y3vec,strict = FALSE,OR = 50,glm.family=binomial(link = "probit"))
yfitout3=predict(out.3m,newdata=tail(x,1)) #where x has one more observation than x3.
我得到了:
predict.bic.glm中的错误(out.3m,newdata = tail(x,1)):
newdata缺少变量
然而,当我使用
时yfitout3=predict(out.3m,newdata=tail(x,2))
有效。
你能告诉我哪里出错了吗?