This function is a constructor for the
representing a general positive-definite matrix. If the matrix
object is of dimension n, it is
represented by n*(n+1)/2 unrestricted parameters,
using the log-Cholesky parametrization described in Pinheiro and
numeric(0), an uninitialized
object, a one-sided formula, or a character vector,
returned as an uninitialized
pdLogChol object (with just
some of its attributes and its class defined) and needs to have its
coefficients assigned later, generally using the
matrix replacement functions.
value is an initialized
object will be constructed from
value is a numeric vector, it is assumed to
represent the unrestricted coefficients of the matrix-logarithm
parametrization of the underlying positive-definite matrix.
pdLogChol(value, form, nam, data)
an optional initialization value, which can be any of the
an optional one-sided linear formula specifying the
row/column names for the matrix represented by
an optional character vector specifying the row/column names
for the matrix represented by object. It must have length equal to
the dimension of the underlying positive-definite matrix and
unreplicated elements. This argument is ignored when
an optional data frame in which to evaluate the variables
pdLogChol representation of a symmetric
positive definite matrix is a vector starting with the logarithms of
the diagonal of the Choleski factorization of that matrix followed by
its upper triangular portion.
pdLogChol object representing a general positive-definite
matrix, also inheriting from class
José Pinheiro and Douglas Bates firstname.lastname@example.org
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.
(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
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.