Function to calculate Cohen's D value when testing effect size

cohens_d(x, y, na.rm = TRUE)