R/ANOVA.R

In cwendorf/easi-dev: Estimation Approach to Statistical Inference

# Estimation Approach to Statistical Inference
## ANOVA Functions

### ANOVA Omnibus Function

nhstANOVA <- function(x,...) 
  UseMethod("nhstANOVA")

nhstANOVA.default <- function(...){
  data <- data.frame(...)
  columns <- dim(data)[2]
  dataLong <- reshape(data,varying=1:columns,v.names="Outcome",timevar="Variables",idvar="Subjects",direction="long")
  dataLong$Subjects <- as.factor(dataLong$Subjects)
  dataLong$Variables <- as.factor(dataLong$Variables)
  model <- lm(Outcome~Variables+Subjects,data=dataLong)
  return(anova(model))
}

nhstANOVA.formula <- function(formula,...){
  model <- lm(formula,...)
  return(anova(model))
}