要检索呼叫，请执行以下操作：

as.list(mdl$call)$weights
# varIdent(form = ~1 | g * h)

要获取权重，可以使用nlme::varWeights()函数：
返回与对象表示的方差函数结构对应的标准差的倒数。

varWeights(mdl$modelStruct)
#       1*1       1*2       1*1       1*2       1*1       1*2       1*1       1*2       1*1       1*2       1*1       1*2       1*1       1*2       1*1       1*2 
# 1.0000000 0.8064094 1.0000000 0.8064094 1.0000000 0.8064094 1.0000000 0.8064094 1.0000000 0.8064094 1.0000000 0.8064094 1.0000000 0.8064094 1.0000000 0.8064094 
#       1*1       1*2       1*1       1*2       1*1       1*2       1*1       1*2       1*1       1*2       1*1       1*2       1*1       1*2       1*1       1*2 
# 1.0000000 0.8064094 1.0000000 0.8064094 1.0000000 0.8064094 1.0000000 0.8064094 1.0000000 0.8064094 1.0000000 0.8064094 1.0000000 0.8064094 1.0000000 0.8064094 
# etc

它们被存储为attr字节。

attr(mdl$modelStruct$varStruct, 'weights')
#       1*1       1*2       1*1       1*2       1*1       1*2       1*1       1*2 
# 1.0000000 0.8064094 1.0000000 0.8064094 1.0000000 0.8064094 1.0000000 0.8064094 
#       1*1       1*2       1*1       1*2       1*1       1*2       1*1       1*2 
# 1.0000000 0.8064094 1.0000000 0.8064094 1.0000000 0.8064094 1.0000000 0.8064094 
# ...

str(mdl)在这种情况下有很大的帮助。
实际上，我们可以给予"gls"一种新的方法。

weights.gls <- function (object, ...) {
  wts <- attr(object$modelStruct$varStruct, 'weights')
  if (is.null(wts)) 
    wts
  else napredict(object$na.action, wts)
}

weights(mdl)
#       1*1       1*2       1*1       1*2       1*1       1*2  ...
# 1.0000000 0.8064094 1.0000000 0.8064094 1.0000000 0.8064094  ...
# ...