A robust heteroscedastic procedure using trimmed means.
Usage
trim.test(Y, X, tr = 0.2)
Arguments
Y
A vector of responses. A quantitative vector
X
A vector of factor levels.
tr
The degree of trimming. A value from 0-0.5.
Value
Returns a dataframe with numerator and denominator degrees of freedom, a test statistic, and a p-value based on the F-distribution.
Details
The method utilized here is based on the simple idea of replacing means with trimmed means and standard error estimates, based on all the data, with the standard error of the trimmed mean (Wilcox 2005). The method has the additional benefit of being resistant to heteroscedasticity due to the use of the Welch method for calculating degrees of freedom. With no trimming the degrees of freedom reduce to those of the one way Welch procedure in oneway.test.
References
Wilcox, R. R. 2005. Introduction to Robust Estimation and Hypothesis Testing, Second
Edition. Elsevier, Burlington, MA.