aov.car: Convenience wrappers for car::Anova using either a formula or factor based interface.

Description

These functions allow convenient access to Anova (from the car package) for data in the long format (i.e., one observation per row), possibly aggregating the data if there is more than one obersvation per individuum and cell. Hence, mixed between-within ANOVAs can be calculated conveniently without using the rather unhandy format of car::Anova. aov.car can be called using a formula similar to aov specifying an error strata for the within-subject factor(s). ez.glm is called specifying the factors as character vectors.

Usage

aov.car(formula, data, fun.aggregate = NULL, type = 3,
    return = "nice", observed = NULL, args.return = list(),
    ...)

  ez.glm(id, dv, data, between = NULL, within = NULL,
    covariate = NULL, observed = NULL, fun.aggregate =
    NULL, type = 3, return = "nice", args.return = list(),
    ..., print.formula = FALSE)

  univ(object)

Arguments

formula

A formula specifying the ANOVA model similar to aov. Should include an error term (i.e., Error( / )). Note that the within-subject factors do not need to be outside the Error term (this

character vector (of length 1) indicating the subject identifier column in data.

character vector (of length 1) indicating the column containing the dependent variable in data.

between

character vector indicating the between-subject(s) factor(s)/column(s) in data. Default is NULL indicating no between-subjects factors.

within

character vector indicating the within-subject(s) factor(s)/column(s) in data. Default is NULL indicating no within-subjects factors.

covariate

character vector indicating the between-subject(s) covariate(s) (i.e., column(s)) in data. Default is NULL indicating no covariates.

observed

character vector indicating which of the variables are observed (i.e, measured) as compared to experimentally manipulated. The default behavior is to return a ANOVA table with generalized eta squared effect size for which this inf

data

A data.frame containing the data. Mandatory.

fun.aggregate

The function for aggregating the data before running the ANOVA if there is more than one obervation per individuum and cell of the design. The default NULL issues a warning if aggregation is necessary and uses

type

The type of sums of squares for the ANOVA. Defaults to 3. Passed to Anova. Possible values are "II", "III", 2, or 3.

print.formula

ez.glm is a wrapper for aov.car. This boolean argument indicates whether the formula in the call to car.aov should be printed.

return

What should be returned? If "nice" (the default) will return a nice ANOVA table (produced by nice.anova. Possible values are

c("Anova", "lm", "data", "nice", "full", "all",

args.return

list of further arguments passed to the function which produces the return value. Currently only supports return = "nice" (the default) which then passes arguments to nice.anov

...

Further arguments passed to fun.aggregate.

object

An object of class Anova.mlm as returned by aov.car, ez.glm, or Anova.

Value

aov.car and ez.glm are wrappers to Anova, the return value is dependent on the return argument. When argument return is "nice" (the default) a nice ANOVA table is returnd (nice.anova). If return = "full" or return = "all" a list list with the following elements: [object Object],[object Object],[object Object],[object Object] If return = "univariate" the object returned from univ. univ returns a list of data.frames containing the univariate results (i.e., the classical ANOVA results) from an object of class "Anova.mlm". This is essentially the output from summary.Anova.mlm with multivariate = FALSE, e.g. summary(aov.car(...), multivriate = FALSE), as a list instead of printed to the console. For objects of class "anova" (i.e., the object returned by car::Anova for a purely between-subjects ANOVA) the object is returned unaltered. The elements of the list returned by univ are: anova, mauchly, and spehricity.correction (containing both, Greenhouse-Geisser and Hyundt-Feldt correction).

Details

Type 3 sums of squares are default in afex. Note that type 3 sums of squares are said to be dangerous and/or problematic. On the other side they are the default in in SPSS and SAS and recommended by e.g. Maxwell and Delaney (2004). For a brief discussion see http://stats.stackexchange.com/q/6208/442{here}. However, note that lower order effects (e.g., main effects) in type 3 ANOVAs are only meaningful with http://www.ats.ucla.edu/stat/mult_pkg/faq/general/effect.htm{effects coding}. That is, contrasts should be set to contr.sum via options(contrasts=c('contr.sum','contr.poly')). This should be done automatically when loading afex and afex will issue a warning when running type 3 SS and http://www.ats.ucla.edu/stat/r/library/contrast_coding.htm{other coding schemes}. You can check the coding with options("contrasts"). The formula for aov.car must contain a single Error term specyfying the ID column and potential within-subject factors (you may use mixed with multiple error terms). Factors outside the Error term are treated as between-subject factors (the within-subject factors specified in the Error term are ignored outside the Error term, i.e., it is not necessary to specify them outside the Error term, see Examples). Suppressing the intercept (i.e, via 0 + or

-
  1

) is ignored. Specific specifications of effects (e.g., excluding terms with - or using ^) could be okay but is not tested. Using the I or poly function within the formula is not tested and not supported! For ez.glm either between or within must not be NULL. ez.glm will concatante all between-subject factors using * (i.e., producing all main effects and interactions) and all covariates by + (i.e., adding only the main effects to the existing between-subject factors). The within-subject factors do fully interact with all between-subject factors and covariates. This is essentially identical to the behavior of SPSS's glm function. Note that the default behavior is to return a nice.anova data.frame. This includes calculation of generalized eta squared for which all non manipluated (i.e., observed) variables need to be specified via the observed argument. Changing the effect size to "pes" (partial eta-squared) via args.return or the return value via return removes this necessity.

References

Maxwell, S. E., & Delaney, H. D. (2004). Designing Experiments and Analyzing Data: A Model-Comparisons Perspective. Mahwah, N.J.: Lawrence Erlbaum Associates.

Examples

Run this code

# exampel using obk.long (see ?obk.long), a long version of the OBrienKaiser dataset from car.

data(obk.long, package = "afex")

# run univariate mixed ANCOVA for the full design:
aov.car(value ~ treatment * gender + age + Error(id/phase*hour), data = obk.long, observed = c("gender", "age"))

ez.glm("id", "value", obk.long, between = c("treatment", "gender"), within = c("phase", "hour"), covariate = "age", observed = c("gender", "age"))

# both calls return the same:
##                         Effect          df   MSE         F  ges     p
## 1                    treatment        2, 9 23.96    3.58 +  .20   .07
## 2                       gender        1, 9 23.96    3.95 +  .14   .08
## 3                          age        1, 9 23.96      0.52  .02   .49
## 4             treatment:gender        2, 9 23.96      1.28  .09   .32
## 5                        phase  1.7, 15.28  3.91 20.28 ***  .17 <.001
## 6              treatment:phase 3.39, 15.28  3.91   6.07 **  .11  .005
## 7                 gender:phase  1.7, 15.28  3.91      0.25 .002   .75
## 8                    age:phase  1.7, 15.28  3.91    3.10 +  .03   .08
## 9       treatment:gender:phase 3.39, 15.28  3.91      1.60  .03   .23
## 10                        hour 2.14, 19.23  2.48 20.52 ***  .14 <.001
## 11              treatment:hour 4.27, 19.23  2.48      0.71  .01   .60
## 12                 gender:hour 2.14, 19.23  2.48      0.71 .006   .51
## 13                    age:hour 2.14, 19.23  2.48    2.82 +  .02   .08
## 14       treatment:gender:hour 4.27, 19.23  2.48      0.59 .009   .68
## 15                  phase:hour 3.48, 31.36  2.83      0.99  .01   .42
## 16        treatment:phase:hour 6.97, 31.36  2.83      0.33 .010   .93
## 17           gender:phase:hour 3.48, 31.36  2.83      0.90  .01   .47
## 18              age:phase:hour 3.48, 31.36  2.83      0.77  .01   .54
## 19 treatment:gender:phase:hour 6.97, 31.36  2.83      0.65  .02   .71

# replicating ?Anova using aov.car:
aov.car(value ~ treatment * gender + Error(id/phase*hour), data = obk.long, type = 2, return = "Anova")
# in contrast to aov you do not need the within-subject factors outside Error()

# replicating ?Anova using ez.glm:
ez.glm("id", "value", obk.long, c("treatment", "gender"), c("phase", "hour"), type = 2, return = "Anova")

#both return:
## Type II Repeated Measures MANOVA Tests: Pillai test statistic
##                             Df test stat approx F num Df den Df       Pr(>F)    
## (Intercept)                  1     0.970      318      1     10 0.0000000065 ***
## treatment                    2     0.481        5      2     10      0.03769 *  
## gender                       1     0.204        3      1     10      0.14097    
## treatment:gender             2     0.364        3      2     10      0.10447    
## phase                        1     0.851       26      2      9      0.00019 ***
## treatment:phase              2     0.685        3      4     20      0.06674 .  
## gender:phase                 1     0.043        0      2      9      0.82000    
## treatment:gender:phase       2     0.311        1      4     20      0.47215    
## hour                         1     0.935       25      4      7      0.00030 ***
## treatment:hour               2     0.301        0      8     16      0.92952    
## gender:hour                  1     0.293        1      4      7      0.60237    
## treatment:gender:hour        2     0.570        1      8     16      0.61319    
## phase:hour                   1     0.550        0      8      3      0.83245    
## treatment:phase:hour         2     0.664        0     16      8      0.99144    
## gender:phase:hour            1     0.695        1      8      3      0.62021    
## treatment:gender:phase:hour  2     0.793        0     16      8      0.97237    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 


# aggregating over one within-subjects factor (phase) with warning:

aov.car(value ~ treatment * gender + age + Error(id/hour), data = obk.long, observed = c("gender", "age"))

ez.glm("id", "value", obk.long, c("treatment", "gender"), "hour", "age", observed = c("gender", "age"))


# runs with "numeric" factors
obk.long$hour2 <- as.numeric(as.character(obk.long$hour))

aov.car(value ~ treatment * gender + Error(id/hour2), data = obk.long, type = 2,observed = c("gender"))

# only between
aov.car(value ~ treatment * gender + age + Error(id), data = obk.long, type = 2,observed = c("gender", "age"))
aov.car(value ~ treatment * gender + Error(id), data = obk.long, type = 2, observed = c("gender"))

ez.glm("id", "value", obk.long, c("treatment", "gender"), within = NULL, covariate = "age", type = 2, print.formula = TRUE, observed = c("gender", "age"))

ez.glm("id", "value", obk.long, c("treatment", "gender"), within = NULL, type = 2, print.formula = TRUE, observed = c("gender"))

# only within

aov.car(value ~ Error(id/phase*hour), data = obk.long, type = 2)

ez.glm("id", "value", obk.long,  NULL, c("phase", "hour"), type = 2, print.formula = TRUE)

# using return = "full":

str(aov.car(value ~ Error(id/phase*hour), data = obk.long, return = "full"), 1)

## List of 4
##  $ Anova:List of 14
##   ..- attr(*, "class")= chr "Anova.mlm"
##  $ lm   :List of 11
##   ..- attr(*, "class")= chr [1:2] "mlm" "lm"
##  $ data :'data.frame':  16 obs. of  16 variables:
##  $ idata:'data.frame':  15 obs. of  2 variables:

# use args.return arguments:
aov.car(value ~ treatment * gender + age + Error(id/phase*hour), data = obk.long, args.return = list(correction = "none", es = "pes"))

aov.car(value ~ treatment * gender + age + Error(id/phase*hour), data = obk.long,observed = c("gender", "age"), args.return = list(correction = "none", MSE = FALSE))

Run the code above in your browser using DataLab