singleLineAnova {wrMisc}

2-factorial Anova on single line of data

Description

This function runs 2-factorial Anova on a single line of data (using aov from package stats) using a model with two factors (without factor-interaction) and extracts the correpsonding p-value.

Usage

singleLineAnova(
  dat,
  fac1,
  fac2,
  inclInteraction = TRUE,
  silent = FALSE,
  debug = FALSE,
  callFrom = NULL
)

Arguments

dat

numeric vector

fac1

(character or factor) vector describing grouping elements of dat for first factor, must be of same langth as fac2

fac2

(character or factor) vector describing grouping elements of dat for second factor, must be of same langth as fac1

inclInteraction

(logical) decide if factor-interactions (eg synergy) should be included to model

silent

(logical) suppress messages

debug

(logical) additional messages for debugging

callFrom

(character) allow easier tracking of messages produced

Value

This function returns the (uncorrected) p for factor 'Pr(>F)' (see aov)

See Also

aov, anova; for repeated tests using the package limma including lmFit and eBayes see test2factLimma

Examples

set.seed(2012); dat <- round(runif(8),1)
singleLineAnova(dat, gl(2,4),rep(1:2,4))