我正在尝试在 R 中绘制延续比率模型的部分残差。
我正在使用以下代码：
load(url("http://www.produnis.de/R/OrdinalSample.RData"))
mydata <- ordinalSample
head(mydata)
require(VGAM)
cr <- vglm(Mood ~ JobFamily + JobSatisfaction, data=mydata, family=sratio(parallel=T,reverse=F));cr
par(mfrow = c(2, 1))
plotvgam(cr,residuals=T)
par(mfrow = c(1, 1))
dev.off()
不幸的是，这给了我以下错误消息：
Warning:
1: In vplot.numeric(x, ...) :
ncol(residuals) != ncol(y) so residuals are not plotted
2: In vplot.numeric(x, ...) :
ncol(residuals) != ncol(y) so residuals are not plotted
我不明白为什么会发生这种情况，我很感谢有关如何解决它的任何建议......
rms我知道如何使用带有命令的“”包绘制部分残差图residuals.lrm(fit,type="partial",pl=T)，不幸的是，此命令不适用于 vglm-fitted 模型......