Data Science Day 15

Risk Ratio

Last time, we give a SAS example of Risk Difference to test if two groups are experiencing the same proportion of a certain event. In order to understand the topic better, we will go over Risk Ratio.

Definition:

Risk Ratio or Relative Risk (RR) is the probability that an event occurs in a group 1 relative to the probability that the same event occurs in group 2.
Note: Without Loss of Generality, In clinical trials, we take a group to be actual treatment group and the other group is the placebo group.

Formula:

Risk Ratio (Relative Risk)= P(event in group1)/P(event in group 2)

Interpretation:
RR ≈ 1 ⇒ The proportion of events are similar in group 1 and group 2.
RR » 1 ⇒ Increased probability of events among those in group 1 compared to group 2.
RR « 1 ⇒ Decreased probability of events among those in group 1 compared to group 2.
Note: This is a general interpretation, we can explain in different scenarios.

Example:

#Import Library
import numpy as np
import pandas as pd
import statsmodels.api as sm
import pylab as pl
import matplotlib.pyplot as plt
import seaborn as sns

#Read in Arthritis data
df = sm.datasets.get_rdataset("Arthritis", "vcd").data
#Create New variable Censor:
#Censor=0 means no-improvement
#Censor=1 means some or big improvement 
df["Censor"]=np.where(df["Improved"]=="None", "0", "1")
#table for Treatment and Censor
tab=pd.crosstab(df["Treatment"], df["Censor"])
tab1=sm.stats.Table(tab)
print(tab1.table_orig)
table= np.asarray([[29,14],[13,28]])
t22=sm.stats.Table2x2(table)
print(t22.summary())

#Output
Censor      0   1
Treatment        
Placebo    29  14
Treated    13  28


 Estimate   SE   LCB   UCB   p-value
--------------------------------------------------
Odds ratio        4.462       1.785 11.154   0.001
Log odds ratio    1.495 0.467 0.579  2.412   0.001
Risk ratio        2.127       0.579  2.412   0.003
Log risk ratio    0.755 0.253 0.260  1.250   0.003
--------------------------------------------------

#Visulatization

%matplotlib inline
pd.crosstab(df.Treatment, df.Censor).plot(kind="bar")
plt.title("Treatment vs Placebo")
plt.xlabel("Treatment Type for Arthritis")
plt.ylabel("Censor Status")

Summary
Censor=0, Arthritis problem remains
Censor=1, Arthritis problem improved
Odds Ratio
Odds Ratio= 4.462, this implies in Placebo Group has 4.5 times more likely to remain the Arthritis problem compared to the Treatment Group.

Risk Ratio
Risk Ratio=2.127 means in Placebo Group the probability of patients’ Arthritis problem remains increased compared to the Patients in Treatment Group.

Connections:

It is very easy to mix the concept for Odds Ratio and Risk Ratio.
Odds Ratio represents the odds of an event in group 1 compared to the event of odds in group 2, where odds means the event over non-event.

$Odds\ in \ Placebo= \frac{29}{14}$

$Odds \ in \ Treatment= \frac{13}{28}$

$Risk \ Ratio = \frac{2.0714}{0.464}=4.46$

Risk Ratio means the probability of an event occurring in group 1 compared to the probability of the same event occurring in group 2.

$P(censor=0 \ in \ Placebo)= \frac{29}{43}$

$P(censor=0 \ in \ Treatment)= \frac{13}{41}$

$Risk \ Ratio = \frac{0.674}{0.317}=2.127$

Conclusion:
Risk ratio compares the probability of the occurrence of the same event in two groups. In addition, we can use the Risk Difference to check if the two groups have the same Risk Ratio.

Happy Studying!