R

Copy code

# Create contingency tables for Never Smokers and Current Smokers Groups never_smokers <- matrix(c(150, 71, 360348, 117804), nrow = 2, byrow = TRUE) colnames(never_smokers) <- c("# Cases", "Person-Years") rownames(never_smokers) <- c("Female", "Male") current_smokers <- matrix(c(38, 10, 147476, 17706), nrow = 2, byrow = TRUE) colnames(current_smokers) <- c("# Cases", "Person-Years") rownames(current_smokers) <- c("Female", "Male")Suppose we want to study the effect of cigarette smoking on the incidence rate of primary open-angle glaucoma. The following data were obtained from the Nurses' Health Study (NHS) over the time period 1980-1996 and the Health Professionals' Follow-Up Study (a study of male health professionals) over the time period 1986-1996. ## Part a (6 points) If we assume the age distributions of current smokers and never smokers are comparable within each study, then test the hypothesis that current cigarette smoking is associated with the incidence of open-angle glaucoma. Report a two-tailed $p$-value and a conclusion in the context of the problem.