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 <[email protected]> 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:[email protected]] En nombre de david > santolaria > Enviado el: miércoles, 23 de marzo de 2016 11:16 > Para: r-help-es <[email protected]> > 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 > [email protected] > https://stat.ethz.ch/mailman/listinfo/r-help-es > >
_______________________________________________ R-help-es mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-help-es
