How to Test a Categorical Data With More Than 2 Categories Against Continuous
Often we use the Pearson Correlation Coefficient to calculate the correlation between continuous numerical variables.
However, we must use a different metric to calculate the correlation between categorical variables – that is, variables that take on names or labels such as:
- Marital status (single, married, divorced)
- Smoking status (smoker, non-smoker)
- Eye color (blue, brown, green)
There are three metrics that are commonly used to calculate the correlation between categorical variables:
1. Tetrachoric Correlation: Used to calculate the correlation between binary categorical variables.
2. Polychoric Correlation: Used to calculate the correlation between ordinal categorical variables.
3. Cramer's V: Used to calculate the correlation between nominal categorical variables.
The following sections provide an example of how to calculate each of these three metrics.
Metric 1: Tetrachoric Correlation
Tetrachoric correlation is used to calculate the correlation between binary categorical variables. Recall that binary variables are variables that can only take on one of two possible values.
The value for tetrachoric correlation ranges from -1 to 1 where -1 indicates a strong negative correlation, 0 indicates no correlation, and 1 indicates a strong positive correlation.
For example, suppose want to know whether or not gender is associated with political party preference so we take a simple random sample of 100 voters and survey them on their political party preference.
The following table shows the results of the survey:
We would use tetrachoric correlation in this scenario because each categorical variable is binary – that is, each variable can only take on two possible values.
We can use the following code in R to calculate the tetrachoric correlation between the two variables:
library(psych) #create 2x2 table data = matrix(c(19, 12, 30, 39), nrow=2) #view table data #calculate tetrachoric correlation tetrachoric(data) tetrachoric correlation [1] 0.27
The tetrachoric correlation turns out to be 0.27. This value is fairly low, which indicates that there is a weak association (if any) between gender and political party preference.
Metric 2: Polychoric Correlation
Polychoric correlation is used to calculate the correlation between ordinal categorical variables. Recall that ordinal variables are variables whose possible values have a natural order.
The value for polychoric correlation ranges from -1 to 1 where -1 indicates a strong negative correlation, 0 indicates no correlation, and 1 indicates a strong positive correlation.
For example, suppose want to know whether or not two different movie ratings agencies have a high correlation between their movie ratings.
We ask each agency to rate 20 different movies on a scale of 1 to 3 with 1 indicating "bad", 2 indicating "mediocre", and 3 indicating "good."
The following table shows the results:
We can use the following code in R to calculate the polychoric correlation between the ratings of the two agencies:
library(polycor) #define movie ratings x <- c(1, 1, 2, 2, 3, 2, 2, 3, 2, 3, 3, 2, 1, 2, 2, 1, 1, 1, 2, 2) y <- c(1, 1, 2, 1, 3, 3, 3, 2, 2, 3, 3, 3, 2, 2, 2, 1, 2, 1, 3, 3) #calculate polychoric correlation between ratings polychor(x, y) [1] 0.7828328
The polychoric correlation turns out to be 0.78. This value is quite high, which indicates that there is a strong positive association between the ratings from each agency.
Metric 3: Cramer's V
Cramer's V is used to calculate the correlation between nominal categorical variables. Recall that nominal variables are ones that take on category labels but have no natural ordering.
The value for Cramer's V ranges from 0 to 1, with 0 indicating no association between the variables and 1 indicating a strong association between the variables.
For example, suppose we want to know if there is a correlation between eye color and gender so we survey 50 individuals and obtain the following results:
We can use the following code in R to calculate Cramer's V for these two variables:
library(rcompanion) #create table data = matrix(c(6, 9, 8, 5, 12, 10), nrow=2) #view table data [,1] [,2] [,3] [1,] 6 8 12 [2,] 9 5 10 #calculate Cramer's V cramerV(data) Cramer V 0.1671
Cramer's V turns out to be 0.1671. This value is quite low, which indicates that there is a weak association between gender and eye color.
Additional Resources
Introduction to the Pearson Correlation Coefficient
Introduction to Tetrachoric Correlation
Categorical vs. Quantitative Variables: What's the Difference?
Levels of Measurement: Nominal, Ordinal, Interval and Ratio
Source: https://www.statology.org/correlation-between-categorical-variables/
0 Response to "How to Test a Categorical Data With More Than 2 Categories Against Continuous"
Post a Comment