使用来自R的coefplot包的命令多值来面对错误

时间:2016-01-07 21:34:31

标签: r coefplot

multiplotR运行命令coefplot时,有没有人像我一样?

即使是例子:

data(diamonds)
model1 <- lm(price ~ carat + cut, data=diamonds)
model2 <- lm(price ~ carat + cut + color, data=diamonds)
model3 <- lm(price ~ carat + color, data=diamonds)
multiplot(model1, model2, model3)

我现在收到以下错误:

Error in get(x, envir = this, inherits = inh)(this, ...) : 
attempt to apply non-function

任何提示?

1 个答案:

答案 0 :(得分:4)

这个答案有点切,但我发现broomdotwhisker包的稍微更近的组合是有用的 - broom是后端(将模型转换为&#34; tidy&#34;系数的数据帧)和dotwhisker是前端(通过相当薄的ggplot2层创建绘图)

library(ggplot2)
library(dotwhisker)
library(broom)

更新:重新安装dotwhisker v 0.2.0.3,这似乎有效:

dwplot(list(model1,model2,model3))

如果您想动态选择不同的型号名称,也可以dwplot(list(m1=model1,m2=model2,m3=model3))

或者,为了更好地控制,您可以自己构建完整的数据框:

mList <- list(carat_cut=model1, carat_cut_color=model2,
              carat_color=model3)
library(plyr)
## extract tidy data frames and combine them ...
mFrame <- ldply(mList,tidy,conf.int=TRUE,.id="model")

现在你可以做到

dwplot(mFrame)