在为 R 中的置信区间编码时如何使用 dplyr 获取列的长度答案

【问题标题】：How to obtain length of column using dplyr when coding for confidence interval in R在为 R 中的置信区间编码时如何使用 dplyr 获取列的长度
【发布时间】：2017-12-12 10:59:44
【问题描述】：

我的数据集可以参考here。总而言之，我有一个名为 fit 的列，我需要为其绘制均值和置信区间。

我正在尝试使用 dplyr 函数为我的 ggplot 编写置信区间

data.melt$time = factor(data.melt$time, levels=paste("t", seq(0, 10), sep=""))

这里是代码

summary_dat = data.melt$time  %>%
              group_by(resource, fertilizer, time) %>%
              summarise(mean_predict=mean(fitted),
                        sd_predict = sd(fitted),
                        n_predict = n(fitted)) %>%

  mutate(se = sd_predict / sqrt(n_predict),
         lower_ci = mean_predict - qt(1 - (0.05 / 2), n_predict - 1) * se_predict,
         upper_ci = mean_predict + qt(1 - (0.05 / 2), n_predict - 1) * se_predict)

但是，R 不允许我将 n_predict 编码为 n(fitted)。我也试过长度（合身）但没有运气。有什么想法吗？

【问题讨论】：

如果您使用length() 而不是n()，您的方法应该有效，但是您的代码有很多小错别字。第一行应该以整个数据框data.melt 开头，而不是data.melt$time，有问题的列称为predict，而不是fitted，标准错误由两个不同的名称引用，se 和se_predict.

标签： r dplyr confidence-interval sample-size

【解决方案1】：

dplyr 中的便捷函数 n() 仅计算子集中的行数，而不计算子集中的列之一的长度。它不需要争论。您想使用n_predict = n() 或n_predict = length(predict)。

【讨论】：