haireye <- matrix(c(
688,116,584,188,4,
326,38,241,110,3,
343,84,909,412,26,
98,48,403,681,81
), byrow=T, ncol=5, nrow=4,
dimnames = list(c('Light','Blue','medium.eye','Dark.eye'),
                               c('fair','red','medium.hair','dark.hair','black')))