mixedModel2Fit_multiple_design: Fitting mixed models of one random effect
In jhsiao999/Humanzee: Tools for statistical analysis of genomics data

Description Usage Arguments Examples

Adapted from statmod package's mixedModel2Fit function. This function includes the additional feature of fitting a different design matrix for each gene (row) of the input data matrix.

1	mixedModel2Fit_multiple_design(design, block, yy)

`design`	Design matrix.
`block`	Phentype for the blocking factor (vector).
`yy`	Gene (feature) by sample matrix.

mixedModel2Fit_multiple_design()
data_matrix <- pheno_data
design <- model.matrix(~ 1+ as.factor(individual) + protein, 
                       data = pheno_data)
block <- as.factor(pheno_data$biorep)
yy <- pheno_data$phos