limma_fit: function to retrieve the contrasts stored in a dataset and...

In geantonicelli/GLM: set of wrapper and adapter functions for statistical analysis

Description

this function is a wrapper and an adapter around 'aov', it retrieves the contrasts for each independent variable specified in an aov() model and stored in the original dataset, and fits an ANOVA model on all contrasts for each independent variable, in order to obtained the computed SS and MS for further calculations, otherwise not available from the original ANOVA model as displayed by summary() or summary.lm(), the model fitted for all contrasts of one variable is the same as displayed by summary.lm() on the original ANOVA model

Usage

limma_fit(
  object,
  design = NULL,
  ndups = 1,
  spacing = 1,
  block = NULL,
  correlation,
  weights = NULL,
  method = "ls",
  ...
)

Arguments

anova

an object of class anova data.frame with the ANOVA model created by aov(), in the original dataset each variable must be in only one column, and all contrasts for each variable should be loaded as contrasts for that variable

Value

the function returns the fitted ANOVA models for all the contrasts stored for each variable

Author(s)

gerardo esteban antonicelli

See Also

'check_contrasts' 'omega_factorial'

Examples

data(gogglesData)
data(depressionData)
goggles.model <- aov(attractiveness ~ gender + alcohol + gender:alcohol, data=gogglesData)
simple.model <- aov(attractiveness ~ simple, data=gogglesData)
depression.model <- aov(diff ~ treat, data=depressionData)
fit_aov(goggles.model)
fit_aov(simple.model)
fit_aov(depression.model)

geantonicelli/GLM documentation built on Dec. 15, 2020, 3:05 p.m.