R corSymm 一般相關結構

說明

該函數是corSymm類的構造函數，代表一般的相關結構。該結構的內部表示，就無約束參數而言，使用 Pinheiro 和 Bates (1996) 中定義的球形參數化。使用此構造函數創建的對象稍後必須使用適當的 Initialize 方法進行初始化。

用法

corSymm(value, form, fixed)

參數

value	帶有參數值的可選向量。默認值為 numeric(0) ，這會導致在初始化 object 時將適當維度的零向量分配給參數(對應於身份相關結構)。
form	~ t 或 ~ t | g 形式的單邊公式，指定時間協變量 t 和可選的分組因子 g 。此相關結構的協變量必須是整數值。當 form 中存在分組因子時，假定相關結構僅適用於同一分組級別內的觀測值；假設具有不同分組級別的觀測值是不相關的。默認為 ~ 1 ，這對應於使用數據中的觀察順序作為協變量，並且沒有組。
fixed	一個可選的邏輯值，指示是否應允許係數在優化中變化，或保持固定在其初始值。默認為 FALSE ，在這種情況下允許係數變化。

值

代表一般相關結構的類 corSymm 的對象。

例子

## covariate is observation order and grouping factor is Subject
cs1 <- corSymm(form = ~ 1 | Subject)

# 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)

# example gls(..., corSpher ...)
# Pinheiro and Bates, pp. 261, 263
fm1Wheat2 <- gls(yield ~ variety - 1, Wheat2)
# p. 262 
fm2Wheat2 <- update(fm1Wheat2, corr =
   corSpher(c(28, 0.2),
     form = ~ latitude + longitude, nugget = TRUE))

# example gls(..., corSymm ... )
# Pinheiro and Bates, p. 251
fm1Orth.gls <- gls(distance ~ Sex * I(age - 11), Orthodont,
                   correlation = corSymm(form = ~ 1 | Subject),
                   weights = varIdent(form = ~ 1 | age))

作者

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

參考

Pinheiro, J.C. and Bates., D.M. (1996) "Unconstrained Parametrizations for Variance-Covariance Matrices", Statistics and Computing, 6, 289-296.

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

也可以看看

Initialize.corSymm、summary.corSymm