R version 3.2.2 (2015-08-14) -- "Fire Safety" Copyright (C) 2015 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu (64-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. Natural language support but running in an English locale R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > > library(glmbb) > packageVersion("glmbb") [1] ‘0.2’ > options(width = 256) > > foo <- read.table("http://www.stat.umn.edu/geyer/5421/data/table-10.1.txt", + header = TRUE) > sapply(foo, class) alcohol cigarette marijuana race gender count "factor" "factor" "factor" "factor" "factor" "integer" > > gout <- glmbb(count ~ (alcohol + cigarette + marijuana + race + gender)^4, + little = ~ alcohol + cigarette + marijuana + race + gender, + data = foo, graphical = TRUE) There were 19 warnings (use warnings() to see them) > summary(gout) Results of search for graphical models with lowest AIC. Search was for all models with AIC no larger than min(AIC) + 10 These are shown below. criterion weight formula 183.0 0.1392642 count ~ alcohol*cigarette*marijuana + alcohol*cigarette*race + alcohol*marijuana*gender 184.3 0.0735257 count ~ alcohol*cigarette*marijuana + alcohol*cigarette*race + alcohol*marijuana*gender + alcohol*race*gender 184.4 0.0703402 count ~ alcohol*cigarette*marijuana + alcohol*cigarette*race + cigarette*marijuana*gender 184.4 0.0695551 count ~ alcohol*race + alcohol*cigarette*marijuana + alcohol*marijuana*gender 184.5 0.0654098 count ~ marijuana*gender + alcohol*cigarette*marijuana + alcohol*cigarette*race 184.8 0.0560423 count ~ alcohol*cigarette*marijuana + alcohol*marijuana*race + alcohol*marijuana*gender 185.6 0.0381889 count ~ marijuana*gender + race*gender + alcohol*cigarette*marijuana + alcohol*cigarette*race 185.7 0.0369259 count ~ alcohol*cigarette*marijuana + alcohol*marijuana*gender + alcohol*race*gender 185.8 0.0351312 count ~ alcohol*race + alcohol*cigarette*marijuana + cigarette*marijuana*gender 185.9 0.0326688 count ~ alcohol*race + marijuana*gender + alcohol*cigarette*marijuana 186.2 0.0284964 count ~ marijuana*race + alcohol*cigarette*marijuana + alcohol*marijuana*gender 186.2 0.0283061 count ~ alcohol*cigarette*marijuana + alcohol*marijuana*race + cigarette*marijuana*gender 186.3 0.0263220 count ~ marijuana*gender + alcohol*cigarette*marijuana + alcohol*marijuana*race 186.9 0.0195891 count ~ alcohol*race + race*gender + alcohol*cigarette*marijuana + cigarette*marijuana*gender 187.0 0.0186964 count ~ alcohol*race + marijuana*gender + race*gender + alcohol*cigarette*marijuana 187.3 0.0165222 count ~ alcohol*cigarette*race + alcohol*cigarette*marijuana*gender 187.3 0.0164781 count ~ alcohol*cigarette*marijuana + alcohol*cigarette*race + cigarette*marijuana*gender + cigarette*race*gender 187.6 0.0143931 count ~ marijuana*race + alcohol*cigarette*marijuana + cigarette*marijuana*gender 187.7 0.0133842 count ~ marijuana*race + marijuana*gender + alcohol*cigarette*marijuana 187.9 0.0120283 count ~ alcohol*gender + alcohol*cigarette*marijuana + alcohol*cigarette*race 187.9 0.0118436 count ~ gender + alcohol*cigarette*marijuana + alcohol*cigarette*race 188.1 0.0108903 count ~ alcohol*marijuana*gender + alcohol*cigarette*marijuana*race 188.7 0.0082520 count ~ alcohol*race + alcohol*cigarette*marijuana*gender 188.9 0.0074619 count ~ alcohol*cigarette*marijuana*gender + alcohol*cigarette*race*gender 189.1 0.0066672 count ~ alcohol*cigarette*marijuana + alcohol*marijuana*gender + marijuana*race*gender 189.1 0.0066488 count ~ alcohol*marijuana*race + alcohol*cigarette*marijuana*gender 189.2 0.0064625 count ~ race*gender + alcohol*cigarette*marijuana + alcohol*cigarette*race 189.2 0.0064465 count ~ alcohol*cigarette*marijuana + alcohol*marijuana*race + cigarette*marijuana*gender + marijuana*race*gender 189.2 0.0063856 count ~ alcohol*cigarette*marijuana + alcohol*cigarette*race + alcohol*race*gender 189.3 0.0061585 count ~ alcohol*cigarette*marijuana + alcohol*marijuana*race + marijuana*race*gender 189.3 0.0061389 count ~ alcohol*cigarette*marijuana + alcohol*marijuana*race*gender 189.3 0.0060075 count ~ alcohol*race + alcohol*gender + alcohol*cigarette*marijuana 189.3 0.0059152 count ~ gender + alcohol*race + alcohol*cigarette*marijuana 189.5 0.0055005 count ~ cigarette*marijuana*gender + alcohol*cigarette*marijuana*race 189.6 0.0051150 count ~ marijuana*gender + alcohol*cigarette*marijuana*race 189.7 0.0048404 count ~ alcohol*gender + alcohol*cigarette*marijuana + alcohol*marijuana*race 189.8 0.0047661 count ~ gender + alcohol*cigarette*marijuana + alcohol*marijuana*race 189.9 0.0044497 count ~ cigarette*gender + alcohol*cigarette*marijuana + alcohol*cigarette*race 189.9 0.0043809 count ~ alcohol*race*gender + alcohol*cigarette*marijuana*gender 190.1 0.0039443 count ~ alcohol*cigarette*marijuana + alcohol*marijuana*gender + cigarette*marijuana*race 190.5 0.0033808 count ~ marijuana*race + alcohol*cigarette*marijuana*gender 190.5 0.0033675 count ~ alcohol*cigarette*marijuana + cigarette*marijuana*gender + marijuana*race*gender 190.5 0.0032277 count ~ alcohol*race + race*gender + alcohol*cigarette*marijuana 190.6 0.0031893 count ~ alcohol*cigarette*marijuana + alcohol*race*gender 190.6 0.0031315 count ~ alcohol*cigarette*marijuana + marijuana*race*gender 190.9 0.0027279 count ~ cigarette*race + alcohol*cigarette*marijuana + alcohol*marijuana*gender 191.0 0.0026006 count ~ race*gender + alcohol*cigarette*marijuana + alcohol*marijuana*race 191.0 0.0025697 count ~ alcohol*cigarette*marijuana + alcohol*marijuana*race + alcohol*race*gender 191.1 0.0024612 count ~ alcohol*gender + marijuana*race + alcohol*cigarette*marijuana 191.1 0.0024235 count ~ gender + marijuana*race + alcohol*cigarette*marijuana 191.2 0.0023059 count ~ race + alcohol*cigarette*marijuana + alcohol*marijuana*gender 191.3 0.0022224 count ~ alcohol*race + cigarette*gender + alcohol*cigarette*marijuana 191.4 0.0020739 count ~ alcohol*cigarette*marijuana + alcohol*cigarette*race + alcohol*cigarette*gender 191.5 0.0019922 count ~ alcohol*cigarette*marijuana + cigarette*marijuana*race + cigarette*marijuana*gender 191.7 0.0018526 count ~ marijuana*gender + alcohol*cigarette*marijuana + cigarette*marijuana*race 191.7 0.0017906 count ~ cigarette*gender + alcohol*cigarette*marijuana + alcohol*marijuana*race 192.1 0.0015137 count ~ cigarette*race + race*gender + alcohol*cigarette*marijuana + alcohol*marijuana*gender 192.2 0.0013778 count ~ cigarette*race + alcohol*cigarette*marijuana + cigarette*marijuana*gender 192.3 0.0013224 count ~ marijuana*race + race*gender + alcohol*cigarette*marijuana 192.3 0.0013126 count ~ alcohol*gender + marijuana*race + race*gender + alcohol*cigarette*marijuana 192.4 0.0012920 count ~ alcohol*cigarette*marijuana*race + alcohol*cigarette*marijuana*gender 192.4 0.0012812 count ~ cigarette*race + marijuana*gender + alcohol*cigarette*marijuana 192.4 0.0012582 count ~ race*gender + alcohol*cigarette*marijuana + alcohol*marijuana*gender 192.5 0.0012066 count ~ alcohol*race + cigarette*gender + race*gender + alcohol*cigarette*marijuana 192.5 0.0011967 count ~ marijuana*race*gender + alcohol*cigarette*marijuana*race 192.5 0.0011929 count ~ alcohol*cigarette*marijuana*race + alcohol*marijuana*race*gender 192.6 0.0011647 count ~ race + alcohol*cigarette*marijuana + cigarette*marijuana*gender 192.7 0.0010830 count ~ race + marijuana*gender + alcohol*cigarette*marijuana 192.8 0.0010358 count ~ alcohol*race + alcohol*cigarette*marijuana + alcohol*cigarette*gender 192.9 0.0009903 count ~ alcohol*cigarette*marijuana + alcohol*cigarette*race + cigarette*race*gender 192.9 0.0009701 count ~ cigarette*gender + race*gender + alcohol*cigarette*marijuana + alcohol*marijuana*race 193.0 0.0009406 count ~ alcohol*gender + alcohol*cigarette*marijuana*race > > # as expected, AIC "selects" lots of models, none very strongly > > gout.bic <- glmbb(count ~ (alcohol + cigarette + marijuana + race + gender)^4, + little = ~ alcohol + cigarette + marijuana + race + gender, + data = foo, graphical = TRUE, criterion = "BIC") There were 19 warnings (use warnings() to see them) > summary(gout.bic) Results of search for graphical models with lowest BIC. Search was for all models with BIC no larger than min(BIC) + 10 These are shown below. criterion weight formula 203.5 0.164700 count ~ alcohol*race + marijuana*gender + alcohol*cigarette*marijuana 204.9 0.080971 count ~ alcohol*race + alcohol*cigarette*marijuana + alcohol*marijuana*gender 205.0 0.076145 count ~ marijuana*gender + alcohol*cigarette*marijuana + alcohol*cigarette*race 205.3 0.067477 count ~ marijuana*race + marijuana*gender + alcohol*cigarette*marijuana 205.5 0.062060 count ~ gender + alcohol*race + alcohol*cigarette*marijuana 206.1 0.045294 count ~ alcohol*race + marijuana*gender + race*gender + alcohol*cigarette*marijuana 206.3 0.040897 count ~ alcohol*race + alcohol*cigarette*marijuana + cigarette*marijuana*gender 206.5 0.037435 count ~ alcohol*cigarette*marijuana + alcohol*cigarette*race + alcohol*marijuana*gender 206.7 0.033173 count ~ marijuana*race + alcohol*cigarette*marijuana + alcohol*marijuana*gender 206.9 0.030642 count ~ marijuana*gender + alcohol*cigarette*marijuana + alcohol*marijuana*race 206.9 0.030287 count ~ alcohol*race + alcohol*gender + alcohol*cigarette*marijuana 207.0 0.028692 count ~ gender + alcohol*cigarette*marijuana + alcohol*cigarette*race 207.2 0.025426 count ~ gender + marijuana*race + alcohol*cigarette*marijuana 207.6 0.021363 count ~ marijuana*gender + race*gender + alcohol*cigarette*marijuana + alcohol*cigarette*race 207.8 0.018908 count ~ alcohol*cigarette*marijuana + alcohol*cigarette*race + cigarette*marijuana*gender 208.1 0.016755 count ~ marijuana*race + alcohol*cigarette*marijuana + cigarette*marijuana*gender 208.1 0.016272 count ~ alcohol*race + race*gender + alcohol*cigarette*marijuana 208.3 0.015065 count ~ alcohol*cigarette*marijuana + alcohol*marijuana*race + alcohol*marijuana*gender 208.4 0.014002 count ~ alcohol*gender + alcohol*cigarette*marijuana + alcohol*cigarette*race 208.7 0.012408 count ~ alcohol*gender + marijuana*race + alcohol*cigarette*marijuana 208.8 0.011546 count ~ gender + alcohol*cigarette*marijuana + alcohol*marijuana*race 208.9 0.011363 count ~ race + marijuana*gender + alcohol*cigarette*marijuana 208.9 0.011204 count ~ alcohol*race + cigarette*gender + alcohol*cigarette*marijuana 208.9 0.010958 count ~ alcohol*race + race*gender + alcohol*cigarette*marijuana + cigarette*marijuana*gender 209.1 0.009926 count ~ alcohol*cigarette*marijuana + alcohol*marijuana*gender + alcohol*race*gender 209.7 0.007609 count ~ alcohol*cigarette*marijuana + alcohol*marijuana*race + cigarette*marijuana*gender 209.7 0.007523 count ~ race*gender + alcohol*cigarette*marijuana + alcohol*cigarette*race 209.9 0.006667 count ~ marijuana*race + race*gender + alcohol*cigarette*marijuana 210.0 0.006459 count ~ cigarette*race + marijuana*gender + alcohol*cigarette*marijuana 210.3 0.005635 count ~ alcohol*gender + alcohol*cigarette*marijuana + alcohol*marijuana*race 210.3 0.005586 count ~ race + alcohol*cigarette*marijuana + alcohol*marijuana*gender 210.4 0.005180 count ~ cigarette*gender + alcohol*cigarette*marijuana + alcohol*cigarette*race 210.7 0.004590 count ~ cigarette*gender + marijuana*race + alcohol*cigarette*marijuana 210.7 0.004564 count ~ alcohol*cigarette*marijuana + alcohol*cigarette*race + alcohol*marijuana*gender + alcohol*race*gender 210.8 0.004282 count ~ race + gender + alcohol*cigarette*marijuana 211.1 0.003713 count ~ alcohol*cigarette*marijuana + alcohol*race*gender 211.1 0.003645 count ~ alcohol*cigarette*marijuana + marijuana*race*gender 211.4 0.003180 count ~ alcohol*gender + marijuana*race + race*gender + alcohol*cigarette*marijuana 211.4 0.003176 count ~ cigarette*race + alcohol*cigarette*marijuana + alcohol*marijuana*gender 211.5 0.003027 count ~ race*gender + alcohol*cigarette*marijuana + alcohol*marijuana*race 211.5 0.002979 count ~ marijuana*gender + race*gender + alcohol*cigarette*marijuana 211.6 0.002923 count ~ alcohol*race + cigarette*gender + race*gender + alcohol*cigarette*marijuana 211.6 0.002822 count ~ race + alcohol*cigarette*marijuana + cigarette*marijuana*gender 211.9 0.002434 count ~ gender + cigarette*race + alcohol*cigarette*marijuana 212.2 0.002157 count ~ marijuana*gender + alcohol*cigarette*marijuana + cigarette*marijuana*race 212.2 0.002090 count ~ race + alcohol*gender + alcohol*cigarette*marijuana 212.2 0.002085 count ~ cigarette*gender + alcohol*cigarette*marijuana + alcohol*marijuana*race 212.5 0.001792 count ~ alcohol*cigarette*marijuana + alcohol*marijuana*gender + marijuana*race*gender 212.6 0.001757 count ~ cigarette*race + marijuana*gender + race*gender + alcohol*cigarette*marijuana 212.6 0.001716 count ~ alcohol*cigarette*marijuana + alcohol*cigarette*race + alcohol*race*gender 212.7 0.001655 count ~ alcohol*cigarette*marijuana + alcohol*marijuana*race + marijuana*race*gender 212.8 0.001604 count ~ cigarette*race + alcohol*cigarette*marijuana + cigarette*marijuana*gender 213.0 0.001465 count ~ race*gender + alcohol*cigarette*marijuana + alcohol*marijuana*gender 213.3 0.001206 count ~ alcohol*race + alcohol*cigarette*marijuana + alcohol*cigarette*gender 213.4 0.001198 count ~ cigarette*gender + marijuana*race + race*gender + alcohol*cigarette*marijuana 213.4 0.001188 count ~ alcohol*gender + cigarette*race + alcohol*cigarette*marijuana 213.5 0.001123 count ~ race*gender + alcohol*cigarette*marijuana > > > proc.time() user system elapsed 146.741 0.136 146.983