根据帮助multinom
, 包裹nnet
,“响应应该是一个因子或具有 K 列的矩阵,它将被解释为每个 K 类的计数。”我尝试在第二种情况下使用此函数,但出现错误。
这是我所做的示例代码:
response <- matrix(round(runif(200,0,1)*100),ncol=20) # 10x20 matrix of counts
predictor <- runif(10,0,1)
fit1 <- multinom(response ~ predictor)
weights1 <- predict(fit1, newdata = 0.5, "probs")
这是我得到的:
'newdata' had 1 row but variables found have 10 rows
我怎么解决这个问题?
额外问题:我还注意到我们可以将多项式与因子预测器一起使用,例如predictor <- factor(c(1,2,2,3,1,2,3,3,1,2))
。我无法理解这在数学上是如何可能的,因为多项式线性 Logit 回归应该仅适用于连续或二分预测变量。