问题描述
我想可视化一个混合效应序数回归,使用 r 中 ordinal
包中的 clmm 函数进行拟合。我正在使用 effects
包。但是,我不希望截距显示在图上(响应类别之间的阈值 - 如下面的水平虚线)。关于如何省略它们的任何指导?
示例代码和绘图:
library(ordinal)
library(MASS)
library(effects)
mod <- clmm(poverty ~ age + (1|country),data = WVS)
plot(Effect("age",mod,latent = TRUE))
解决方法
它读取存储在您的 eff
对象中的阈值。plot.eff
中可能还有其他一些选项,但现在您可以编写一个函数将阈值设置为 NULL
和情节:
plot_without_intercept = function(eff){
eff$thresholds = NULL
plot(eff)
}
Eff = Effect("age",mod,latent = TRUE)
plot_without_intercept(Eff)