在 r 中的“效果”包的序数回归图中省略截距线

问题描述

我想可视化一个混合效应序数回归,使用 r 中 ordinal 包中的 clmm 函数进行拟合。我正在使用 effects 包。但是,我不希望截距显示在图上(响应类别之间的阈值 - 如下面的水平虚线)。关于如何省略它们的任何指导?

示例代码和绘图:

library(ordinal)
library(MASS)
library(effects)

mod <- clmm(poverty ~ age + (1|country),data = WVS)
plot(Effect("age",mod,latent = TRUE))

enter image description here

解决方法

它读取存储在您的 eff 对象中的阈值。plot.eff 中可能还有其他一些选项,但现在您可以编写一个函数将阈值设置为 NULL 和情节:

plot_without_intercept = function(eff){
       eff$thresholds = NULL
       plot(eff)
}

Eff = Effect("age",mod,latent = TRUE)
plot_without_intercept(Eff)

enter image description here