En primer lugar, para el calculo de los intervalos de confianza voy a usar las manera más simple, haciendo una regresión lineal y calculando los intervalos sobre ésta (Ver fuente)
ci <- function(x, level=0.95, bound=1) {
l.model <- lm(x ~ 1, list(x))
as.vector(confint(l.model, level=level))[bound]
}
Leemos tus datos:
data <- read.csv('~/Descargas/data_95.csv', check.names = FALSE)
Lo primero que hago es a) ignorar la columna Mean
, vamos a volver a calcularla, b) expandir los 5 valores de Repl
de forma vertical c) con esta estructura larga, podremos agrupar y sumarisar para obtener los coeficientes buscados de manera sencilla:
data %>%
select(-Mean) %>%
pivot_longer(c(-N)) %>%
group_by(N) %>%
summarise(mean = mean(value),
desvio = sd(value),
error_est = desvio / sqrt(n()),
intervalo_sup = ci(value, bound = 2),
intervalo_inf = ci(value, bound = 1),
.groups="keep") -> plot_data
plot_data
# A tibble: 80 × 6
# Groups: N [80]
N mean desvio error_est intervalo_sup intervalo_inf
<int> <dbl> <dbl> <dbl> <dbl> <dbl>
1 1 0 0 0 0 0
2 2 0 0 0 0 0
3 3 40.8 25.7 11.5 72.7 8.91
4 4 208. 23.9 10.7 238. 179.
5 5 195 82.9 37.1 298. 92.0
6 6 181. 102. 45.4 307. 55.1
7 7 156 89.6 40.1 267. 44.8
8 8 149. 84.5 37.8 254. 43.9
9 9 101. 63.5 28.4 180. 22.6
10 10 125. 74.9 33.5 218. 31.8
# … with 70 more rows
Por último, puedes graficar la media y los intervalos haciendo algo así:
plot_data %>%
ggplot(aes(x=N, y=mean)) +
geom_errorbar(aes(ymax = intervalo_sup, ymin = intervalo_inf)) +
geom_point(colour="blue")
Reultado:
