使用 R 中的 GGplot2 绘制带有误差线的趋势图

Question

我有以下数据，是按集团和季度汇总的。我想为每个组绘制三个不同的趋势图，并带有误差线。数据已经汇总，我希望按原样绘制值。我对 ggplot2 相当陌生，正在努力计算绘图值周围的误差线。有人可以帮忙吗？

Group     Quarter    Proportion
1          1         1.03
1          2         4.02
2          1         2.02
2          2         6.01
3          1         2.01
3          2         8.09

Answer 1

数据显示比例。这非常重要，因为我们无法计算与均值相同的比例置信区间。您可以在此处或此处找到有关此主题的小插图。

幸运的是，所有这些都已在名为

ANOPA

的包中实现（使用反正弦变换进行比例分析）。

首先，要拥有一个数据集，让我们使用

GRP

生成成功 (1) 或失败 (0) 数据（生成随机比例）：

library(ANOPA)
BSDesign <- list(Group = c(1,2,3), Quarter = c(1,2) )
thePs    <- c(.01, .04, .02, .06, .02, .08)
res      <- GRP( thePs, 6000, BSDesign )

人口比例

thePs

与问题中给出的相符；根据问题，

BSDesign

中描述的设计是 3 x 2 设计，有 6 组。我生成了一个非常大的样本（每组 n = 6000），因为问题中考虑的比例非常小。

从那里，您可以使用

分析数据

w <- anopa(s ~ Quarter + Group, res)
summarize(w)

###                      MS  df         F pvalue correction    Fcorr pvalcorr
### Quarter        0.008614   1 206.74613      0   1.000028 206.7404        0
### Group          0.001362   2  32.68931      0   1.000056  32.6875        0
### Quarter:Group  0.007651   2 183.63635      0   1.000486 183.5471        0
### Error(between) 0.000042 Inf        

毫不奇怪，一切都非常重要（所有ps< .01) because the sample is so large.

最后一点，有一个情节

anopaPlot(w)

我们很容易看到绿线（第 2 组）的交互结果，其比例与其他两组的比例表现完全不同。

该图是一个 ggplot2 对象，因此您可以通过添加图形指令来自定义它，例如，

anopaPlot(w) + theme_bw()

您可以在一对列 {s, n} 中看到编译后的数据

w$compData
###   Group Quarter   s    n
### 1     1       1  49 6000
### 3     2       1 227 6000
### 5     3       1 116 6000
### 2     1       2 351 6000
### 4     2       2  99 6000
### 6     3       2 459 6000