Estimate Diagonal of H — cpp

Compute a sample estimate of the expected negative Hessian by taking advantage of the second Bartlett's identity at the single pair level

Usage

cpp_DH(
  C_VEC,
  A,
  CONSTRLOGSD,
  LLC,
  THETA,
  FREQ,
  N,
  CORRFLAG,
  NTHR,
  NLOAD,
  NCORR,
  NVAR
)

C_VEC: Vector containing the number of categories for each item
A: Constraint matrix. Loadings free to be estimated are identified by a 1.
CONSTRLOGSD: \(q\)-dimensional vector. Elements set to NA refers to free latent log standard deviations parameters. Elements set to numerical values denote fixed values constraints.
LLC: Linear loadings constraints. Expects a list of constraints. See fit_plFA documentation.
THETA: Parameter vector
FREQ: output from pairs_freq()
N: sample size
CORRFLAG: 1 to estimate latent correlations. 0 for orthogonal latent factors.
NTHR: Number of thresholds parameters.
NLOAD: Number of free loadings parameters
NCORR: Number of free latent correlations parameters.
NVAR: Number of free latent variance parameters.