Hello,

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:

m3.m495<-arima(Rate,order=c(1,0,1),seasonal=list(order=c(1,0,1),period=5),
  xreg=data.frame(t74=1*(seq(Rate)==74),t107=1*(seq(Rate)==107)))

'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),
 xreg=data.frame(t74=1*(seq(Rate)==74),t107=1*(seq(Rate)==107)))

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.

Regards,

Matt Koopmans

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to