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

【问题标题】：Getting confidence intervals on probabilities when doing ordinal logistic regression在进行序数逻辑回归时获得概率的置信区间
【发布时间】：2016-07-22 16:32:41
【问题描述】：

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

       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

【问题讨论】：

是的，它似乎奏效了。

标签： r regression logistic-regression ordinal

【解决方案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。感谢您的帮助。

【讨论】：