R help - Jun 2009 - Restrict AIC comparison to succesful models?

Hello list,

I'm doing a bootstrap analysis where some models occasionally fail to
converge. I'd like to automate the process of restricting AIC to the
models that do converge. A contrived example of what I'd like to do is
below:

resp <- c(1,1,2)
pred <- c(1,2,3)

m1 <- lm(resp~pred)
m2 <- lm(resp~poly(pred,2))
m3 <- lm(resp~poly(pred,3)) # Fails, obviously

## Some test to see which models were successful
models <- test(m1,m2,m3)

## AIC on these models only
AIC(models)

-- 
http://mutualism.williams.edu

Manuel Morales wrote:
> 
> Hello list,
> 
> I'm doing a bootstrap analysis where some models occasionally fail to
> converge. I'd like to automate the process of restricting AIC to the
> models that do converge. A contrived example of what I'd like to do is
> below:
> 
> resp <- c(1,1,2)
> pred <- c(1,2,3)
> 
> m1 <- lm(resp~pred)
> m2 <- lm(resp~poly(pred,2))
> m3 <- lm(resp~poly(pred,3)) # Fails, obviously
> 
> ## Some test to see which models were successful
> models <- test(m1,m2,m3)
> 
    How about 

models <- list(m1,m2,m3)
is.OK <- sapply(models,test)
do.call(AIC,models[is.OK])

  ?

-- 
View this message in context:
http://www.nabble.com/Restrict-AIC-comparison-to-succesful-models--tp23985797p23987185.html
Sent from the R help mailing list archive at Nabble.com.