# problem-9.3

problem-9.3 The data set gives percentages, but chisq.test() wants proportions. We first define a quick function to compute proportions, then apply the two tests.

    > prop = function(x) x/sum(x)
    > bagfull = c(15,34,7,19,29,24)
    > names(bagfull) = c("blue","brown","green","orange","red","yellow")
    > chisq.test(bagfull,p = prop(mandms["milk chocolate",]))

            Chi-squared test for given probabilities

    data:  bagfull
    X-squared = 7.065, df = 5, p-value = 0.2158

    > chisq.test(bagfull,p = prop(mandms["Peanut",]))

            Chi-squared test for given probabilities

    data:  bagfull
    X-squared = 13.33, df = 5, p-value = 0.02049

It appears that the bag is milk chocolate not peanut.