Как получить гофстат из подогнанных моделей по группам в R

У меня есть большой набор данных, состоящий из тысяч уникальных тестов (TestNum) и связанных с ними ответов (Response), причем ответы различаются по длине в разных тестах. Тесты отбрасываются по некоторым критериям, отсюда и отсутствие последовательности в значениях TestNum. Вот упрощенный пример:

dat <- data.frame(Response=c(rlnorm(10, 2.9, 0.3), rlnorm(14, 2.88, 0.38), rlnorm(19, 2.44, 0.08)),TestNum=rep(c(1,4,9), times=c(10,14,19)))
dat$TestNum<-factor(dat$TestNum)
dat

Я подгоняю распределение lnorm к каждому TestNum и извлекаю коэффициенты

dat_fit1 <- with(dat,
        by(dat[,1], TestNum, fitdist, "lnorm"))
dat_fit2 <-t(sapply(dat_fit1, coef))

Я хочу протестировать другие дистрибутивы, но мне потребуется статистика согласия (gofstat; например, «chi», «chisqpvalue», «cvm», «ad», «ks», «aic», «bic») для каждого подобранная кривая TestNum Я могу получить «aic» и «bic» с кодом ниже, но не остальную статистику.

gof_f<-do.call(rbind, dat_fit1) 
gof_f<-gof_f[,7:8]

Мы ценим любые предложения!