coord2K: A coordinate descent algorithm to find the maximum likelihood...

Description Usage Arguments Value Examples

View source: R/coord2K.R

Description

This function implements a simple coordinate descent algorithm to find the maximum likelihood estimator over Gaussian MTP2 distributions. For details see Lauritzen, Uhler, Zwiernik (2017).

Usage

1
coord2K(S, n = 1, K0 = solve(diag(diag(S))), tol = 1e-08)

Arguments

S

the sample covariance matrix

n

the sample size (default 1), relevant for testing

K0

the starting point (default K0=solve(diag(diag(S))))

tol

the convergence tolerance (default tol=1e-8)

Value

the optimal value of the concentration matrix

the number of iterations the algorithm needed to converge

the corresponding value of the log-likelihood

Examples

1

pzwiernik/mtp2 documentation built on Aug. 9, 2020, 12:34 p.m.