R corMatrix.corStruct corStruct 对象的矩阵

说明

此方法函数提取与 covariate 和 object 对应的相关矩阵(或其转置逆平方根因子)或相关矩阵列表(或其转置逆平方根因子)。令\Sigma表示相关矩阵，\Sigma的平方根因子是任何方阵L使得\Sigma = L'L。当 corr = FALSE 时，此方法提取 L^{-t} 。

用法

## S3 method for class 'corStruct'
corMatrix(object, covariate, corr, ...)

参数

object	继承自类 "corStruct" 的对象，表示相关结构。
covariate	可选的协变量向量(矩阵)或协变量向量(矩阵)列表，将在其值处评估相关矩阵或相关矩阵列表。默认为 getCovariate(object) 。
corr	一个逻辑值。如果 TRUE 该函数返回相关矩阵或相关矩阵列表，由 object 表示。如果FALSE，该函数返回相关矩阵的转置逆平方根，或相关矩阵的转置逆平方根因子列表。
...	该泛型的某些方法需要额外的参数。此方法中没有使用任何内容。

值

如果covariate是向量(矩阵)，则返回值将是具有相应相关矩阵(或其转置逆平方根因子)的数组。如果 covariate 是向量(矩阵)列表，则返回值将是与 covariate 的每个分量相对应的相关矩阵(或其转置逆平方根因子)的列表。

例子

cs1 <- corAR1(0.3)
corMatrix(cs1, covariate = 1:4)
corMatrix(cs1, covariate = 1:4, corr = FALSE)

# Pinheiro and Bates, p. 225
cs1CompSymm <- corCompSymm(value = 0.3, form = ~ 1 | Subject)
cs1CompSymm <- Initialize(cs1CompSymm, data = Orthodont)
corMatrix(cs1CompSymm)

# Pinheiro and Bates, p. 226
cs1Symm <- corSymm(value = c(0.2, 0.1, -0.1, 0, 0.2, 0),
                   form = ~ 1 | Subject)
cs1Symm <- Initialize(cs1Symm, data = Orthodont)
corMatrix(cs1Symm)

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

# Pinheiro and Bates, p. 237 
cs1ARMA <- corARMA(0.4, form = ~ 1 | Subject, q = 1)
cs1ARMA <- Initialize(cs1ARMA, data = Orthodont)
corMatrix(cs1ARMA)

# Pinheiro and Bates, p. 238 
spatDat <- data.frame(x = (0:4)/4, y = (0:4)/4)
cs1Exp <- corExp(1, form = ~ x + y)
cs1Exp <- Initialize(cs1Exp, spatDat)
corMatrix(cs1Exp)

作者

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

参考

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

也可以看看

corFactor.corStruct , Initialize.corStruct