R Sl.setup 设置表示块对角惩罚矩阵的列表

说明

内部函数，用于设置表示来自 gam.setup 生成的对象的块对角罚矩阵的列表。

用法

Sl.setup(G,cholesky=FALSE,no.repara=FALSE,sparse=FALSE)

参数

G	gam.setup 的输出。
cholesky	仅使用 Cholesky 重新参数化。
no.repara	设置为 TRUE 以关闭所有初始重新参数化。
sparse	稀疏设置？

值

每个块都有一个元素的列表。对于块 b，Sl[[b]] 是包含以下元素的列表

• repara：是否应该将重新参数化应用于模型矩阵等？如果系数是非线性的，通常为FALSE。

• start, stop ：这样start:stop就是该块的参数索引。

• S：块的惩罚矩阵列表(dim = stop-start+1)如果length(S)==1，那么这将是身份惩罚。否则，它是多重惩罚，并且将添加平方根惩罚矩阵的 rS 列表。 S (如果 repara==TRUE )和 rS (始终)将被投影到总惩罚矩阵的范围空间中。

• rS：如果使用多个惩罚，则惩罚矩阵的平方根。

• D：块的重新参数化矩阵。适用于 start:stop 中的 cols/params 。如果是数字，则 X[,start:stop]%*%diag(D) 是 X[,start:stop] 和 b.orig = D*b.repara 的重新参数化(其中 b.orig 是原始参数向量)。如果是矩阵，则 X[,start:stop]%*%D 是 X[,start:stop] 和 b.orig = D%*%b.repara 的重新参数化(其中 b.orig 是原始参数向量)。

作者

Simon N. Wood <simon.wood@r-project.org>.