Observed-data loglikelihood for normal data
Evaluates the observed-data loglikelihood function at a user-supplied value of the parameter. This function is useful for monitoring the progress of EM and data augmentation.
loglik.norm(s, theta)
s |
summary list of an incomplete normal data matrix created by the
function |
theta |
vector of normal parameters expressed on transformed scale in packed
storage, such as one produced by the function |
value of the observed-data loglikelihood
See Section 5.3.5 of Schafer (1996)
data(mdata) s <- prelim.norm(mdata) #do preliminary manipulations thetahat <- em.norm(s) #compute MLE loglik.norm(s,thetahat) #loglikelihood at the MLE
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.