Hola, Estoy visualizando una serie temporal para determinar sus órdenes ARIMA y no consigo lo siguiente: ¿Cómo puedo sacar la ACF de los retardos múltiplos del periodo? Es decir, sólo ver en el gráfico ACF los retardos 12, 24, 36... Gracias!! David [[alternative HTML version deleted]]
Puedes extraer los acf de la "serie" (los tienes en serie$acf para los retardos serie$lag), los seleccionas con seq() y los pones en un gráfico... Isidro Hidalgo Arellano Observatorio del Mercado de Trabajo Consejería de Economía, Empresas y Empleo http://www.castillalamancha.es/ -----Mensaje original----- De: R-help-es [mailto:r-help-es-bounces en r-project.org] En nombre de david santolaria Enviado el: miércoles, 23 de marzo de 2016 11:16 Para: r-help-es <r-help-es en r-project.org> Asunto: [R-es] ACF retardos múltiplos del periodo Hola, Estoy visualizando una serie temporal para determinar sus órdenes ARIMA y no consigo lo siguiente: ¿Cómo puedo sacar la ACF de los retardos múltiplos del periodo? Es decir, sólo ver en el gráfico ACF los retardos 12, 24, 36... Gracias!! David [[alternative HTML version deleted]] _______________________________________________ R-help-es mailing list R-help-es en r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es
Por si no me he explicado bien: datos <- c(123152, 124017, 122268, 118182, 114263, 113023, 108888, 108001, 107205, 109355, 110795, 114385 117895, 118650, 116934, 112564, 107998, 106336, 101285, 100865, 100753, 101670, 103811, 110674 110601, 109846, 107855, 103080, 98199, 95838, 90920, 90735, 88838, 90299, 90180, 94923 , 95462, 93472, 94240, 93365, 90809, 87286, 82626, 83058, 83072, 84751, 85931, 89274 , 93055, 93320, 92301, 89590, 86411, 84892, 82579, 82211, 82149, 83838, 85063, 88302 , 93001, 93581, 93645, 88520, 88157, 81964, 79467, 81933, 82111, 83271, 85485, 88037 , 93366, 96485, 95470, 91052, 87301, 85016, 83026, 83756, 84414, 85971, 88297, 92219 , 96183, 96484, 95855, 91768, 88677, 88018, 85578, 86911, 87816, 90414, 92071, 97191 , 99259, 99393, 99994, 99858) serie <- ts(datos, start= c(1996, 1), frequency= 12) acf_serie <- acf(serie, lag.max= 36) seleccionar <- acf_serie$lag == as.integer(acf_serie$lag) plot(acf_serie$lag[seleccionar], acf_serie$acf[seleccionar]) Un saludo Isidro Hidalgo Arellano Observatorio del Mercado de Trabajo Consejería de Economía, Empresas y Empleo http://www.castillalamancha.es/ -----Mensaje original----- De: R-help-es [mailto:r-help-es-bounces en r-project.org] En nombre de david santolaria Enviado el: miércoles, 23 de marzo de 2016 11:16 Para: r-help-es <r-help-es en r-project.org> Asunto: [R-es] ACF retardos múltiplos del periodo Hola, Estoy visualizando una serie temporal para determinar sus órdenes ARIMA y no consigo lo siguiente: ¿Cómo puedo sacar la ACF de los retardos múltiplos del periodo? Es decir, sólo ver en el gráfico ACF los retardos 12, 24, 36... Gracias!! David [[alternative HTML version deleted]] _______________________________________________ R-help-es mailing list R-help-es en r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es
Lo primero de todo, muchas gracias Isidro por tus respuestas. He probado tu segunda solución y parece correcta, pero pierdo el formato de gráfico ACF. Me refiero a que pierdo la banda de confianza del cero, las lineas... (ver imagenes). ¿Alguna otra solución? Dado que es algo, relativamente común, no se si existe algo implementado en algún paquete. Yo estoy usando "Forecast", pero no si hay algún paquete más completo o mejor. David [image: Imágenes integradas 1] El 23 de marzo de 2016, 12:20, Isidro Hidalgo Arellano <ihidalgo en jccm.es> escribió:> Por si no me he explicado bien: > > datos <- c(123152, 124017, 122268, 118182, 114263, 113023, 108888, 108001, > 107205, 109355, 110795, 114385 > 117895, 118650, 116934, 112564, 107998, 106336, 101285, 100865, 100753, > 101670, 103811, 110674 > 110601, 109846, 107855, 103080, 98199, 95838, 90920, 90735, 88838, 90299, > 90180, 94923 > , 95462, 93472, 94240, 93365, 90809, 87286, 82626, 83058, 83072, 84751, > 85931, 89274 > , 93055, 93320, 92301, 89590, 86411, 84892, 82579, 82211, 82149, 83838, > 85063, 88302 > , 93001, 93581, 93645, 88520, 88157, 81964, 79467, 81933, 82111, 83271, > 85485, 88037 > , 93366, 96485, 95470, 91052, 87301, 85016, 83026, 83756, 84414, 85971, > 88297, 92219 > , 96183, 96484, 95855, 91768, 88677, 88018, 85578, 86911, 87816, 90414, > 92071, 97191 > , 99259, 99393, 99994, 99858) > > serie <- ts(datos, start= c(1996, 1), frequency= 12) > acf_serie <- acf(serie, lag.max= 36) > seleccionar <- acf_serie$lag == as.integer(acf_serie$lag) > plot(acf_serie$lag[seleccionar], acf_serie$acf[seleccionar]) > > Un saludo > > Isidro Hidalgo Arellano > Observatorio del Mercado de Trabajo > Consejería de Economía, Empresas y Empleo > http://www.castillalamancha.es/ > > > > > -----Mensaje original----- > De: R-help-es [mailto:r-help-es-bounces en r-project.org] En nombre de david > santolaria > Enviado el: miércoles, 23 de marzo de 2016 11:16 > Para: r-help-es <r-help-es en r-project.org> > Asunto: [R-es] ACF retardos múltiplos del periodo > > Hola, > > Estoy visualizando una serie temporal para determinar sus órdenes ARIMA y > no > consigo lo siguiente: ¿Cómo puedo sacar la ACF de los retardos múltiplos > del > periodo? Es decir, sólo ver en el gráfico ACF los retardos 12, 24, 36... > > Gracias!! > David > > [[alternative HTML version deleted]] > > _______________________________________________ > R-help-es mailing list > R-help-es en r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es > >------------ próxima parte ------------ Se ha borrado un adjunto en formato HTML... URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20160323/6a7b620a/attachment-0001.html> ------------ próxima parte ------------ A non-text attachment was scrubbed... Name: image.png Type: image/png Size: 22298 bytes Desc: no disponible URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20160323/6a7b620a/attachment-0001.png>