R gamlss.gH 計算回歸係數的對數似然導數

R語言 gamlss.gH 位於 mgcv 包(package)。

說明

主要供位置比例模型係列內部使用。給定對數似然相對於線性預測器的導數，該函數獲得導數和 Hessian 相對於回歸係數和 Hessian 相對於回歸係數的導數。平滑參數。有關輸入導數數組打包約定，請參閱 trind.generator 。

用法

gamlss.gH(X, jj, l1, l2, i2, l3 = 0, i3 = 0, l4 = 0, i4 = 0, d1b = 0,
  d2b = 0, deriv = 0, fh = NULL, D = NULL,sandwich=FALSE)

參數

`X`	包含所有線性預測變量的模型矩陣的矩陣。
`jj`	索引向量列表，使得 `X[,jj[[i]]]` 是 i-th 線性預測器的模型矩陣。
`l1`	每個參數的對數似然的每個元素的一階導數數組。
`l2`	每個參數的對數似然的每個元素的二階導數數組。
`i2`	二維索引數組，使得 `l2[,i2[i,j]]` 包含部分 w.r.t. params 由 i,j 索引，對索引值沒有限製(除了它們位於 1,...,ncol(l1) 中)。
`l3`	每個參數的對數似然的每個元素的三階導數數組。
`i3`	third-dimensional 索引數組，使得 `l3[,i3[i,j,k]]` 包含部分 w.r.t.參數由 i,j,k 索引。
`l4`	每個參數的對數似然的每個元素的四階導數數組。
`i4`	third-dimensional 索引數組，使得 `l4[,i4[i,j,k,l]]` 包含部分 w.r.t.參數由 i,j,k,l 索引。
`d1b`	回歸係數與平滑參數的一階導數。
`d2b`	回歸係數關於平滑參數的二階導數。
`deriv`	如果`deriv==0` 將僅計算一階和二階導數。如果`deriv==1`，該函數還返回Hessian 矩陣的一階導數的對角線，如果`deriv==2`，它返回完整的三階導數，如果`deriv==3`，它還提供四階導數。
`fh`	懲罰 Hessian 矩陣的特征分解或 Cholesky 因子。
`D`	對角矩陣，用於提供一些縮放。
`sandwich`	設置為 `TRUE` 以返回三明治估計器 'filling'，而不是 `l2` 中的 Hessian 矩陣。

值

包含 lb 的列表 - 梯度向量 w.r.t.係數； lbb - Hessian 矩陣 w.r.t.係數； d1H - Hessian 矩陣的導數列表平滑參數，或單個矩陣，其列是這些導數矩陣的前導對角線； trHid2H - 逆 Hessian 矩陣的跡乘以 Hessian 矩陣的二階導數。平滑參數的所有組合。

作者

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

也可以看看

trind.generator

相關用法

注：本文由純淨天空篩選整理自R-devel大神的英文原創作品 Calculating derivatives of log-likelihood wrt regression coefficients。非經特殊聲明，原始代碼版權歸原作者所有，本譯文未經允許或授權，請勿轉載或複製。