[R] arma: what is the meaning of Pr(>|t|)?

Prof Brian Ripley ripley at stats.ox.ac.uk
Wed Aug 20 16:24:28 CEST 2008


On Wed, 20 Aug 2008, Prof Brian Ripley wrote:

> On Wed, 20 Aug 2008, Alberto Monteiro wrote:
>
>> In the summary of the output of arma, there's a number Pr(>|t|), however, I
>> don't know what is its meaning - at least, it doesn't _seem_ to be a
>> Student's t distribution.
>
> It is using asymptotic normality.  There is no exact theory.  Who mentioned 
> Student's t?
>
>> Reproducible test case:
>>  x <- c(0.5, sin(1:9))
>>  reg <- arma(x, c(1,0))
>>  summary(reg)
>> 
>> <output>
>> Call:
>> arma(x = x, order = c(1, 0))
>> 
>> Model:
>> ARMA(1,0)
>> 
>> Residuals:
>>    Min      1Q  Median      3Q     Max
>> -0.9217 -0.4915  0.2254  0.4580  0.7481
>> 
>> Coefficient(s):
>>           Estimate  Std. Error  t value Pr(>|t|)
>> ar1          0.6089      0.2490    2.446   0.0145 *
>> intercept    0.0790      0.1815    0.435   0.6634
>> ---
>> Signif. codes:  0 ?***? 0.001 ?**? 0.01 ?*? 0.05 ?.? 0.1 ? ? 1
>> 
>> Fit:
>> sigma^2 estimated as 0.3348,  Conditional Sum-of-Squares = 2.68,  AIC = 
>> 21.44
>> </output>
>> 
>> Now, 2.446 is 0.6089 / 0.2490, but 0.0145 is not
>> 2 * (1 - pt(2.446, df = 7))
>> 
>> (I think there are seven degrees of freedom: the first value of
>> the series x is deterministic, and two degrees are lost in the
>> estimation of ar1 and intercept)
>
> Why is the first value deterministic?  This is not a conditional mle (see the 
> help page).

It is a conditional SSq, so perhaps you meant 'held constant'?

>> What am I misunderstanding?
>> 
>> BTW, a similar example:
>> x <- 1:10
>> y <- sin(x)
>> reg <- lm(y ~ x)
>> summary(reg)
>> 
>> will give a t-value for 'x' = 0.704 with P(>|t|) = 0.501,
>> which is 2 * (1 - pt(0.704, df=8))
>> 
>> Alberto Monteiro
>> 
>> ______________________________________________
>> R-help at 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.
>> 
>
> -- 
> Brian D. Ripley,                  ripley at stats.ox.ac.uk
> Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
> University of Oxford,             Tel:  +44 1865 272861 (self)
> 1 South Parks Road,                     +44 1865 272866 (PA)
> Oxford OX1 3TG, UK                Fax:  +44 1865 272595
>

-- 
Brian D. Ripley,                  ripley at stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595



More information about the R-help mailing list