corAR1
位于 nlme
包(package)。 说明
此函数是 corAR1
类的构造函数,表示 1 阶自相关结构。使用此构造函数创建的对象稍后必须使用适当的 Initialize
方法进行初始化。
用法
corAR1(value, form, fixed)
参数
value |
滞后 1 自相关的值,必须介于 -1 和 1 之间。默认为 0(无自相关)。 |
form |
|
fixed |
一个可选的逻辑值,指示是否应允许系数在优化中变化,或保持固定在其初始值。默认为 |
值
类 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.lme
、corARMA
、corClasses
、Dim.corSpatial
、Initialize.corStruct
、summary.corStruct
相关用法
- R corARMA ARMA(p,q) 相关结构
- R corRatio 有理二次相关结构
- R corLin 线性相关结构
- R corSymm 一般相关结构
- R corMatrix.pdMat 从 pdMat 对象中提取相关矩阵
- R corMatrix 提取相关矩阵
- R corFactor.corStruct corStruct 对象矩阵的因子
- R corExp 指数相关结构
- R corSpher 球面相关结构
- R corFactor 相关矩阵的因子
- R corMatrix.reStruct 从 reStructure 对象的组件中提取相关矩阵
- R corGaus 高斯相关结构
- R corNatural 自然参数化中的一般相关性
- R corCAR1 连续 AR(1) 相关结构
- R corCompSymm 复合对称相关结构
- R corSpatial 空间相关结构
- R corMatrix.corStruct corStruct 对象的矩阵
- R coef.modelStruct 提取模型结构对象系数
- R collapse.groupedData 折叠分组数据对象
- R comparePred 比较预测
- R coef.gnls 提取 glns 系数
- R coef.lmList 提取 lmList 系数
- R coef.pdMat pdMat 对象系数
- R coef.varFunc varFunc 对象系数
- R coef.lme 提取lme系数
注:本文由纯净天空筛选整理自R-devel大神的英文原创作品 AR(1) Correlation Structure。非经特殊声明,原始代码版权归原作者所有,本译文未经允许或授权,请勿转载或复制。