我试图在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拟合模型...