我想用二项分布族进行非线性回归,其中包含二次项,主要有两个问题:
1.formula该如何表达呢,我这样写
> y1<-Xp[,5]
> x1<-Xp[,2]
> x2<-Xp[,3]
> x3<-Xp[,4]
> lm.form1<-y1~(x1+x2+x3)^2
> lm.solve1<-glm(lm.form1,family="binomial")
结果里面不显示二次项,是不是只能这么写
> lm.form1<-y1~(x1+x2+x3)^2+I(x1^2)+I(x2^2)+I(x3^2)1.这样感觉很麻烦的说。
2.虽然按照上面一行的表达式可以求解,但是总是出现:
警告信息:
In eval(expr, envir, enclos) : 二项通用线性模型里出现了非整数的#成功!
这是否说明拟合不成功呢?
初学R,不知道该如何查询错误列表。请多指教