Estoy tratando datos climáticos de precipitación en R y tengo una serie de datos donde hay un valor de lluvia por cada día desde 1964.
En algunas ocasiones existen datos faltantes que los transformo en NA al importar el excel. Hay una columna (día) con el día del año, otra (mes) con el numero del mes, otra (año, o ano) con el numero de año.
Necesito generar la media de cada mes para todos los años. Por ejemplo la media de precipitación de todos los Julios o de todos los Eneros...
Estoy usando este código, pero no parece ser el adecuado:
if(pira_cacchan_recu_99_9_p_1_2_1964$mes == 1){is.na(media_pira <- mean(pira_cacchan_recu_99_9_p_1_2_1964$ppira),
media_cacchan <- mean(pira_cacchan_recu_99_9_p_1_2_1964$pcacchan),
media_recuay <- mean(pira_cacchan_recu_99_9_p_1_2_1964$precuay)}
print(media_pira)
print(media_cacchan)
print(media_recuay)
Alguna idea sobre qué puedo hacer? o como llevar la media a cabo?
Para poneros en situación:
- pira_cacchan_recu_99_9_p_1_2_1964 es el excel
- ppira, cacchan y recuay, son columnas en el excel que corresponden con cada estación meteoroógica.
- la columna pira_cacchan_recu_99_9_p_1_2_1964$mes se corresponde con el numero del mes a tratar. Como si fuese un código o identificador.