R ldetS 獲取懲罰矩陣的對數廣義行列式

說明

內部函數計算按塊存儲在 Sl 列表中的懲罰矩陣 S 的對數廣義行列式(這是 Sl.setup 的輸出)。

用法

ldetS(Sl, rho, fixed, np, root = FALSE, repara = TRUE,
      nt = 1,deriv=2,sparse=FALSE)

參數

Sl	Sl.setup 的輸出。
rho	對數平滑參數。
fixed	指示平滑參數是否固定(或自由)的數組。
np	係數的數量。
root	指示是否返回總懲罰 S_tot 的矩陣平方根 E 。
repara	如果為 TRUE，multi-term 塊將使用 gam.reparam 重新參數化，並在返回的對象中提供重新參數化對象。
nt	要使用的並行線程數。
deriv	導數的使用順序
sparse	E 應該稀疏嗎？

值

包含以下內容的列表：

• ldetS：S 的對數行列式。

• ldetS1：S 的對數行列式的梯度。

• ldetS2：S 的對數行列式的 Hessian 矩陣。

• Sl：如果需要，修改 rS 項並將 rho 添加到每個塊

• rp：重新參數化列表。

• rp ：E 總懲罰平方根，使得 t(E)%*%E = S_tot (如果 root==TRUE )。

作者

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