包mgcv中sp.vcov的元素?

时间:2017-05-02 09:12:49

标签: r

考虑以下代码来拟合广义加性模型:

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矩阵,而模型中有两个平滑参数,所以我希望它是一个2x2矩阵?

1 个答案:

答案 0 :(得分:1)

请参阅select * from database_name.table_name;

  

如果比例参数已经(RE)ML估计(即方法是" ML"或" REML"并且比例参数未知)则最后一行和一列与对数比例参数有关。