[R-es] Análisis de una serie temporal diaria

[Paco Pastor]

Hola a todos

Aunque parezca mentira sigo enzarzado con esto de las series diarias. 
Además de hallar una regresión lineal simple quiero utilizar stl o 
decompose para encontrar la estacionalidad y la tendencia de una serie 
diaria de temperatura pero no hay manera. En la creación del objeto ts 
(desde uno zoo) de la serie diaria creo que estoy haciendo algo mal ya 
que no obtengo resultados. El motivo de querer untilizar stl es que ya 
lo había empleado para una serie de datos mensuales.

Si intento leer directamente como ts

temperatura=read.csv("tmax.dat",header=F,sep=",")

fecha=as.Date(temperatura[,1],"%d/%m/%Y")

temp=temperatura[order(fecha),]

media=rowMeans(temp[2:length(temp)], na.rm=TRUE)

temp=cbind(temp,media)

tsdata=ts(temp$media,freq=7,start=c(1968,1))

a=decompose(tsdata)

plot(a)


se obtiene la gráfica que os adjunto y que creo que no es correcta. Mi 
duda es si al poner freq=t lo que estoy buscando es la componente 
estacional en base periodos semanales. Tengo claro que para datos 
mensuales la frecuencia es 12 ¿pero cual es para datos diarios de 40 
años? ¿365? En la ayuda de ts dice que para una serie de datos diarios 
con periodo natural semanal freq=7. Pero yo no quiero analizar una serie 
de datos semanales sino diarios, no tengo muy claro lo que estoy haciendo.

Si primero creo un objeto zoo

library(zoo)

temperatura=read.csv("tmax.dat",header=F,sep=",")

fecha=as.Date(temperatura[,1],"%d/%m/%Y")

temp=temperatura[order(fecha),]

media=rowMeans(temp[2:length(temp)], na.rm=TRUE)

zoodata=zoo(temp$media,fecha)

zz=as.ts(zoodata,frequency=7,start=c(1968,1))     # los datos comienzan el 1 de enero de 1968

decompose(zz)

Error in decompose(zz) : time series has no or less than 2 periods


Entonces aparece el "famoso" error de que la serie no es periódica.

Estoy ya un poco desesperado, ¿existe alguna posibilidad por esta 
vía?¿hay otras funciones que puedan analizar los datos?¿alguien tiene un 
ejemplo que pueda adaptar a mis datos? No siendo muy ducho en 
estadística quería (necesitaba) solamente algún tipo de análisis 
relativamente simple o no muy elaborado, algo que me dé una idea de la 
evolución de los datos.

Gracias de nuevo a todos

Paco

-----------

Francisco Pastor
Meteorology department, Instituto Universitario CEAM-UMH
http://www.ceam.es
-----------
mail: paco en ceam.es
skype: paco.pastor.guzman
Researcher ID: http://www.researcherid.com/rid/B-8331-2008
Cosis profile: http://www.cosis.net/profile/francisco.pastor
-----------
Parque Tecnologico, C/ Charles R. Darwin, 14
46980 PATERNA (Valencia), Spain
Tlf. 96 131 82 27 - Fax. 96 131 81 90


---------------------------------------------------------------------
Este mensaje y los ficheros anexos son confidenciales. Los mismos contienen información reservada de la empresa que no puede ser difundida. Si usted ha recibido este correo por error, tenga la amabilidad de eliminarlo de su sistema y avisar al remitente mediante reenvío a su dirección electrónica; no deberá copiar el mensaje ni divulgar su contenido a ninguna persona.

Su dirección de correo electrónico junto a sus datos personales forman parte de un fichero titularidad de la Fundación de la Comunidad Valenciana Centro de Estudios Ambientales del Mediterráneo - CEAM, con CIF: G-46957213, cuya finalidad es la de mantener el contacto con Ud. De acuerdo con la Ley Orgánica 15/1999, usted puede ejercitar sus derechos de acceso, rectificación, cancelación y, en su caso, oposición enviando una solicitud por escrito, acompañada de una fotocopia de su DNI dirigida a: Fundación de la Comunidad Valenciana Centro de Estudios Ambientales del Mediterráneo - CEAM. C/ Charles R. Darwin, 14. Parque Tecnológico.46980 PATERNA (Valencia).

This message and the attached files are confidential. They contain reserved information belonging to our centre and are not to be broadcast. If you have received this email by mistake, please delete it from your system and alert the sender by returning it to his/her email address. You must not copy or divulge the contents of the message to anyone.

Your email address and personal data are included in a file belonging to the Fundación de la Comunidad Valenciana Centro de Estudios Ambientales del Mediterráneo - CEAM, con CIF: G-46957213. The purpose of this file is to allow us to keep in contact with you. In accordance with Organic Law 15/1999, you are permitted to access, rectify, cancel or oppose the contents of this file by submitting a written request, accompanied by a photocopy of your DNI, to: Fundación de la Comunidad Valenciana Centro de Estudios Ambientales del Mediterráneo - CEAM. C/ Charles R. Darwin, 14. Parque Tecnológico.46980 PATERNA (Valencia).





------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: plot_decompose_tsdata.png
Type: image/png
Size: 56362 bytes
Desc: no disponible
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20110412/c1fff41d/attachment-0001.png>