phase1 function

Core function for Phase 1

phase1(
  startingestimates,
  inv.zcov,
  inv.scaling,
  z.phase1,
  z.obs,
  nodes,
  effects,
  objects,
  r.truncation.p1,
  length.p1,
  fixed.estimates,
  verbose = FALSE
)

Arguments

• startingestimates: vector containing initial parameter values
• inv.zcov: inverted covariance matrix
• inv.scaling: scaling matrix
• z.phase1: statistics retrieved from phase 1
• z.obs: observed statistics
• nodes: node set (data frame)
• effects: effects/sufficient statistics (list with a vector "names", and a vector "objects")
• objects: objects used for statistics calculation (list with a vector "name", and a vector "object")
• r.truncation.p1: numeric used to limit extreme values in the covariance matrix (for stability)
• length.p1: number of samples in phase 1
• fixed.estimates: if some parameters are fixed, list with as many elements as effects, these elements equal a fixed value if needed, or NULL if they should be estimated
• verbose: logical: should intermediate results during the estimation be printed or not? Defaults to FALSE.

Returns

estimated parameters after phase 1