考虑以下代码来拟合广义加性模型:
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矩阵?
答案 0 :(得分:1)
请参阅select * from database_name.table_name;
:
如果比例参数已经(RE)ML估计(即方法是" ML"或" REML"并且比例参数未知)则最后一行和一列与对数比例参数有关。