R getVarCov 提取方差-協方差矩陣

說明

從擬合模型(例如混合效應模型)中提取方差-協方差矩陣。

用法

getVarCov(obj, ...)
## S3 method for class 'lme'
getVarCov(obj, individuals,
    type = c("random.effects", "conditional", "marginal"), ...)
## S3 method for class 'gls'
getVarCov(obj, individual = 1, ...)

參數

obj	合身的模型。方法可用於 lme 和 gls 擬合的模型
individuals	對於由 lme 擬合的模型，可以為條件矩陣或邊際方差-協方差矩陣指定分組因子級別的向量。
individual	對於 gls 擬合的模型，提供的唯一類型的方差-協方差矩陣是按組響應的邊際方差-協方差。可選參數individual指定響應組。
type	對於由 lme 擬合的模型，type 參數指定方差-協方差矩陣的類型，"random.effects" 表示 random-effects 方差-協方差(默認)，或 "conditional" 表示條件。響應的方差-協方差或"marginal" 響應的邊際方差-協方差。
...	某些方法的可選參數，如上所述

值

方差-協方差矩陣或方差-協方差矩陣列表。

例子

fm1 <- lme(distance ~ age, data = Orthodont, subset = Sex == "Female")
getVarCov(fm1)
getVarCov(fm1, individual = "F01", type = "marginal")
getVarCov(fm1, type = "conditional")
fm2 <- gls(follicles ~ sin(2*pi*Time) + cos(2*pi*Time), Ovary,
           correlation = corAR1(form = ~ 1 | Mare))
getVarCov(fm2)

作者

Mary Lindstrom lindstro@biostat.wisc.edu

也可以看看

lme , gls