Вы можете использовать класс "lmList"
, который возвращается функцией lmList
из пакета nlme:
fit1 <- lm(Sepal.Length ~ Sepal.Width, data = iris)
fit2 <- lm(Petal.Length ~ Petal.Width, data = iris)
fits <- list(fit1, fit2)
class(fits) <- "lmList"
library(nlme)
confint(fits)
#[[1]]
# 2.5 % 97.5 %
#(Intercept) 5.579865 7.47258038
#Sepal.Width -0.529820 0.08309785
#
#[[2]]
# 2.5 % 97.5 %
#(Intercept) 0.9393664 1.227750
#Petal.Width 2.1283752 2.331506
См. help("lmList")
для списка функций с методами для lmList
объектов , (Вам необходимо проявлять особую осторожность, если у подгонки нет идентичных формул, как в случае с примером выше.)