You have a very nice graph of a dose-response function here.
library(vcd)
library(RColorBrewer)
Pen <- array(c(0, 0, 6, 5,
3, 0, 3, 6,
6, 2, 0, 4,
5, 6, 1, 0,
2, 5, 0, 0),
dim = c(2, 2, 5),
dimnames = list(
Delay = c("None", "1.5h"),
Response = c("Cured", "Died"),
Penicillin.Level = c("1/8", "1/4", "1/2", "1", "4")))
mosaic(structable(aperm(Pen, c(3, 1, 2)), split_vertical=c(TRUE, TRUE,
FALSE)),
highlighting=3,
highlighting_fill=brewer.pal.ext(2,'RdYlGn',reverse=TRUE))
See our paper in the 2011 UseR conference for more information on the graph.
Heiberger, R. (Presenter & Author), Neuwirth, E. (Presenter & Author),
UserR! 2011, "Graphical Syntax for Structables and their Mosaic Plots,"
University of Warwick and R Foundation for Statistical Computing, Coventry,
England. (August 2011).
http://web.warwick.ac.uk/statsdept/user-2011/TalkSlides/Contributed/16Aug_1600_FocusII_3-Visual_2-Heiberger.pdf
We have a MH graph on pages 15-16.
I don't think MH applies for this penicillin example, since the five dose
levels are not giving the same story. The dose-response relationship is
the story here.
Rich
On Wed, Jun 20, 2012 at 1:05 PM, francogrex <francog...@mail.com> wrote:
> If we take the matel-haenszel test on these data of five 2x2 tables
> stratified along Penicillin.Levels
> array(c(0, 0, 6, 5,
> 3, 0, 3, 6,
> 6, 2, 0, 4,
> 5, 6, 1, 0,
> 2, 5, 0, 0),
> dim = c(2, 2, 5),
> dimnames = list(
> Delay = c("None", "1.5h"),
> Response = c("Cured", "Died"),
> Penicillin.Level = c("1/8", "1/4", "1/2", "1", "4")))
>
> The test works fine. How would one proceed for example if within each
> Penicillin level we have two sub 2x2 tables split along the gender ...?
> I mean practically how would that array be and when you have multiple
> stratification levels how practical it is to work that into that array
> structure above? Can one not have a different data.frame structure that the
> test likes as well, or should we use another package (like the meta
> package)?
>
> --
> View this message in context:
> http://r.789695.n4.nabble.com/array-complexity-for-the-MH-test-tp4633999.html
> Sent from the R help mailing list archive at Nabble.com.
>
> ______________________________________________
> R-help@r-project.org mailing list
> 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.
>
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
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.
