A new R package is presented for dealing with non-normality and variance heterogeneity of sample data when conducting hypothesis tests of main effects and interactions in mixed models. The proposal departs from an existing SAS program which implements Johansen’s general formulation of Welch-James’s statistic with approximate degrees of freedom, which makes it suitable for testing any linear hypothesis concerning cell means in univariate and multivariate mixed model designs when the data pose non-normality and non-homogeneous variance. Improved type I error rate control is obtained using bootstrapping for calculating an empirical critical value, whereas robustness against non-normality is achieved through trimmed means and Winsorized variances. A wrapper function eases the application of the test in common situations, such as performing omnibus tests on all effects and interactions, pairwise contrasts, and tetrad contrasts of two-way interactions. The package is demonstrated in several problems including unbalanced univariate and multivariate designs.
Supplementary materials are available in addition to this article. It can be downloaded at RJ-2017-049.zip
ART, WRS2, robustbase, robust, robustlmm, nlme, lme4, welchADF, gamm4, mgcv
Robust, Econometrics, Environmetrics, SocialSciences, Bayesian, OfficialStatistics, Psychometrics, SpatioTemporal, ChemPhys, Finance, Multivariate, Spatial
Text and figures are licensed under Creative Commons Attribution CC BY 4.0. The figures that have been reused from other sources don't fall under this license and can be recognized by a note in their caption: "Figure from ...".
For attribution, please cite this work as
Villacorta, "The welchADF Package for Robust Hypothesis Testing in Unbalanced Multivariate Mixed Models with Heteroscedastic and Non-normal Data", The R Journal, 2017
BibTeX citation
@article{RJ-2017-049,
author = {Villacorta, Pablo J.},
title = {The welchADF Package for Robust Hypothesis Testing in Unbalanced Multivariate Mixed Models with Heteroscedastic and Non-normal Data},
journal = {The R Journal},
year = {2017},
note = {https://doi.org/10.32614/RJ-2017-049},
doi = {10.32614/RJ-2017-049},
volume = {9},
issue = {2},
issn = {2073-4859},
pages = {309-328}
}