R vcov.gam 從 GAM 擬合中提取參數(估計器)協方差矩陣

說明

從擬合的 gam 對象中提取參數的貝葉斯後驗協方差矩陣或參數估計量的頻率協方差矩陣。

用法

## S3 method for class 'gam'
vcov(object, sandwich=FALSE, freq = FALSE, dispersion = NULL,unconditional=FALSE, ...)

參數

object	由 gam() 生成的 gam 類的擬合模型對象。
sandwich	計算協方差矩陣的三明治估計。目前，離散式 bam 配合的價格昂貴。
freq	TRUE 返回參數估計量的頻率協方差矩陣，FALSE 返回參數的貝葉斯後驗協方差矩陣。後一個選項包括根據貝葉斯平滑先驗的預期平方偏差。
dispersion	色散參數的值：通常不使用。
unconditional	如果 TRUE (和 freq==FALSE )，則返回貝葉斯平滑參數不確定性校正協方差矩陣(如果可用)。
...	其他參數，目前被忽略。

細節

本質上，隻是從 gamObject 中提取 object$Ve 、 object$Vp 或 object$Vc (如果可用)，除非 sandwich==TRUE 在這種情況下計算夾心估計(有或沒有平方偏差分量)。

值

與模型參數估計量/係數的估計頻率協方差矩陣相對應的矩陣，或參數的估計後驗協方差矩陣，具體取決於參數 freq 。

例子

 
require(mgcv)
n <- 100
x <- runif(n)
y <- sin(x*2*pi) + rnorm(n)*.2
mod <- gam(y~s(x,bs="cc",k=10),knots=list(x=seq(0,1,length=10)))
diag(vcov(mod))

作者

Henric Nilsson. Maintained by Simon N. Wood simon.wood@r-project.org

參考

Wood, S.N. (2017) Generalized Additive Models: An Introductio with R (2nd ed) CRC Press

也可以看看

gam