Tengo un conjunto de datos que más o menos son así (nótese que esto es un ejemplo mínimo y que lo limitado de los datos puede afectar a la apariencia visual del gráfico obtenido):
y x z g
1 0 0 1
2 1 0 1
2 0 0.5 1
3 1 0.5 1
1.5 0 1 1
2 1 1 1
2 0 0 2
2 1 0 2
3 0 0.5 2
3 1 0.5 2
0.5 0 1 2
2 1 1 2
2 0 0 3
2 1 0 3
1 0 0.5 3
1 1 0.5 3
0.5 0 1 3
0.5 1 1 3
Me gustaría representar en un gráfico la media de y para cada posible combinación de x & z. Representando y en el eje vertical y g en el eje horizontal.
Hasta ahora he utilizaro el siguiente código:
means <- tapply(y,g,mean)
plot(means, col="red",pch=18, ylim=c(0,3), type = 'l', ylab='y', xlab="g")
Después, para cada conjunto de datos (para cada posible combinación de x y z que realizo manualmente con subset
), dibujo una nueva línea en el gráfico, con un color diferente. Uso este código:
lines(means, col="black",pch=18)
Me gustaría poder realizar el gráfico de una manera menos engorrosa, utilizando ggplot. También me gstaría implementar los intervalos de confianza del 95%.
Muchas gracias.