R gam.reparam 寻找平方根惩罚的稳定正交重新参数化。

说明

用于查找正交重新参数化的内部函数，避免平方根惩罚分量之间的“主导机器零泄漏”。

用法

gam.reparam(rS, lsp, deriv)

参数

rS	平方根惩罚列表：最后一项是固定惩罚的根，如果 fixed.penalty==TRUE (即 length(rS)>length(sp) )。这里的假设是 rS[[i]] 已经处于总惩罚的零空间中；参见例如totalPenaltySpace 和 mini.roots 。
lsp	对数平滑参数的向量。
deriv	如果 deriv==1 还返回惩罚矩阵的对数行列式的一阶导数，如果 deriv>1 还返回二阶导数。

值

包含的列表

• S：为了稳定性而变换的总惩罚矩阵相似度。

• rS：分量平方根，以同样的方式转换。

• Qs ：正交变换矩阵 S = t(Qs)%*%S0%*%Qs ，其中 S0 是 sp 和 rS 对输入隐含的未变换总惩罚。

• det：日志|S|。

• det1 ：dlog|S|/dlog(sp) 如果deriv >0 。

• det2：log|S| 的粗麻布矩阵写日志(sp)如果deriv>1。

作者

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