在进行序数逻辑回归时获得概率的置信区间

Question

我有一个包含以下列的数据框（培训）：

       GPAGrpd    Subj1.Grade
1       [0,2)          6
2       [0,2)          6
3       [3,4.3]        1
4       [0,2)          6
5       [2,3)          2
6       [3,3.7)        3

我运行了以下命令：

 Res.clm <- clm(GPAGrpd ~ Subj1.Grade,data=training)
 summary(Res.clm)
 newdat <- data.frame(
      Subj1.Grade = rep(1:6, each = 200),
      GPAGrpd = factor(rep(1:4, each = 300))
 )
 newdat <- cbind(newdat, predict(Res.clm, newdata=newdat, se.fit=TRUE,
          interval=TRUE, type="prob"))

但是，我收到一条错误消息，指出：Error in predict.clm(Res.clm, newdata = newdat, se.fit = TRUE, interval = TRUE, : response factor 'GPAGrpd' has new levels

Answer 1

我能够通过定义newdat来解决我的错误，如下所示：

    newdat <- data.frame(
      Subj1.Grade = rep(1:6, 200),
      GPAGrpd = factor(rep(1:4, each=300), 
                   labels = c("[0,2)","[2,3)","[3,3.7)","[3.7,4.3]"))
    )

事实证明，当我试图绘制预测电话的结果时，我确实需要GPAGrpd。谢谢你的帮助。