我正在尝试使用coefplot(),如本文所述:http://www.r-bloggers.com/visualization-of-regression-coefficients-in-r/
然而,当我运行那些确切的代码时,我只得到一个回归图,而不是3.这是一个屏幕截图,显示我运行的确切代码,以及输出图。 http://i.imgur.com/YtDND.png
我真的不确定还能做什么。非常感谢您的帮助。
答案 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)