Chris,

Here is an example.  Please ignore the particulars of the models - they were
built only to demonstrate sapply():


library(lme4)

statedata <- data.frame(state.x77, state.region)

#create empty list object
mods <- list()

#run the models, adding each one to the list
mods$m1 <- glm(Income ~ Life.Exp, data=statedata, family =
gaussian(link="log"))
mods$m2 <- update(mods$m1, . ~ . + HS.Grad)
mods$m3 <- update(mods$m2, . ~ . + Population)
mods$m4 <- lmer(Income ~ Life.Exp + (1|state.region), data=statedata)

#collect the model fit info with sapply() and store as data frame
dfmi <- data.frame(
  Formula=sapply(mods, FUN=function(x){ paste(formula(x)[2], formula(x)[1],
formula(x)[-c(1,2)], collapse=" ") } ),
  AIC=sapply(mods, FUN="AIC") )

#write the data frame to a comma-separated file
write.table(dfmi, file="dfmi.txt", col.names=TRUE, sep=",")

Nathan Pellegrin

	[[alternative HTML version deleted]]

