无法让coefplot()在彼此之上绘制两个回归

时间:2011-07-01 19:09:16

标签: r

我正在尝试使用coefplot(),如本文所述:http://www.r-bloggers.com/visualization-of-regression-coefficients-in-r/

然而,当我运行那些确切的代码时,我只得到一个回归图,而不是3.这是一个屏幕截图,显示我运行的确切代码,以及输出图。 http://i.imgur.com/YtDND.png

我真的不确定还能做什么。非常感谢您的帮助。

2 个答案:

答案 0 :(得分:2)

coefplot的代码不再接受偏移量的参数。它不在文档中,也不在正式列表中。您可以通过修改coefplot的代码来创建一个版本:

键入coefplot2 return 。将函数复制粘贴到命令行并在其前面加上

coefplot2 <-  # the rest of the pasted function should follow

然后将voffset=0添加到formals列表并更改此行:

arrows(ci1, (1:k), ci2, (1:k), lty = lty[1], lwd = lwd[1], col = col,

到此

arrows(ci1, (1:k)+voffset, ci2, (1:k)+voffset, lty = lty[1], lwd = lwd[1], col = col,

并将分数线更改为:

points(cf , (1:k)+voffset, pch = pch, col = col)

然后按回车键,你应该有一个新的coefplot2功能。然后这应该工作

coefplot2(M2, xlim=c(-2, 6) )
par(new=TRUE)    # could not get the add=TRUE argument to work either.
coefplot2(M3, col="blue", xlim=c(-2, 6), voffset=0.4)
par(new=TRUE)
coefplot2(M1, col="red", xlim=c(-2, 6) , voffset=0.2)

答案 1 :(得分:2)

我会指出coefplot2(我的coefplot的扩展,它有一些很好的功能,但仍然需要更多的工作),可以从R-forge获得,这样做:

library(coefplot2)
coefplot2(list(M2,M1,M3),col=c("black","red","blue"),legend=TRUE)