R help - Dec 2006 - organizing stats from a list of models

Hi there,
   
  I have a list of models (about 600 glm models) and I included the prefix
"mod_" on their name. Now I would like retrieve the list from the R
environment and save their AICs (and other info) on a table. I´m trying
something like:
   
  -------------------------
  y<-runif(20,min=0,max=1)
  x1<-runif(20,min=0,max=1)
  x2<-runif(20,min=0,max=1)
   
  mymod_1 <-glm(y~x1)
  mymod_2<-glm(y~x2)
  mymod_3<-glm(y~x1+x2)
   
  model.list<-ls(pat="mymod_")
  model.list
  [1] "mymod_1" "mymod_2" "mymod_3"
-------------------------
  Now I would like have something like
   
  model_name call AIC
  mymod_1 glm(y~x1) 13.11942
mymod_2 glm(y~x2) 13.11942
mymod_3 glm(y~x1+x2) 13.11942

  How can I do that?
   
  All the best,
   
  Miltinho
   

 		
---------------------------------

	[[alternative HTML version deleted]]

>   I have a list of models (about 600 glm models) and I included the prefix
"mod_" on their name. Now I would like retrieve the list from the R
environment and save their AICs (and other info) on a table. I?m trying
something like:
I have some (rather un)tested code:

findmodels <- function(modeltype = "lm", dataset, pattern) {
	ls <- ls(".GlobalEnv", pattern=pattern)
	mods <- ls[sapply(ls, function(x) inherits(get(x), modeltype))]
	if (!missing(dataset)) {
		data.name <- function(x) as.character(x$call[["data"]])
		mods <- mods[sapply(mods, function(x) data.name == dataset)]
	}
	
	models <- lapply(mods, get)
	class(models) <- c("ensemble", class(models))
	models
}

summary.ensemble <- function(object, ...) {
	fits <- data.frame(t(sapply(object, function(mod) {
		sum <- summary(mod)
		c(
			df = .df(mod),
			logL = logLik(mod),
			AIC = -AIC(mod),
			BIC = -AIC(mod, k=log(length(fitted(mod)))),
			R2 = sum$r.squared,
			adjR2 = sum$adj.r.squared
		)
	})))
	fits$model <- factor(names(object))
	rownames(fits) <- paste("m", fits$model, sep="")
	fits
}

Hopefully that should give you an idea on how to proceed.

Regards,

Hadley