[R] Re strict AIC comparison to succesful models?

[Ben Bolker]

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.