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