A N by p outcome data from N clusters/batches/experiments; and p is the number of samples within each cluster.
The first sample within each cluster is a reference sample with a different error variance than other samples. Missing values are coded as NAs. Note the model allows unbalanced data.
X
A covariate array of dimension of N by k by p, where k is the number of covariates.