[R] MICE data analysis with glmulti

[Andras Farkas]

Dear All,

wonder if you have some thoughts on running the with() function (and perhaps including the pool() function to get the results?) in glmulti? In other words, how to run glmulti with a data set that is produced by mice()?

publicly available code:

data <- airquality
data[4:10,3] <- rep(NA,7)
data[1:5,4] <- NA
data <- data[-c(5,6)]
library(mice)
library(glmulti)

the following line will compute the missing data:
tempData <- mice(data,m=5,maxit=50,meth='pmm',seed=500)

and the following 2 lines will run the regression on the mice output and pool the results to establish the final result of interest for the model specified...
modelFit1 <- with(tempData,glm(Temp~ Ozone+Solar.R+Wind))
summary(pool(modelFit1))


with glmulti I am trying to establish the "best" model by evaluating combinations of all predictors and interactions in different models and would like to force the variable "Ozone" into all models with the following code:

glm.redefined = function(formula, data, always="", ...) 
{glm(as.formula(paste(deparse(formula), always)), data=data, ...)}

then run glmulti:


output<-glmulti(with(tempData,Temp~Solar.R+Wind), 
                fitfunc=glm.redefined, 
                level=1, 
                crit=aic, 
                method="h", 
                always= "+Ozone")


which will obviously fail once you give it a try... any thoughts on how to identify the best model using glmulti in this fashion  that would fit the different combination of predictors with interactions on the mice() output of tempData?

much appreciate the help...

Andras