respiratory {multcomp} | R Documentation |
Sub-group and whole-group analysis of a respiratory therapy drug; data set taken from Westfall et al. (1999, p. 195). Three-way ANOVA with weighted contrasts.
data(respiratory)
This data frame contains the following variables
Active
and
Placebo
.Younger
and Older
.Good
and Poor
.See Westfall et al. (1999, p. 195)
P. H. Westfall, R. D. Tobias, D. Rom, R. D. Wolfinger, Y. Hochberg (1999). Multiple Comparisons and Multiple Tests Using the SAS System. Cary, NC: SAS Institute Inc.
data(respiratory) # compute the contrast matrix in several steps # overall active vs. placebo CA <- c(0, 13, 0, 11, 0, 13, 0, 17, 0) CP <- c(0, 0, 14, 0, 12, 0, 19, 0, 12) CA <- CA/sum(CA) CP <- CP/sum(CP) C1 <- CP-CA # for older subgroup only CAO <- c( 0, 13, 0, 0, 0, 13, 0, 0, 0 ) CPO <- c( 0, 0, 14, 0, 0, 0, 19, 0, 0 ) CAO <- CAO/sum(CAO) CPO <- CPO/sum(CPO) C2 <- CPO - CAO # for younger subgroup only CAY <- c( 0, 0, 0, 11, 0, 0, 0, 17, 0 ) CPY <- c( 0, 0, 0, 0, 12, 0, 0, 0, 12 ) CAY <- CAY/sum(CAY) CPY <- CPY/sum(CPY) C3 <- CPY - CAY # subgroup with inital good health CAG <- c( 0, 13, 0, 11, 0, 0, 0, 0, 0 ) CPG <- c( 0, 0, 14, 0, 12, 0, 0, 0, 0 ) CAG <- CAG/sum(CAG) CPG <- CPG/sum(CPG) C4 <- CPG - CAG # subgroup with inital poor health CAP <- c( 0, 0, 0, 0, 0, 13, 0, 17, 0 ) CPP <- c( 0, 0, 0, 0, 0, 0, 19, 0, 12 ) CAP <- CAP/sum(CAP) CPP <- CPP/sum(CPP) C5 <- CPP - CAP # all 4 subgroup combinations of age and initial health condition C6 <- c( 0, -1, 1, 0, 0, 0, 0, 0, 0 ) C7 <- c( 0, 0, 0, 0, 0, -1, 1, 0, 0 ) C8 <- c( 0, 0, 0, -1, 1, 0, 0, 0, 0 ) C9 <- c( 0, 0, 0, 0, 0, 0, 0, -1, 1 ) # contrast matrix, note: first column is zero and corresponds to the # intercept implicitly given in the formula C <- rbind(C1, C2, C3, C4, C5, C6, C7, C8, C9) # numerate the contrasts colnames(C) <- NULL rownames(C) <- c("Overall", "Older", "Younger", "Good Init", "Poor Init", "Old x Good", "Old x Poor", "Young x Good", "Young x Poor") # remove the intercept (not needed, simtest can deal with contrast matrices # with and without a column of zeros for the intercept C <- C[,-1] summary(simtest(Score ~ Treatment:AgeGroup:InitHealth, data=respiratory, ttype="logical", alternative="greater", cmatrix=C))