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.
array complexity for the MH test
5 messages · francogrex, David Winsemius, Richard M. Heiberger
An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20120620/0b36d483/attachment.pl>
3 days later
Thanks for your answer. The answer advertises your VCD package, which by the way is a very nice package that I use and recommend for everyone doing such kind of data analysis. However if you really examine the answer you gave me, it does not really or specifically answer my question. -- View this message in context: http://r.789695.n4.nabble.com/array-complexity-for-the-MH-test-tp4633999p4634318.html Sent from the R help mailing list archive at Nabble.com.
On Jun 24, 2012, at 5:30 AM, francogrex wrote:
Thanks for your answer. The answer advertises your VCD package, which by the way is a very nice package that I use and recommend for everyone doing such kind of data analysis. However if you really examine the answer you gave me, it does not really or specifically answer my question.
That's true but neither does your question offer a reproducible example, nor does you response offer the context of either the original question or the response. The thought I had looking at the original was to wonder why you were not using logistic regression or Poisson regression.
Sent from Nabble.com.
As always >>>>>>>>>
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
<<<<<<<<<<<<
David Winsemius, MD West Hartford, CT
1 day later
An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20120625/53947ee0/attachment.pl>