Finds the Maximum Likelihood (ML) Estimates of the mean vector and variance-covariance matrix for multivariate normal data with (potentially) missing values.
mlest(data, ...)Maximum Likelihood Estimation (MLE) of the mean vector.
MLE of the variance-covariance matrix.
The objective function that is minimized by nlm.
Is is proportional to twice the negative log-likelihood.
The curvature of the likelihood surface at the MLE, in the parameterization used internally by the optimization algorithm. This parameterization is: mean vector first, followed by the log of the diagonal elements of the inverse of the Cholesky factor, and then the elements of the inverse of the Cholesky factor above the main diagonal. These off-diagonal elements are ordered by column (left to right), and then by row within column (top to bottom).
The stop code returned by nlm.
The number of iterations used by nlm.
A data frame or matrix containing multivariate normal data. Each row should correspond to an observation, and each column to a component of the multivariate vector. Missing values should be coded by 'NA'.
Optional arguments to be passed to the nlm optimization routine.
The estimate of the variance-covariance matrix returned by
mlest is necessarily positive semi-definite. Internally,
nlm is used to minimize the negative log-likelihood, so
optional arguments mayh be passed to nlm which modify the
details of the minimization algorithm, such as iterlim. The
likelihood is specified in terms of the inverse of the Cholesky factor
of the variance-covariance matrix (see Pinheiro and Bates (2000, ISBN:1441903178)).
mlest cannot handle data matrices with more than 50 variables.
Each varaible must also be observed at least once.
Little, R. J. A., and Rubin, D. B. (1987) Statistical Analysis with Missing Data. New York: Wiley, ISBN:0471802549.
Pinheiro, J. C., and Bates, D. M. (1996) Unconstrained parametrizations for variance-covariance matrices. Statistics and Computing 6, 289--296, tools:::Rd_expr_doi("10.1007/BF00140873").
Pinheiro, J. C., and Bates, D. M. (2000) Mixed-effects models in S and S-PLUS. New York: Springer, ISBN:1441903178.
library(mvnmle)
data(apple)
mlest(apple)
data(missvals)
mlest(missvals, iterlim = 400)
Run the code above in your browser using DataLab