Foro de debate

Agregar datos mensualmente pero teniendo en cuenta el nombre del sensor

 
Imagen de Nuria Naranjo
Agregar datos mensualmente pero teniendo en cuenta el nombre del sensor
de Nuria Naranjo - jueves, 7 de febrero de 2019, 09:22
 

Buenos días, 

Tengo este grupo de datos. La primera columna es el nombre del sensor, las 3 siguientes sus coordenadas, la quinta el dia medido contado desde el 1-1-1975, la sexta el nivel (es la variable de la que quiero obtener una media mensual), la siguiente un contador de datos y la última es la fecha real de medida.

1 1 2 20 105 1 1 16-04-1975

1 1 2 20 173 2 1 23-04-1975

1 1 2 20 201 1 1 25-04-1975

2 2 2 20 232 3 1 21-08-1975

2 2 2 20 271 1 1 29-09-1975

2 2 2 20 300 4 1 30-09-1975

3 3 2 20 313 1 1 10-11-1975

3 3 2 20 338 5 1 05-12-1975

3 3 2 20 388 1 1 24-01-1976

4 4 2 20 419 6 1 24-02-1976

4 4 2 20 447 1 1 24-02-1976

4 4 2 20 500 7 1 24-02-1976

5 5 2 20 502 1 1 17-05-1976

5 5 2 20 532 8 1 16-05-1976

5 5 2 20 574 1 1 28-07-1976


Me gustaría saber la media mensual de las medidas de LEVEL pero que R me de la salida por nombre (columna 1)

data <- read.table("C:/Rmapa/Rdatos/prueba.txt", stringsAsFactors=FALSE)

colnames(data)<- c("name","x","y","z","time","level","contador","fecha")

head

str(data)

name <- data[,1]

date <- data[,5]

level <- as.numeric(as.character(data[,6]))

date_format <- as.Date(date,origin="1975-01-01")

install.packages("zoo")

library(xts)

ts <- xts(data$level, as.Date(data$fecha,"%d-%m-%Y"))

ts_m = apply.monthly(ts,FUN=mean)

time(ts_m)<-as.yearmon(time(ts_m))

ts_m


He aplicado este script pero me da como salida la media mensual pero de todos los datos en general, es decir la media mensual de todo enero de 1975 para todos los sensores, y yo necesito que sea para el sensor 1, la media del nivel es X en enero de 1975 y asi sucesivamente para los meses en los que existan medidas, ya que la serie temporal real es irregular y tiene gaps en algunos meses



Preguntas:

1. ¿Como puedo conseguir que me lo haga mensual para cada sensor ?  xts de zoo es el paquete adecuado?

2. ¿Como puedo importar la salida que obtenga a un txt en el que conserve todos los datos asociados al sensor (name, x, y, z, time(month-year), media de las medidas mensual)?


Muchisimas gracias de antemano