similar to: Mantel-Haenszel relative risk with Greenland-Robins variance estimate

Hi all, Does anyone know if an R function for the Mantel-Haenszel chi-square for a 2x2 table exists? I’ve also seen it called the randomization Q statistic. Note that I’m not looking for the Cochran-Mantel-Haenszel…I did see that out there as cmh.test. Thanks in advance, Jeff Internal Virus Database is out-of-date. Checked by AVG Free Edition. 11:27 AM [[alternative HTML version

Hello, I have some satellite tag time-at-depth (TAD) frequency data that I would like some help with. The data was transmitted via satellite as percent time spent in each of 7 depth bins (0m, 0-1m, 1-10m, 10-50m etc.), binned over 6-hour intervals. I categorized each row of data corresponding to a date and time into summer vs. winter, and day vs. night, and then summed and averaged the given

In SAS, for a two-way (or 3-way, stratified) table, the CMH option in SAS PROC FREQ gives 3 tests that take ordinality of the factors into account, for both variables, just the column variable or neither. Is there an equivalent in R? The mantelhaen.test in stats gives something quite different (a test of conditional independence for *nominal* factors in a 3-way table). e.g. I'd like to

Dear R users, In metafor 1.6-0, the Mantel-Haenszel method is implemented by the rma.mh() function. I have observed that the sum of the weights computed by weights(x) doesn't add to 100% when x is an object of class rma.mh. The consequences of this fact can be clearly seen when a forest diagram is drawn with forest(x), which calls weights(x) (or more precisely, the method weights.rma.mh()

Dear List, I am looking for what B.S.Everitt refers to as Cochrane Method for testing independence in combined 2x2 contingency tables. Is it the same method as the Cochran-Mantel-Haenszel Chi-Squared Test for Count Data in R? Thanks, Serguei [[alternative HTML version deleted]]

Hello, I am searching for an R-Package that does an exentsion of the Mantel-Haenszel test for ordinal data as described in Liu and Agresti (1996) "A Mantel-Haenszel type inference for cummulative odds ratios". in Biometrics. I see packages such as Epi that perform it for binary data and derives a varaince for it using the Robbins and Breslow variance method. As well as another pacakge

Dear R'ers, When I used S-plus i wrote a small program for a Mantel-Haenszel test for trend (I think it worked). Unfortunately I can't get it working in R. It appears as if my use of 'el' is the problem but I can't sort it out. Error in apply(array, c(, 2, 3), function(el) el * 1:s) : argument is missing, with no default Further down in the program I use 'el'

Dear list I would like to compute a Mantel-Haenszel chi-square in which the matching variable is a continuous variable. The MH chi-square is used to assess the relationship between two categorical variables at each level or strata defined by a third variable. Specifically I would like to know if there is a straightforward way to divide the matching variable into levels, in which each level has a

Hi there, is there any way to compare 2 odds ratios? I have two tests that are supposed to detect a disease presence. So for each test, I can compute an odds ratio. My problem is how can I compare the 2 tests by testing whether the 2 odds ratios are the same? Appreciate

Has anyone written R code for extended Mantel-Haenszel Tests- as described in Section 8.4.5 (eq. 8.18) of Agresti's book "Categorical Data Analysis"? In other words, I would like to use column scores to test for a location shift for different groups as defined by the row variable in a set of rxs tables. In addition, I would like to use row and column scores to test for a linear

Dear R users, I am happy to announce the new package `samplesizeCMH` is now available on CRAN <https://cran.r-project.org/package=samplesizeCMH>. The `samplesizeCMH` package will compute power and sample size for the Cochran-Mantel-Haenszel test for stratified 2x2 tables. There are also several helper functions for working with 2x2 tables, such as converting a proportion into an odds

Dear R users, I am happy to announce the new package `samplesizeCMH` is now available on CRAN <https://cran.r-project.org/package=samplesizeCMH>. The `samplesizeCMH` package will compute power and sample size for the Cochran-Mantel-Haenszel test for stratified 2x2 tables. There are also several helper functions for working with 2x2 tables, such as converting a proportion into an odds

Hello, I've tried to analyze some data with a CMH test. My 3 dimensional contingency tables are 2x2xN where N is usually between 10 and 100. The problem is that there may be 2 strata with opposite counts (the 2x2 contigency table for these are reversed), producing opposite odds ratios that cancle out in the overall statistics. These opposite counts are very important for my analysis, since

Dear R users: In the example of meta-analysis (cochrane, package rmeta), I can not found the p-value of Test for overall effect, and some other indices (Z, I, weight and et al). How can I get the these indices listed? > library(rmeta) > data(cochrane) > cochrane name ev.trt n.trt ev.ctrl n.ctrl 1 Auckland 36 532 60 538 2 Block 1 69 5

I have a print method for a set of statistical tests, vcdExtra::CMHtest, for which I'd like to have more sensible printing of pvalues, as in print.anova(). [Testing this requires the latest version of vcdExtra, from R-Forge **|install.packages("vcdExtra", repos="http://R-Forge.R-project.org")|** ] With my current print method, I get results like this, but all Prob values

Hi R users I like to extract (or collect) a numeric element with a name from a list. Is there any way to extract just a numeric element without the name attached to the element. For example, >mylist Mantel-Haenszel chi-squared test with continuity correction data: table(mydata[, x]) Mantel-Haenszel X-squared = 8.3832, df = 1, p-value = 0.003787 alternative hypothesis: true

Hi all: For stratified count data,how to perform regression analysis? My data: age case oc count 1 1 1 21 1 1 2 26 1 2 1 17 1 2 2 59 2 1 1 18 2 1 2 88 2 2 1 7 2 2 2 95 age: 1:<40y 2:>40y case: 1:patient 2:health oc: 1:use drug 2:not use drug My purpose: Anaysis whether case and

Full_Name: Chien-yu Peng Version: 2.0.1 OS: Windows XP Professional Submission from: (NULL) (140.109.72.181) Dear all: Although I don't know you, I am thankful for your help. When I use the function mantelhaen.test for R x C x K (R, C > 2) table, the output is not the same as SAS's. I don't know that the result consist with one of SAS's. But it works correctly for 2

I would like to perform partial mantel tests on only within group values, with "between group" values assigned to NA. This is possible in package ncf partial.mantel.test, however this sues a different permutation to that used in ecodist.ecodist will not accept data with NA values, returning a "matrix is not square error. is it possible to perform this test in ecodist? many thanks

Dear R users, I am having a problem performing the mantel test (both with functions mantel {vegan} and mantel.test {ape}) due to I believe the presence of NAs on my distance matrices, which look like e.g.: NA 1 2 1 2 3 NA 4 5 6 and 1 1 2 2 3 4 5 6 7 8 Would any of you have a solution for that? Thank you for much for your help! Marcia Rocha -- View this message in