Learn R Programming

WRS2 (version 0.3-1)

rmanova: A heteroscedastic one-way repeated measures ANOVA for trimmed means.

Description

The rmanova function computes a one-way repeated measures ANOVA for the trimmed means. Homoscedasticity assumption not required. Corresponding post hoc tests can be performed using rmmcp.

Usage

rmanova(y, groups, blocks, tr = 0.2)
rmmcp(y, groups, blocks, tr = 0.2)

Arguments

y
a numeric vector of data values (response).
groups
a vector giving the group of the corresponding elements of y.
blocks
a vector giving the block of the corresponding elements of y.
tr
trim level for the mean.

Value

  • rmanova an object of class "t1way" containing:
  • testvalue of the test statistic
  • df1degrees of freedom
  • df2degrees of freedom
  • p.valuep-value
  • callfunction call
  • rmmcp returns an object of class "mcp1" containing:
  • compinference for all pairwise comparisons
  • fnamesnames of the factor levels

References

Wilcox, R. (2012). Introduction to Robust Estimation and Hypothesis Testing (3rd ed.). Elsevier.

See Also

med1way,t1way

Examples

Run this code
head(WineTasting)
rmanova(WineTasting$Taste, WineTasting$Wine, WineTasting$Taster)

## post hoc
rmmcp(WineTasting$Taste, WineTasting$Wine, WineTasting$Taster)

head(bush)
require(reshape)
bushLong <- melt(bush, id.var = "participant", variable_name = "food")
rmanova(bushLong$value, bushLong$food, bushLong$participant)

## post hoc
rmmcp(bushLong$value, bushLong$food, bushLong$participant)

Run the code above in your browser using DataLab