[R-es] Análisis de una serie temporal diaria
Paco Pastor
paco en ceam.es
Mar Abr 12 15:30:39 CEST 2011
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>
Más información sobre la lista de distribución R-help-es