As the Title says. Implemented as a method of the predict generic, for objects created by the up-to-date workflow using gaussian mixture modelling of a joint distribution of parameters and statistics (hence the newdata argument, shared by many predict methods; but these newdata should be parameter values, not data).
# S3 method for SLik_j
predict(
object, newdata, log = TRUE, which = "lik",
tstat = t(attr(object$logLs, "stat.obs")),
solve_t_chol_sigma_lists = object$clu_params$solve_t_chol_sigma_lists,
...)Numeric: a single value, or a vector of (log-)likelihoods for different rows of the input newdata.
an object of class SLik_j, as produced by infer_SLik_joint.
A matrix, whose rows each contain a full vector of the fitted parameters; or a single vector. If parameter names are not provided (as column names in the matrix case), then the vector is assumed to be ordered as object$colTypes$fittedPars.
Boolean: whether to return log-likelihood or likelihood.
"lik" or "safe". The latter protects against some artefacts of predictions beyond the regions of parameter space well sampled by the inference procedure.
The data (as projected summary statistics). Defaults to the data input in the inference procedure (i.e., the projected statistics used as stat.obs argument of infer_SLik_joint).
For programming purposes. Do not change this argument.
For consistency with the generic. Currently ignored.
## see help("example_reftable")
Run the code above in your browser using DataLab