在multiplot
包R
运行命令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
任何提示?
答案 0 :(得分:4)
这个答案有点切,但我发现broom
和dotwhisker
包的稍微更近的组合是有用的 - 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)