Estimates the model parameters for all or a subset of the units.
# S3 method for SingleArrayUnitModel
fit(this, arrays=NULL, units="remaining", ..., force=FALSE, verbose=FALSE)Returns nothing.
The arrays to be fitted.
If NULL, all arrays are considered.
If remaining, only non-fitted arrays are considered.
The units to be fitted.
If NULL, all units are considered.
If remaining, only non-fitted units are considered.
Arguments passed to *readUnits().
If TRUE, already fitted units are re-fitted, and
cached data is re-read.
See Verbose.
Henrik Bengtsson
All estimates are stored to file.
The parameter estimates specific to each array, typically "chip effects", are stored in array specific files.
Array-specific estimates [K = nbr of arrays]: theta [K doubles] (chip effects), sd(theta) [K doubles], isOutlier(theta) [K logicals]
For each array and each unit group, we store: 1 theta, 1 sd(theta), 1 isOutlier(theta), i.e. (float, float, bit) => For each array and each unit (with \(G_j\) groups), we store: \(G_j\) theta, \(G_j\) sd(theta), \(G_j\) isOutlier(theta), i.e. \(G_j\)*(float, float, bit).
For more information see SingleArrayUnitModel.