我正在尝试使用一组变量预测一个三级因子。这个组别是A、B和C。

m<-glm(as.factor(Group)~Sex+BloodType+Pressure,data=Hel,family = "binomial")newdata <- data.frame(Sex="M",BloodType="A+", Pressure=80)predict(m,newdata)

这返回的结果是：

        1 0.7133324 

我希望得到的结果是这样的：

A        B        C20.00   40.00    40.000

我该怎么做呢？谢谢。

回答：

我想你得使用支持多项式回归的工具，比如nnet？

library(nnet)trn = sample(1:nrow(iris),100)fit_nnet <- multinom(Species ~ ., iris[trn,])head(predict(fit_nnet,iris[-trn,],type="prob"))      setosa   versicolor    virginica7  1.0000000 2.348893e-11 2.166176e-589  0.9999323 6.765017e-05 2.438304e-5210 0.9999940 6.035319e-06 1.530614e-5612 1.0000000 1.355834e-09 1.315066e-57