# Chi-square GOF test: pea flower color
color_obs <- c(705,224)
color_probs <- c(.75,.25)
color_test <- chisq.test(color_obs, p=color_probs)
attach(color_test)
color_test
chisqdif <- observed - expected
chisqterms <- residuals^2
details <- cbind(observed,expected,chisqdif,chisqterms)
rownames(details) <- c("red","white")
details
detach(color_test)
# Chi-square GOF test: pea seed shape/color
shape_obs <- c(315,101,108,32)
shape_probs <- c(.5625, .1875, .1875, .0625)
shape_test <- chisq.test(shape_obs, p=shape_probs)
attach(shape_test)
shape_test
chisqdif <- observed - expected
chisqterms <- residuals^2
details <- cbind(observed,expected,chisqdif,chisqterms)
rownames(details) <- c("RY","WY","RG","WG")
details
detach(shape_test)
# Chi-square GOF test: maize leaf type
maize_obs <- c(773, 231, 238, 59)
maize_probs <- c(.5625, .1875, .1875, .0625)
maize_test <- chisq.test(maize_obs, p=maize_probs)
attach(maize_test)
maize_test
chisqdif <- observed - expected
chisqterms <- residuals^2
details <- cbind(observed,expected,chisqdif,chisqterms)
rownames(details) <- c("green","golden","grst","grgost")
details
detach(maize_test)
# Chi-square GOF test: maize leaf type (omitting green-golden-striped)
maize2_obs <- c(773, 231, 238)
maize2_probs <- c(.6, .2, .2)
maize2_test <- chisq.test(maize2_obs, p=maize2_probs)
attach(maize2_test)
maize2_test
chisqdif <- observed - expected
chisqterms <- residuals^2
details <- cbind(observed,expected,chisqdif,chisqterms)
rownames(details) <- c("green","golden","grst")
details
detach(maize2_test)


-------- results ---------------------------------------

> # Chi-square GOF test: pea flower color
> color_obs <- c(705,224)
> color_probs <- c(.75,.25)
> color_test <- chisq.test(color_obs, p=color_probs)
> attach(color_test)
> color_test

	Chi-squared test for given probabilities

data:  color_obs
X-squared = 0.3907, df = 1, p-value = 0.5319

> chisqdif <- observed - expected
> chisqterms <- residuals^2
> details <- cbind(observed,expected,chisqdif,chisqterms)
> rownames(details) <- c("red","white")
> details
      observed expected chisqdif chisqterms
red        705   696.75     8.25 0.09768568
white      224   232.25    -8.25 0.29305705
> detach(color_test)
> # Chi-square GOF test: pea seed shape/color
> shape_obs <- c(315,101,108,32)
> shape_probs <- c(.5625, .1875, .1875, .0625)
> shape_test <- chisq.test(shape_obs, p=shape_probs)
> attach(shape_test)
> shape_test

	Chi-squared test for given probabilities

data:  shape_obs
X-squared = 0.47, df = 3, p-value = 0.9254

> chisqdif <- observed - expected
> chisqterms <- residuals^2
> details <- cbind(observed,expected,chisqdif,chisqterms)
> rownames(details) <- c("RY","WY","RG","WG")
> details
   observed expected chisqdif chisqterms
RY      315   312.75     2.25 0.01618705
WY      101   104.25    -3.25 0.10131894
RG      108   104.25     3.75 0.13489209
WG       32    34.75    -2.75 0.21762590
> detach(shape_test)
> # Chi-square GOF test: maize leaf type
> maize_obs <- c(773, 231, 238, 59)
> maize_probs <- c(.5625, .1875, .1875, .0625)
> maize_test <- chisq.test(maize_obs, p=maize_probs)
> attach(maize_test)
> maize_test

	Chi-squared test for given probabilities

data:  maize_obs
X-squared = 9.2714, df = 3, p-value = 0.02589

> chisqdif <- observed - expected
> chisqterms <- residuals^2
> details <- cbind(observed,expected,chisqdif,chisqterms)
> rownames(details) <- c("green","golden","grst","grgost")
> details
       observed expected chisqdif chisqterms
green       773 731.8125  41.1875  2.3180940
golden      231 243.9375 -12.9375  0.6861549
grst        238 243.9375  -5.9375  0.1445202
grgost       59  81.3125 -22.3125  6.1226460
> detach(maize_test)
> # Chi-square GOF test: maize leaf type (omitting green-golden-striped)
> maize2_obs <- c(773, 231, 238)
> maize2_probs <- c(.6, .2, .2)
> maize2_test <- chisq.test(maize2_obs, p=maize2_probs)
> attach(maize2_test)
> maize2_test

	Chi-squared test for given probabilities

data:  maize2_obs
X-squared = 2.6914, df = 2, p-value = 0.2604

> chisqdif <- observed - expected
> chisqterms <- residuals^2
> details <- cbind(observed,expected,chisqdif,chisqterms)
> rownames(details) <- c("green","golden","grst")
> details
       observed expected chisqdif chisqterms
green       773    745.2     27.8  1.0370907
golden      231    248.4    -17.4  1.2188406
grst        238    248.4    -10.4  0.4354267
> detach(maize2_test)
>