>
> Kway <- function(formula, family, data, ..., order=nt) {
>  models <- list()
>  mod <- glm(formula, family, data, ...)
>  terms <- terms(formula)
>  tl <- attr(terms, "term.labels")
>  nt <- length(tl)
>  models[[1]] <- mod      for(i in 2:order) {
>      models[[i]] <- update(mod, .~.^i)
>  }      # null model
>

debug shows the problem is here.  Examining that little section separately,
it seems that update() and the for loop are not working together properly,
but I am not certain how to fix it.


>  mod0 <- update(mod, .~1)
>  models <- c(mod0, models)      class(models) <- "glmlist"
>  models
> }
>
>



-- 
Joshua Wiley
Senior in Psychology
University of California, Riverside
http://www.joshuawiley.com/

	[[alternative HTML version deleted]]

