在编码R中的置信区间时，如何使用DPLYR获得列的长度

Question

我的数据集可以在此处引用。在GIST中，我有一个名为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（拟合）。我还尝试了长度（安装），但没有运气。有什么想法吗？

Answer 1

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