我正在尝试从R中的Stata复制图表。我有几个变量,并希望在每个治疗组中显示其平均值,其中有两个。 Stata图形如下:
该系数图实际上不是系数图,而是每个单独变量的每种处理方法的平均值。 df基本上看起来像。
答案 0 :(得分:1)
如果没有可重复的数据,很难回答您的问题。
但是,这可能会满足您的期望:
library(dplyr)
mpg %>%
select(manufacturer, cty, trans) %>%
group_by(manufacturer, trans) %>%
summarize(cty_mean = mean(cty)) %>%
ggplot(aes(x=cty_mean, y=reorder(manufacturer, cty_mean), color=trans)) +
geom_point()
如果您还希望包括系数或标准误差,则可以通过在summarize()
中包含一个函数来实现。
答案 1 :(得分:0)