在 R 中合并两个 "lm" 个对象

Question

我有两个来自线性拟合的 class "lm" 对象；让我们将它们记为 fit1 和 fit2

我想创建一个唯一变量，例如 "fits" of class object "lm" 包含两个对象作为元素。这是为了在拟合的单个分量上应用 "confint" 或 "summary" 函数。

有什么建议吗？

Answer 1

您可以使用 "lmList" class，它由 nlme 包中的 lmList 函数返回：

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

有关 lmList 对象的函数列表，请参阅 help("lmList")。（如果拟合没有与上例中相同的公式，则需要特别小心。）

在 R 中合并两个 "lm" 个对象

Merging two "lm" objects in one in R

merge

r

lm