em.control function

Control Parameters for the EM algorithm

Input parameters for the EM algorithm.

em.control(pi.initial = NULL, max.iter = 10000, tol = 1e-12,
           nr.threads = 0, verbose = TRUE)

Arguments

• pi.initial: Initial guess for the probabilities of the vectors of associations status. If NULL then 0.9 is assigned for the c(0,...,0) configuration and 0.1 is distributed uniformly for all other configurations.
• max.iter: Maximum number of EM iterations.
• tol: Tolerance (in maximum absolute difference between two EM iterations in estimated probabilities) before declaring convergence and stopping.
• nr.threads: Number of processing threads to use. If zero (the default), will automatically detect the number of compute cores available and spawn one thread per core.
• verbose: An indicator of whether to report progress (running iteration number) during computation.

Details

The function is used inside the control argument in repfdr and piem.

Returns

A list with the input values.

See Also

repfdr piem

Examples

## Not run:

download.file('http://www.math.tau.ac.il/~ruheller/repfdr_RData/binned_zmat.RData',
  destfile = "binned_zmat.RData")
load(file = "binned_zmat.RData")
out <- repfdr(pbz,bz,"replication",
              control = em.control(pi.initial = c(0.48,rep(0.02,26)),
              verbose = TRUE, nr.threads = 1))
# iterations are printed; run bit slower (1 thread)
## End(Not run)