当前位置: 首页>>代码示例 >>用法及示例精选 >>正文


R corAR1 AR(1) 相关结构


R语言 corAR1 位于 nlme 包(package)。

说明

此函数是 corAR1 类的构造函数,表示 1 阶自相关结构。使用此构造函数创建的对象稍后必须使用适当的 Initialize 方法进行初始化。

用法

corAR1(value, form, fixed)

参数

value

滞后 1 自相关的值,必须介于 -1 和 1 之间。默认为 0(无自相关)。

form

~ t~ t | g 形式的单边公式,指定时间协变量 t 和可选的分组因子 g 。此相关结构的协变量必须是整数值。当 form 中存在分组因子时,假定相关结构仅适用于同一分组级别内的观测值;假设具有不同分组级别的观测值是不相关的。默认为 ~ 1 ,这对应于使用数据中的观察顺序作为协变量,并且没有组。

fixed

一个可选的逻辑值,指示是否应允许系数在优化中变化,或保持固定在其初始值。默认为 FALSE ,在这种情况下允许系数变化。

corAR1 的对象,表示 1 阶自相关结构。

例子

## covariate is observation order and grouping factor is Mare
cs1 <- corAR1(0.2, form = ~ 1 | Mare)

# Pinheiro and Bates, p. 236
cs1AR1 <- corAR1(0.8, form = ~ 1 | Subject)
cs1AR1. <- Initialize(cs1AR1, data = Orthodont)
corMatrix(cs1AR1.)

# Pinheiro and Bates, p. 240
fm1Ovar.lme <- lme(follicles ~ sin(2*pi*Time) + cos(2*pi*Time),
                   data = Ovary, random = pdDiag(~sin(2*pi*Time)))
fm2Ovar.lme <- update(fm1Ovar.lme, correlation = corAR1())

# Pinheiro and Bates, pp. 255-258:  use in gls
fm1Dial.gls <-
  gls(rate ~(pressure + I(pressure^2) + I(pressure^3) + I(pressure^4))*QB,
      Dialyzer)
fm2Dial.gls <- update(fm1Dial.gls,
                 weights = varPower(form = ~ pressure))
fm3Dial.gls <- update(fm2Dial.gls,
                    corr = corAR1(0.771, form = ~ 1 | Subject))

# Pinheiro and Bates use in nlme:  
# from p. 240 needed on p. 396
fm1Ovar.lme <- lme(follicles ~ sin(2*pi*Time) + cos(2*pi*Time),
                   data = Ovary, random = pdDiag(~sin(2*pi*Time)))
fm5Ovar.lme <- update(fm1Ovar.lme,
                corr = corARMA(p = 1, q = 1))
# p. 396
fm1Ovar.nlme <- nlme(follicles~
     A+B*sin(2*pi*w*Time)+C*cos(2*pi*w*Time),
   data=Ovary, fixed=A+B+C+w~1,
   random=pdDiag(A+B+w~1),
   start=c(fixef(fm5Ovar.lme), 1) )
# p. 397
fm2Ovar.nlme <- update(fm1Ovar.nlme,
         corr=corAR1(0.311) )

作者

José Pinheiro and Douglas Bates bates@stat.wisc.edu

参考

Box, G.E.P., Jenkins, G.M., and Reinsel G.C. (1994) "Time Series Analysis: Forecasting and Control", 3rd Edition, Holden-Day.

Pinheiro, J.C., and Bates, D.M. (2000) "Mixed-Effects Models in S and S-PLUS", Springer, esp. pp. 235, 397.

也可以看看

ACF.lmecorARMAcorClassesDim.corSpatialInitialize.corStructsummary.corStruct

相关用法


注:本文由纯净天空筛选整理自R-devel大神的英文原创作品 AR(1) Correlation Structure。非经特殊声明,原始代码版权归原作者所有,本译文未经允许或授权,请勿转载或复制。