Re: [R] Maximum likelihood acf

From: Prof Brian Ripley <ripley_at_stats.ox.ac.uk>
Date: Fri 12 Jan 2007 - 15:40:13 GMT

On Fri, 12 Jan 2007, Alain Guillet wrote:

> Prof. Brian Ripley,
>
> You are right, my question was not clear.
>
> In fact, I want to estimate the k first components of the acf, i.e. I
> want to estimate the k parameters (c(0),c(1),...c(k-1)), where c is the
> autocorrelation function, by a maximum likelihood estimator.

And does ARMAacf applied to the result of ar.mle not do just that? An accessible reference would help us, if not.

>
> Alain
>
>
>
> Prof Brian Ripley a écrit :
>> You will need to give us a reference, as the acf is not a parameter in
>> a model in your description and MLEs apply to model parameters.
>>
>> Just possibly ar.mle is what you are looking for, perhaps plus ARMAacf?
>>
>> On Fri, 12 Jan 2007, Alain Guillet wrote:
>>
>>> Hello!
>>>
>>> I am looking for a function which computes the maximum likelihood
>>> estimator of the autocorrelation function for a gaussian time series.
>>> Does a such function already exist in R?
>>> The estimator by default in R, acf(), uses the method of moments.
>>>
>>> Thanks a lot,
>>> Alain
>>>
>>>
>>>
>>
>
>

-- 
Brian D. Ripley,                  ripley@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


______________________________________________ R-help@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.

Received on Sat Jan 13 02:49:16 2007

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Sat 13 Jan 2007 - 02:30:25 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.