[R] Default lag.max in ACF
Prof Brian Ripley
ripley at stats.ox.ac.uk
Tue Mar 28 08:19:19 CEST 2006
The default is taken from S-PLUS, so the reference is the S-PLUS manual.
It is pretty similar to the recommendation of Brockwell & Davis.
On Mon, 27 Mar 2006, Spencer Graves wrote:
> I don't know why the default lag.max is 10*log10(N/m) for m series.
> The "acf" help page includes the following:
>
> Author(s):
>
> Original: Paul Gilbert, Martyn Plummer. Extensive modifications
> and univariate case of 'pacf' by B.D. Ripley.
>
> I've copied these three contributors on this email. With luck, one
> of them will enlighten us both.
>
> Best Wishes,
> spencer graves
>
> Rafal Stankiewicz wrote:
>
>> Hi,
>>
>> The default value for lag.max in ACF implementation is 10*log10(N)
>>
>> There several publications recommending setting lag.max to:
>> - N/4 (Box and Jenkins, 1970; Chatfield, 1975; Anderson, 1976; Pankratz,
>> 1983; Davis, 1986; etc.)
>> - sqrt(N)+10 (Cryer, 1986)
>> - 20<=N<=40 (Brockwell and Davis)
>>
>> Why R uses 10*log10(N) as a default?
>>
>> Please, give me a reference to a book or article where the recommendation
>> for using lag.max=10*log10(N) is proposed and explained.
>>
>> Thanks
>> Rafal
--
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