mixedModel {rADA}

Mixed model wrapper for assay dataframe

Description

This function is a wrapper for the lmer() function to produce a table of results. Here, the sample ID is included as a random intercept effect, then the results of the fixed effect model estimates are reported together with 95% confidence intervals and t statistics.

Usage

mixedModel(assay.df.melted, var = c("Operator", "Day"))

Arguments

assay.df.melted

A data.frame produced by assayMelt()

var

Variable to look at. Either "Day" or "Operator".

Value

A data.frame with the following columns: "Parameter", "Estimate", "LowerCI", "UpperCI", "Tstat"

Author(s)

Emma Gail

Lidija Turkovic

Examples

assay.df.melted <- assayMelt(assay.df = lognormAssay, exp.name = 'Experiment1')
mixedModel(assay.df.melted, var = 'Day')