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