[R-es] ACF retardos múltiplos del periodo

[Isidro Hidalgo Arellano]

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