结果〜1 +速度+ Reltailsize + tailcontrst + relsvl:速度+ Reltailsize:速度+ Reltailsize:relsvl +策略:relsvl + tailcontrst:Reltailsize
我用它装了它 需要(MASS) polr(结果~1 + Speed + Relailsize + tailcontrst + relsvl:Speed + Relailsize:Speed + Relailsize:relsvl + strategy:relsvl + tailcontrst:Relailsize,data = results)
得到系数: 我现在的问题是我现在该如何重新计算“结果”的值来代替因子的值?
我对此有两个具体的疑问, 1)我用什么代替“:”?也许 ”/” ? 2)我是否需要重新使用probit链接功能?我该怎么做?
答案 0 :(得分:0)
我提出了我的问题。 :不能被/代替,因为它意味着因子向量应该逐个倍增(Venables& Ripley,2002)。 仍然不知道如何使用值1,2和3进行有序的概率转换。 干杯 Agus的
答案 1 :(得分:0)
我提前了。现在我使用来自nnet包的multinom来拟合模型,并使用predict.nnet从具有实际数据的数据帧中获取预测值。
M=multinom(result~1+Speed+relsvl, data=results)
summary(M)
require(nnet)
d=predict (M, realdata)
干杯