counts <- c(18,17,15,20,10,20,25,13,12)
outcome <- gl(3,1,9)
treatment <- gl(3,3)
model1 <- glm(counts ~ outcome + treatment, family = poisson())
model2 <- glm(counts ~ outcome, family = poisson())
print_params <- function(model){
cat("\nAIC:", model$aic,
"\nDeviance:", model$deviance)}
list_models <- list(model1,model2)
do.call("rbind", sapply(list_models, print_params, simplify = FALSE))
AIC: 56.76132
Deviance: 5.129141
AIC: 52.76132
Deviance: 5.129141NULL
我不知道爲什麼我得到這個。有沒有另外一種方法呢?
'print_params'回報'NULL':什麼是'的的rbind'它的結果呢? – nicola