Skip to content
Snippets Groups Projects
Normal view History Permalink
create-categorical-data.R 1.07 KiB
Newer Older
Matthew K Defenderfer's avatar
Push to git
Matthew K Defenderfer committed
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 
library(gmodels)

good = FALSE
while(!good) {
bacteria <- data.frame(drug.tx = factor(c(rep('Placebo',20),
                                          rep('Treatment',30))),
                       diet = factor(rbinom(50,1,0.2),
                                     levels = c(0,1),
                                     labels = c('Veg','Meat')),
                       bac.presence = factor(c(rbinom(20,1,.7),
                                               rbinom(30,1,.2)),
                                             levels = c(0,1),
                                             labels = c('Not Present','Present')))

bac.ct <- xtabs(~ drug.tx+diet+bac.presence, bacteria)
bac.full <- loglm(~ drug.tx*diet*bac.presence, bac.ct, fit = TRUE)
bac.no3 <- update(bac.full, . ~ . - drug.tx:diet:bac.presence, fit = TRUE)
bac.anova <- anova(bac.full, bac.no3)
bac.list <- as.list(bac.anova)

if(bac.list[[14]] < 0.01){ good = TRUE}
}

CrossTable(bacteria$drug.tx, bacteria$bac.presence, fisher = TRUE, chisq = TRUE, expected = TRUE, sresid = TRUE, format = 'SPSS')

newTable <- xtabs(~., data = bacteria)