[R-es] Conversión de objeto temporal (TS) a matriz (o data.frame)

Hola:

Ante todo, gracias por las respuestas.

El 03/06/13 01:23, Oscar Perpiñán Lamigueiro escribió:

Hola,

Si no te entiendo mal, necesitas la conversión a matriz como paso
intermedio para agregar por meses (o cualquier unidad temporal).

¡Hum! No estoy seguro de que agregar sea lo que quiero hacer.

Si es
así, te aconsejo que uses el paquete zoo y el método aggregate para la
clase zoo.

Por ejemplo:

library(zoo)
dats <- data.frame(date=seq(as.Date('2012-01-01'), by='day', length=3*365),
                         A=rnorm(3*365), B=rnorm(3*365), C=rnorm(3*365))

z <- zoo(dats[,-1], dats$date)

Zyearmon <- aggregate(z, by=as.yearmon)

month <- function(x) as.numeric(format(x, "%m"))
year <- function(x) as.numeric(format(x, "%Y"))

Zmon <- aggregate(z, by=month)
Zyear <- aggregate(z, by=year)

Posiblemente sea por desconocimiento de la clase zoo pero... no veo el 
camino para realizar las distintas operaciones que estoy realizando, por 
ejemplo:

# Creación de la matriz de salida.
   matriz.salida<-data.frame(
     ndatos= apply(x.matriz,2,function(x) length(which(!is.na(x)))),
     vminimo= formatC(
       apply(x.matriz,2,function(x) min(x,na.rm=T)),
       dig=2,format="f"),
     vmax=formatC(
       apply(x.matriz,2,function(x) max(x,na.rm=T)),
       dig=2,format="f"),
     media=formatC(
       apply(x.matriz,2,function(x) mean(x,na.rm=T)),
       dig=2,format="f"),
     mediana=formatC(
       apply(x.matriz,2,function(x) median(x,na.rm=T)),
       dig=2,format="f"),
     desvtip=formatC(
       apply(x.matriz,2,function(x) sd(x,na.rm=T)),
       dig=2,format="f"),
     varianza=formatC(
       apply(x.matriz,2,function(x) var(x,na.rm=T)),
       dig=2,format="f"),
     shapiro.valor=formatC(
       aux.shapiro[1,],dig=4,format="f"),
     shapiro.valp=formatC(
       aux.shapiro[2,],dig=4,format="f")
     )

Siendo x.matriz un objeto TS convertido a matriz. Aux.shapiro la calculo 
así:

# Prueba de Shapiro y Wilk
   aux.shapiro<-matrix(unlist(apply(x.matriz,2,function(x)
     shapiro.test(x)[1:2])),
     ncol=12,nrow=2)

¡Ah, vale!:

aggregate(z,by=month,FUN=shapiro.test)

(Aunque esto me fuerza a tirar a la basura la función que había 
construido -con mi buen rato que me había llevado- y hacerla de nuevo)

Última pregunta, con un objeto ts es muy fácil llevarlo a un documento 
como una tabla: xtable(objeto.ts), ¿y con un objeto zoo? (Estoy buscando 
y nada)

Y otra, ¿convertir objetos ts a zoo y viceversa? Estoy probando con 
as.ts() y agua también, pongo unos datos de ejemplo que los que has 
puesto son en exceso numerosos:

lere<-data.frame(date=seq(as.Date("2013-01-01"),by="month",length=3*12),A=rnorm(3*12))
lara<-zoo(lere[,-1],lere$date)
as.ts(lara)

El resultado no es ni remotamente parecido a lo que debería salir.

Saludos.

Oscar

Gracias por todo.

Salud y Revolución.

Lobo.

Pd: Acabo de encontrar como resolver el tema de xtable:

xtable(aggregate(lara,as.yearmon))

Lo curioso es que xtable convierte el nombre de los meses del perfecto 
español en que lo muestra aggregate a la lengua de la pérfida Albión. 
Habrá que pelear con eso.

Libertad es poder elegir en cualquier momento. Ahora yo elijo GNU/Linux,
para no atar mis manos con las cadenas del soft propietario.
Porque la libertad no es tu derecho, es tu responsabilidad.
http://www.mucharuina.com
---------
Desde El Ejido, en Almería, usuario registrado Linux #294013
http://www.counter.li.org