added na.rm = T, because some observations are missing data - NA