パッケージmgcvのsp.vcovの要素は？

Question

は、一般化加法モデルをフィッティングするための以下のコードを検討：

require(mgcv) 
attach(data) 
mod = gam(y~s(x1)+s(x2), method="REML") 

パラメータ推定値を平滑化するログの推定共分散行列は次式で与えられる。出力される

sp.vcov(model) 

：

  [,1]   [,2]  [,3] 
[1,] 0.45465578 1.551113e-01 0.02868156 
[2,] 0.15511135 4.122903e+05 0.07369381 
[3,] 0.02868156 7.369381e-02 0.02221753 

なぜ3x3行列ですか、モデルには2つの平滑化パラメータがありますので、2x2行列になると思いますか？

Answer 1

?sp.vcovを参照してください：

スケールパラメータは、（RE）ML推定されている場合（つまり、メソッドは「ML」または「REML」と尺度パラメータが不明であった場合）、最後の行と列の後、ログスケールパラメータに関連します。