Мой набор данных можно сослаться на здесь. В сущности, у меня есть столбец, называемый приспособленным, для которого мне нужно построить средние и доверительные интервалы.
Я пытаюсь закодировать доверительный интервал, используя функцию 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(подходящий). Я также пробовал длину (установлено), но не повезло. Любые идеи?