[R] modeling outliers in arfima

Koopmans, Matthijs mkoopmans at mercy.edu
Fri May 24 21:48:14 CEST 2013


I am analyzing a time series trajectory with many outliers.

When I model the outliers using arima (TSA library),  the following program statement runs like a charm:


'Rate' is the outcome variable. The program generates the regular arima estimates as well as those associated with the outliers (here t=74 and 107).

The corresponding arfima statement (arfima library) follows:

m495.m4a<-arfima(Rate,order=c(1,0,1),lmodel="g", numeach=c(1,1),
 seasonal=list(order=c(1,0,1), period=5,numeach=c(1,1),

The program accepts this language and produces the regular arfima estimates, but not those associated with the outliers.

I' be greatly helped by any insights as to why this might be the case.


Matt Koopmans

More information about the R-help mailing list