Skip to content

array complexity for the MH test

5 messages · francogrex, David Winsemius, Richard M. Heiberger

#
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.
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:

            
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.
As always  >>>>>>>>>
<<<<<<<<<<<<
1 day later