mppe {bgms}  R Documentation 
The function mppe
estimates the parameters for the ordinal MRF
by optimizing the pseudoposterior with the NewtonRaphson method.
mppe(
x,
interaction_prior = c("Cauchy", "UnitInfo"),
cauchy_scale = 2.5,
threshold_alpha = 0.5,
threshold_beta = 0.5,
convergence_criterion = sqrt(.Machine$double.eps),
maximum_iterations = 1000,
thresholds,
interactions
)
x 
A dataframe or matrix with 
interaction_prior 
The prior distribution for the interaction effects.
Currently, two prior densities are implemented: The Unit Information prior
( 
cauchy_scale 
The scale of the Cauchy prior for interactions. Defaults
to 
threshold_alpha , threshold_beta 
The shape parameters of the Betaprime
prior for the thresholds. Default to 
convergence_criterion 
The convergence criterion for the
pseudoposterior values in the EM algorithm. Defaults to

maximum_iterations 
The maximum number of EM iterations used. Defaults
to 
thresholds 
A matrix with 
interactions 
A matrix with 
A list containing:
interactions
: A matrix with p
rows and p
columns,
containing the maximum pseudoposterior estimates of the pairwise
associations in the offdiagonal elements.
thresholds
: A matrix with p
rows and max(m)
columns, containing the maximum pseudoposterior estimates of the category
thresholds for each node.
hessian
: A square matrix with sum(m) + p(p1)/2
rows and
columns, evaluated at the maximum pseudoposterior estimates. The topleft
square contains the thresholds, the bottomright square the associations (of
the form (1,2), (1, 3), ..., (2, 1), ...
).
In the case that interaction_prior = "UnitInfo"
, the list also
contains the p
by p
matrix unit_info
, which contains the
asymptotic variances that are used to set the unit information prior for the
association effects in the bgms
function.