1 / 35

Comparing Groups – Part 2

Comparing Groups – Part 2. Wilcoxon Signed Rank (again). H 0 : θ = 0 H A : θ ≠ 0. θ represents the population mean or median. or in colors… add up the rank_AbsDiff for the green negative scores.

kimo
Download Presentation

Comparing Groups – Part 2

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. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Comparing Groups – Part 2

  2. Wilcoxon Signed Rank (again)

  3. H0: θ = 0 HA: θ≠ 0 θ represents the population mean or median

  4. or in colors… add up the rank_AbsDiff for the green negative scores or in colors… add up the rank_AbsDiff for the blue positive scores

  5. or in colors… add up the rank_AbsDiff for the green negative scores R- = 1.5 + 3 + 4.5 + 4.5+ 7.5 + 7.5+7.5+14.5+18.5 = 64.5 or in colors… add up the rank_AbsDiff for the blue positive scores R+ = 1.5 + 7.5 + 11 +11 +11 +14.5 +14.5 +14.5+18.5+18.5+18.5+21+22= 188.5 188.5-64.5/2=62

  6. One-way ANOVA in R • In the last lecture I gave you a demonstration of using SAS to conduct a one-way ANOVA. Many people think that R is a superior tool for doing analyses like ANOVA. While R is not critical for simple models like you have seen, it is invaluable for graphics to describe more complicated models.

  7. R and Anxiety Disorders • Start R and then type library(Rcmdr) • Load the dataset using the Import Data option on the Data menu. • Load the Generalized Anxiety Disorder (GAD) data. It has Hamilton Rating Scale for Anxiety data for people taking placebo, high or low dose drug.

  8. Summarize the Data

  9. Show me the means. • You want to see if there is a difference in means in the 3rd variable for the levels of the 2nd variable. You can request the mean plot:

  10. The Plot of the Design I really like boxplots to see the variability around these means.

  11. Set Factor Order

  12. Side-by-side Boxplot

  13. Boxplot Reordered

  14. Gad <- na.omit(Gad) plot.design(Gad[, c(2,3)]) attach(Gad) thePlot = (boxplot(HAMA~DOSEGRP, ylab="HAMA", xlab="DOSEGRP")) where = seq(thePlot$n) theMeans = tapply(HAMA, DOSEGRP, mean) points(where, theMeans, col = "red", pch = 18) I want means.

  15. Testing for Differences • Once the data is set with the placebo group as the baseline, it is easy to ask for tests for the differences vs. the baseline level.

  16. General Linear Models • It is important for you to know that ANOVAs are more than stand-alone analyses that are not related to modeling. ANOVA fits into a class of statistics called General Linear Models. They are well implemented in R and Rcmdr.

  17. Notice you now have a working model.

  18. Checking the Model • With the working model you can get lots of summary information, including diagnostic plots.

  19. More Advanced Contrasts • There is a plethora of methods for dealing with multiple comparisons. You can search CRAN for specific methods you see in textbooks. For example, Walker uses Dunnett’s T test for contrasting the placebo vs. the other levels. • I start search CRAN with RSiteSearch(“blah”)

  20. Nasty code…. • I eventually found the method implemented in a package called multcomp: library(multcomp) gad.aov = aov(HAMA ~ DOSEGRP, data = Gad) summary(glht(gad.aov, linfct = mcp(DOSEGRP="Dunnett"))) confint(glht(gad.aov, linfct = mcp(DOSEGRP="Dunnett")))

  21. Two-Way ANOVA • When you have two or more predictors, you want to know if the variables impact the outcome means by themselves and also if they interact. If variables interact, it means that together they do things to the outcome variable beyond what you would expect from looking at each variable alone. For example, smoking and eating too much both hurt longevity but the combination of the two factors may not be as bad as expected or the combination may be especially lethal.

  22. Anemia • People with Cervical (C), Prostate(P) or Colorectal (R) cancers with chemotherapy-induced anemia were treated either with a drug or placebo and the changes in their hemoglobin levels were assessed. The question of interest is, does the drug reduce anemia?

  23. Box Plots of Course • With code or R commander get boxplots for the levels of the predictors:

  24. Design Plots • You have analysis variables in columns 1, 2 and 4. So specify the design plot like this:

  25. Interaction Plots • In addition to the main effects for the predictors, you want to see if the drug behaves differently in the people with the three different cancer groups. Perhaps it helps increase hemoglobin in one group and decreases it another. Type your outcome variable last.

  26. It looks like there are differences between the cancer types and the drug seems to increase the hemoglobin relative to the placebo.

  27. Modeling Anova with a capital A is part of the car package

  28. Results To turn off the * stuff, type this code before you model: options(show.signif.stars=FALSE)

  29. ANOVA as a Model • You can build this ANOVA as a linear model.

  30. used for nesting y ~ A/B y ~ A + A:B y ~ A + B % in % A interaction only main effect and interaction Set interaction limits (A+B+C) ^ 2 is equal to A*B*C - A:B:C Main effect Remove an explanatory effect.

  31. SAS EG example • See the parallel information in the SAS EG project.

More Related