R pdLogChol 廣義正定矩陣

說明

該函數是pdLogChol 類的構造函數，表示一般正定矩陣。如果與 object 關聯的矩陣的維度為 n ，則它由 n(n+1)/2 無限製參數表示，使用 Pinheiro 和 Bates (1996) 中說明的 log-Cholesky 參數化。

• 當 value 為 numeric(0) 、未初始化的 pdMat 對象、單邊公式或字符向量時，object 將作為未初始化的 pdLogChol 對象返回(僅定義其部分屬性和類) )並且需要稍後分配其係數，通常使用coef或matrix替換函數。

• 如果 value 是已初始化的 pdMat 對象，則 object 將從 as.matrix(value) 構造。

• 最後，如果 value 是數值向量，則假定表示基礎正定矩陣的 matrix-logarithm 參數化的無限製係數。

用法

pdLogChol(value, form, nam, data)

參數

value	可選的初始化值，可以是以下任意值：pdMat 對象、正定矩陣、單邊線性公式(變量由 + 分隔)、字符串向量或數值向量。默認為 numeric(0) ，對應於未初始化的對象。
form	可選的單邊線性公式，指定 object 表示的矩陣的行/列名稱。由於 form 中可能存在因子，因此需要在數據幀上計算公式以解析其定義的名稱。當 value 是單邊公式時，該參數將被忽略。默認為 NULL 。
nam	一個可選的字符向量，指定對象表示的矩陣的行/列名稱。它的長度必須等於基礎正定矩陣和未複製元素的維度。當 value 是字符向量時，將忽略此參數。默認為 NULL 。
data	一個可選 DataFrame ，用於評估 value 和 form 中命名的變量。它用於獲取 factors 的級別，這會影響底層矩陣的維度和行/列名稱。如果是 NULL ，則不會嘗試獲取有關公式中出現的 factor 的信息。默認為調用該函數的父框架。

細節

在內部，對稱正定矩陣的 pdLogChol 表示是一個向量，該向量以該矩陣的 Choleski 分解對角線的對數開始，後跟其上三角部分。

值

表示一般正定矩陣的 pdLogChol 對象，也繼承自類 pdMat 。

例子

(pd1 <- pdLogChol(diag(1:3), nam = c("A","B","C")))

(pd4 <- pdLogChol(1:6))
(pd4c <- chol(pd4)) # -> upper-tri matrix with off-diagonals  4 5 6
pd4c[upper.tri(pd4c)]
log(diag(pd4c)) # 1 2 3

作者

José Pinheiro and Douglas Bates bates@stat.wisc.edu

參考

Pinheiro, J.C. and Bates., D.M. (1996) Unconstrained Parametrizations for Variance-Covariance Matrices, Statistics and Computing 6, 289-296.

Pinheiro, J.C., and Bates, D.M. (2000) Mixed-Effects Models in S and S-PLUS, Springer.

也可以看看

as.matrix.pdMat , coef.pdMat , pdClasses , matrix<-.pdMat