data(forams,package='geostats')
abundant <- forams[,c('quinqueloba','pachyderma','incompta',
'glutinata','bulloides')]
other <- rowSums(forams[,c('uvula','scitula')])
dat <- cbind(abundant,other)
chisq.test(dat)
Run the code above in your browser using DataLab