pdf("ex0508.pdf")

# Read in data from ASCII formatted (Fixed With File)

ex0508dat <- read.fwf("http://www.stat.ufl.edu/~winner/data/biostat/ex0508.dat", width=c(8,8,8), 
        col.names=c("marital","cancer","cancer_r"))


# Create a dataset (frame) from input data

ex0508 <- data.frame(ex0508dat)

# Attach the dataset for analysis 

attach(ex0508)

# Obtain a cross-tabulation of marital status and cancer

freqmc <- table(marital,cancer)

freqmc

# Obtain the Chi-square statistic and p-value

chisq.test(freqmc)






dev.off()