> > 3) Is there a way to assign the numbers to the treatment group that they
> were
> > originally in? For example, "1" in the result should be "3", indicating that
> > that was the temperature the experiment took place at. Pairwise.prop.test
> > doesn't seem to accept matrices with more than 2 columns.
>
> No. The test is to do the pairwise comparisons between the samples.
A possible wrinkle, there: If you know the treatment labels for your group
numbers (1...n in a simple htest p-value matrix), you could relabel the htest
p-value matrix via its dimnames property. For example, in the smokers example
under ?pairwise.prop.test:
smk.pwtst <- pairwise.prop.test(smokers, patients)
groupnames <- paste("group", 1:4) #This would have to be manual for your example
new.dims <- list(groupnames[1:3], #we know this is from 1:(n-1)
groupnames[2:4]) #and this is 2:n
dimnames(smk.pwtst$p.value) <- new.dims
smk.pwtst
# group 1 group 2 group 3
# group 2 1.000 - -
# group 3 1.000 1.000 -
# group 4 0.119 0.093 0.124
Steve Ellison
*******************************************************************
This email and any attachments are confidential. Any use...{{dropped:8}}
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
