MICE data analysis with glmulti

Previous Topic Next Topic
 
classic Classic list List threaded Threaded
1 message Options
Reply | Threaded
Open this post in threaded view
|

MICE data analysis with glmulti

R help mailing list-2
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 

______________________________________________
[hidden email] mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.