Description Usage Format Examples
Exactly the same data explained in "adhdData" but reshaped as follows.
1 |
A data frame with 120 rows and 4 variables:
whether the child has ADHD or is healty (normal)
.
the stimulus to which the reaction time in this row corresponds.
an integer ID that identifies the subject to which the reaction time corresponds.
reaction time (milliseconds) to of the aforementioned subject to the aforementioned stimulus.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 | # Omnibus test of a mixed between x within subjects model,
# using trimmed means and Winsorized variances
omnibus_trimmed <- welchADF.test(adhdData2, response = "Milliseconds",between.s = "Group",
within.s = "Stimulus", subject = "Subject", contrast = "omnibus", trimming = TRUE)
# The same using S3 method for class formula. The data can be in separate variables in
# the scope of the call, not necessarily in a data.frame
millis <- adhdData2$Milliseconds
gr <- adhdData2$Group
st <- adhdData2$Stimulus
sbj <- adhdData2$Subject
omnibus_trimmed_formula <- welchADF.test(millis ~ gr*st + (st|sbj),
contrast = "omnibus", trimming = TRUE)
summary(omnibus_trimmed_formula)
# Pairwise contrasts of the effects
pairwise_LSM <- welchADF.test(adhdData2, response = "Milliseconds", between.s = "Group",
within.s = "Stimulus", subject = "Subject", contrast = "all.pairwise", effect = "Stimulus")
pairwise_trimmed <- update(pairwise_LSM, trimming = TRUE)
# Bootstrapping to obtain an empirical critical value
## Not run:
pairwise_trimmed_boot <- update(pairwise_trimmed, bootstrap =TRUE, seed = 123456)
summary(pairwise_trimmed_boot)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.