我这样做了:
reg1_v1 <- multinom(Y ~ X1 + X2, data = y1, model = TRUE)
r1_v1 = summary(reg1_v1)
mod = r1_v1$model
fit1_v1 = fitted(reg1_v1)
y1_predict = cbind(mod, fit1_v1)
所以我所做的是尝试将模型中使用的原始数据与拟合概率相结合。 fit1_v1中有3个拟合概率列,称为prob1,prob2,prob3。然后我尝试使用这些概率进行一些计算。然而,不能使用prob3 - 它在我加载数据集并查看数字时显示,但我不能通过它的列标题“prob3”来调用它。这就像我的数据中的“幽灵列”。例如,如果我
ratio = y1_predict$prob1/y1_predict$prob3
它会说:
Error in `$<-.data.frame`(`*tmp*`, "ratio", value = numeric(0)) :
replacement has 0 rows, data has 563778
但我可以像往常一样使用prob1和prob2。知道为什么会这样,以及如何解决?谢谢。