如何从 gls 对象中提取 delta 和 sigma 的 SE 值？答案

【问题标题】：How I extract the SE values of delta and sigma from gls object?如何从 gls 对象中提取 delta 和 sigma 的 SE 值？
【发布时间】：2017-06-13 17:49:45
【问题描述】：

如何提取 gls 目标的 delta（幂函数）和 sigma 的 SE？如果我做总结，我有这个输出：

> summary(l3s)
Generalized least squares fit by REML
  Model: biom.kg ~ I(dbh^2 * h) 
  Data: ssf 
      AIC      BIC    logLik
  4495.75 4512.982 -2243.875

Variance function:
 Structure: Power of variance covariate
 Formula: ~I(dbh^2 * h) 
 Parameter estimates:
    power 
0.9383058 

Coefficients:
                  Value  Std.Error  t-value p-value
(Intercept)  -0.7649061 0.19920153 -3.83986   1e-04
I(dbh^2 * h)  0.0312859 0.00066244 47.22858   0e+00


 Correlation: 
             (Intr)
I(dbh^2 * h) -0.541


Standardized residuals:
        Min          Q1         Med          Q3         Max 
-1.98250304 -0.82199784 -0.07727941  0.69239313  3.05748659 

Residual standard error: 0.02026094 
Degrees of freedom: 551 total; 549 residual

谢谢。

【问题讨论】：

标签： r nlme

【解决方案1】：

无法提取方差协变量的标准误差，因为它不是由 gls 计算的。考虑使用蒙特卡罗模拟并提取fit$modelStruct[2]。您可以在 l3s$sigma 中找到您的模型的 sigma。

【讨论】：