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

Question

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

> 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

谢谢。

Answer 1

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