robust between groups factorial and robust rm n.
Download
Skip this Video
Loading SlideShow in 5 Seconds..
Robust Between Groups Factorial and Robust RM PowerPoint Presentation
Download Presentation
Robust Between Groups Factorial and Robust RM

Loading in 2 Seconds...

play fullscreen
1 / 11

Robust Between Groups Factorial and Robust RM - PowerPoint PPT Presentation


  • 89 Views
  • Uploaded on

Robust Between Groups Factorial and Robust RM. Brief examples. Robustitude. As in most statistical situations there would be a more robust method for going about factorial anova

loader
I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.
capcha
Download Presentation

PowerPoint Slideshow about 'Robust Between Groups Factorial and Robust RM' - tasmine


An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript
robustitude
Robustitude
  • As in most statistical situations there would be a more robust method for going about factorial anova
  • Instead of using means, we might prefer trimmed means or medians so as to have tests based on estimates not so heavily influenced by outliers.
  • And, there’s nothing to it.
robustitude1
Robustitude
  • Or using R you type in something along the lines of
    • t2way(A, B, x, grp=c(1:p), tr=.2, alpha=.05)*
  • Again, don’t be afraid to try robust methods as they are often easily implemented with appropriate software.
robust rm
Robust RM
  • One can perform robust procedures for repeated measures designs
  • Consider our previous dataset regarding stress and test taking
  • I changed the values for the first subject to create an outlier case
    • This person was stressed before and after but not the week of the exam

Before During After

36.00 15.00 35.00

34.00 33.00 22.00

26.00 28.00 25.00

30.00 33.00 17.00

35.00 33.00 30.00

16.00 31.00 23.00

27.00 28.00 28.00

28.00 33.00 18.00

28.00 34.00 27.00

25.00 27.00 25.00

compare results
Compare results
  • SPSS output for this data

 Not significant

compare results1
Compare results
  • Wilcox’s libraries provide a simple way to perform an RM analysis for trimmed means
  • TITLE:

Function rmanova

  • USAGE:

rmanova(x, tr, grp)

  • x refers to the dataset, tr the amount of trimming, and grp can be used to test a subset of the variables
r output
R output
  • rmanova(stressdata, tr=.2)
    • 20% off both ends
  • Results (generalized Hunyh-Feldt correction, provides the epsilon used)
  • Trimmed means (non-trimmed in parentheses):
    • 28.8 (28.5)
    • 31.0 (29.5)
    • 25.0 (25.2)
  • F(2,10) = 3.94, p = .05  sig outcome unlike before
  • To perform a bootstrapped F:
    • rmanovab(stressdata, tr=.2, alpha = .05, nboot = 1000)
multiple comparisons
Multiple comparisons
  • Again, the omnibus F in RM is most likely uninteresting, so perhaps we want to test to see what the specifics of the situation are
  • Function rmmcp can be used
    • rmmcp(x, con=0,tr=.2, alpha=.05, dif=F)
    • dif=T uses difference scores, dif=F the marginal means
    • To test a specific contrast, one can add the con function with appropriate weights or specify a matrix of weights
  • The procedure for testing pairwise comparisons uses Rom’s correction (a modified Bonferonni) for multiple comparisons
    • Order p-values in descending order
    • Test largest at .05, if less than .05 reject and all subsequent
    • If not test next at corrected critical. If less than .05 reject and all subsequent. If not continue until rejection or all hypes have been tested
  • But as mentioned previously, one could use an FDR method also that may be more powerful
r output1
R output
  • rmmcp(stressdata, tr=.1, con=c(0,1,-1),alpha=.05, dif=F)
    • Compares ‘during’ and ‘after’ times
  • Results from the contrast
    • t= 2.7, p = .03
  • Once again, do not be afraid to implement robust procedures when necessary