从 regsubsets 中提取详细信息

【问题标题】：Extract details info from regsubsets从 regsubsets 中提取详细信息
【发布时间】：2018-03-14 08:12:41
【问题描述】：

给定一个数据集，如果我们在 R 中使用 regsubsets 函数，我们只能使用像 Adjusted R Squared 这样的比例来绘制结果：

    data(iris)
    a = regsubsets(Petal.Width~.,data=mmat,nbest=5)
    plot(a,scale="adjr2")

有没有办法获得像跳跃功能这样的详细信息？变量“which”和“adjr2”等细节？

    data(iris)
    a = leaps(x=iris[,-5], y=iris[,5], method=c("adjr2"), names=names(iris[-5]))
    print(a$which)
        print(a$adjr2)

【问题讨论】：

【解决方案1】：

您可以通过将summary() 函数应用于regsubsets() 来获取此信息。该方法允许计算$which、$rsq、$adjr2、$cp、$bic等。

library(leaps)
data(iris)
a <- summary(regsubsets(Petal.Width~.,data=iris,nbest=5))
summary(a)
a$which
a$rsq
a$adjr2

【讨论】：

【解决方案2】：

我不确定您到底在问什么，但查看regsubsets()https://cran.r-project.org/web/packages/leaps/leaps.pdf 上的文档，我们可以看到该函数还支持which 以及调整后 r 平方以外的指标（例如 Cp , r 平方, BIC)。

【讨论】：