Module 11

> df <- data.frame("Choose the same hotel again" = c("Yes", "No", "Total"), "Beachcomber" = c(163, 64, 227), "Windsurfer" = c(154, 108, 262), "Total" = c(317, 172, 489))
> chisq.test(df$Beachcomber)

Chi-squared test for given probabilities

data:  df$Beachcomber
X-squared = 89.132, df = 2, p-value <
2.2e-16

> chisq.test(df$Windsurfer)                    

Chi-squared test for given probabilities

data:  df$Windsurfer

X-squared = 71.557, df = 2, p-value =

2.894e-16

> chisq <- sum((df$Beachcomber-df$Windsurfer)^2/df$Windsurfer)

> chisq

[1] 23.12747

> pchisq(c(23.127), df=4, lower.tail = FALSE)

[1] 0.0001194391

>  qplot(Beachcomber, Windsurfer, data = df)

>