Chi Square for Independence-Mantel–Haenszel test in R

Chi-Square for Independence, It is simple to evaluate using the chi-square test or Fisher exact test when we have two categorical variables.

Consider the case when we have three categorical variables and wish to discover how they are related.

The Cochran–Mantel–Haenszel test, which is an extension of the chi-square test of association, is useful in this scenario.

What will the Cochran–Mantel–Haenszel test reveal?

The Cochran–Mantel–Haenszel test assesses the significance of categorical variables’ associations.

We’ll use the VCD package and the mantelhean.test() function in this example.

Repeated Measures of ANOVA in R Complete Tutorial »

Hypothesis

Let’s look at the test hypotheses.

Ho:-The two inner variables do not have any relationship.

H1: There is a link between the two inner variables.

Mantel–Haenszel test in R

First, we need to install vcd package from cran library,

#install.packages("vcd")

Let’s load the library

library(vcd)
head(Arthritis)
  ID Treatment  Sex Age Improved
1 57   Treated Male  27     Some
2 46   Treated Male  29     None
3 77   Treated Male  30     None
4 17   Treated Male  32   Marked
5 36   Treated Male  46   Marked
6 23   Treated Male  58   Marked

Treatment, Sex, Age, and Improved are the four factors here. The goal is to figure out how these factors are related.

Deep Neural Network in R » Keras & Tensor Flow

Chi Square for Independence

Consider the following three variables: Treatment, Sex, and Improved.

mantelhaen.test(Arthritis$Treatment,Arthritis$Improved,Arthritis$Sex)
        Cochran-Mantel-Haenszel test
data:  Arthritis$Treatment and Arthritis$Improved and Arthritis$Sex
Cochran-Mantel-Haenszel M^2 = 14.632, df = 2, p-value = 0.0006647

Conclusion

Because the p-value of 0.0006647 is less than 0.05, the null hypothesis is rejected and the alternate hypothesis is accepted.

This suggests that at each level of sex, the treatment received and the improvement claimed are not independent.

Similarly, we can divide the age column into two or more components and measure the relationship.

Principal component analysis (PCA) in R »

You may also like...

Leave a Reply

Your email address will not be published. Required fields are marked *

two − two =