[R] ARIMA standard error

Andrew Robinson A.Robinson at ms.unimelb.edu.au
Fri Mar 16 04:06:05 CET 2007


Hi Gad,

try:


> class(a)
[1] "Arima"
> getAnywhere(print.Arima)

...

Cheers,

Andrew


On Fri, Mar 16, 2007 at 01:47:25PM +1100, Gad Abraham wrote:
> Hi,
> 
> Can anyone explain how the standard error in arima() is calculated?
> 
> Also, how can I extract it from the Arima object? I don't see it in there.
> 
>  > x <- rnorm(1000)
>  > a <- arima(x, order = c(4, 0, 0))
>  > a
> 
> Call:
> arima(x = x, order = c(4, 0, 0))
> 
> Coefficients:
>            ar1     ar2     ar3      ar4  intercept
>        -0.0451  0.0448  0.0139  -0.0688     0.0010
> s.e.   0.0316  0.0316  0.0317   0.0316     0.0296
> 
> sigma^2 estimated as 0.9775:  log likelihood = -1407.56,  aic = 2827.12
>  > names(a)
>   [1] "coef"      "sigma2"    "var.coef"  "mask"      "loglik"    "aic" 
>        "arma"      "residuals" "call"      "series"    "code" 
> "n.cond"    "model"
> 
> 
> Thanks,
> Gad
> 
> -- 
> Gad Abraham
> Department of Mathematics and Statistics
> The University of Melbourne
> Parkville 3010, Victoria, Australia
> email: g.abraham at ms.unimelb.edu.au
> web: http://www.ms.unimelb.edu.au/~gabraham
> 
> ______________________________________________
> R-help at stat.math.ethz.ch 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.

-- 
Andrew Robinson  
Department of Mathematics and Statistics            Tel: +61-3-8344-9763
University of Melbourne, VIC 3010 Australia         Fax: +61-3-8344-4599
http://www.ms.unimelb.edu.au/~andrewpr
http://blogs.mbs.edu/fishing-in-the-bay/



More information about the R-help mailing list