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>.