search for: haenszel

Does anyone know of an R function for computing the Greenland-Robins variance for Mantel-Haenszel relative risks? Thanks Frank -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University

...51.6 2.0640000 31 wint night 300 411.4 16.4560000 32 wint night 500 1981.2 79.2480000 I wanted to test whether significant differences existed between depth in summer vs. winter, and day vs. night, controlling first for season and then for time of day. I carried out a Cochran-Mantel-Haenszel test, using Average Frequency (Avrg) as the dependent variable (2x2x8 contingency table). > ct<-xtabs(Avrg~Time+Depth+Season,data=HG03.dat) > mantelhaen.test(ct) Cochran-Mantel-Haenszel test data: ct Cochran-Mantel-Haenszel M^2 = 28.4548, df = 7, p-value = 0.0001818 > c...

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

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 an...

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]]

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

...ay table). e.g. I'd like to reproduce: *-- CMH tests; proc freq data=sexfun order=data; weight count; tables husband * wife / cmh chisq nocol norow; run; The FREQ Procedure Summary Statistics for Husband by Wife Cochran-Mantel-Haenszel Statistics (Based on Table Scores) Statistic Alternative Hypothesis DF Value Prob 1 Nonzero Correlation 1 10.0142 0.0016 2 Row Mean Scores Differ 3 12.5681 0.0057 3 General Association...

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

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 estimate and converting two row/column proportions into an odds ratio estimate. Please see package vignettes for more information <https://cran.r-pro...

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

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 estimate and converting two row/column proportions into an odds ratio estimate. Please see package vignettes for more information <https://cran.r-pro...

...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 should better be reported as something like '<0.0001' . > CMHtest(MSPatients[,,1]) Cochran-Mantel-Haenszel Statistics for New Orleans Neurologist by Winnipeg Neurologist AltHypothesis Chisq Df Prob cor Nonzero correlation 51.424 1 7.4426e-13 cmeans Col mean scores differ 55.393 3 5.6601e-12 rmeans Row mean scores differ 53.631 3 1.3450e-11 general General associ...

...81 11 63 4 Gamsu 14 131 20 137 5 Morrison 3 67 7 59 6 Papageorgiou 1 71 7 75 7 Tauesch 8 56 10 71 > a=meta.MH(n.trt,n.ctrl,ev.trt,ev.ctrl,names=name,data=cochrane) > summary(a) Fixed effects ( Mantel-Haenszel ) meta-analysis Call: meta.MH(ntrt = n.trt, nctrl = n.ctrl, ptrt = ev.trt, pctrl = ev.ctrl, names = name, data = cochrane) ------------------------------------ OR (lower 95% upper) Auckland 0.58 0.38 0.89 Block 0.16 0.02 1.45 Doran 0.25 0....

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'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...

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 common odds ratio is not equal to 1 95 percent confidence interval: 0.2596963 0.7647255 sample estimates: common odds ratio 0.44...

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

...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 oc are correlated, and age is a stratified variable. My solution: 1,Mantel-Haenszel test by using function "mantelhaen.test" 2,loglinear regression by using function glm(count~case*oc,family=poisson).But I don't know how to handle variable "age",which is the stratified variable. Many thanks for your help. My best. [[alternative HTML version deleted]]

...42 14 28 84 if i do: y<- matrix(c(5,6,16, 19,7,6,7,2,5,10,0,1), nc=3) and chisq.test(y) i found that X-squared=20.3176 but the true value is 15.075 (http://www.math.yorku.ca/SCS/Courses/grcat/grc2.html#H2_43:Sample) Now i tried: the Cochran-Mantel-Haenszel Chi-Squared Test for Count Data: x <- array(c(6, 19, 7, 10, 5, 7, 2, 0, 16, 6, 5, 1), dim = c(2, 2, 3), dimnames = list( Active = c("Female", "Male"), Placebo = c("Female", "Male"), Outcome.Level = c(&q...

...package is now available via CRAN, called metafor. The metafor package consists of a collection of functions for conducting meta-analyses in R. Fixed- and random-effects models (with and without moderators) can be fitted via the general linear (mixed-effects) model. For 2x2 table data, the Mantel-Haenszel and Peto's method are also implemented. The package also provides various plot functions (for example, for forest, funnel, and radial plots) and functions for assessing the model fit and obtaining case diagnostics. Best, -- Wolfgang Viechtbauer Department of Methodology and Statistics Unive...