Re: [R] continuous time AR(1)

From: peter dalgaard <pdalgd_at_gmail.com>
Date: Sat, 28 May 2011 14:01:38 +0200

On May 27, 2011, at 18:00 , frossard victor wrote:

> Dear R helpers,
> I would like to model temporal trend of biological remains in sediment cores.
> All samples are temporally auto-correlated and I would like to take this effect
> into account.
> Initially I thought that I could use AR(1) or ARIMA functions but these
> functions only work with regular temporal intervals between samples.
> Hence I would like to use a continuous time AR(1) that allow irregular time
> intervals between samples.
> Unfortunalty I don't find this function in any R package.
> Doesn't someone know if this function has already be implemented for R?

library(nlme)
help(corCAR1)
help(corExp)

You use these corStructs with gls/lme/nlme.

As far as I remember, the corCAR1 is equivalent to corExp without nugget effect (modulo parametrization?), but you often want the nugget effect if some data are very close in time without actually being identical (since corCAR1 implies that the correlation is one for simultaneous observations).

I wouldn't know about continuous time AR of higher orders (I can't even guess what MA(n) might mean in continuous time).

-pd

> Many thanks in advance.
> Victor Frossard
> Phd. student.
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> 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.

-- 
Peter Dalgaard
Center for Statistics, Copenhagen Business School
Solbjerg Plads 3, 2000 Frederiksberg, Denmark
Phone: (+45)38153501
Email: pd.mes_at_cbs.dk  Priv: PDalgd_at_gmail.com

______________________________________________
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.
Received on Sat 28 May 2011 - 12:04:31 GMT

This quarter's messages: by month, or sorted: [ by date ] [ by thread ] [ by subject ] [ by author ]

All messages

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Sat 28 May 2011 - 13:30:10 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.

list of date sections of archive