regsubsets 我在做什么错？答案

【问题标题】：What am I doing wrong with regsubsets?regsubsets 我在做什么错？
【发布时间】：2015-09-25 04:45:57
【问题描述】：

我的任务是：

使用 jumps 包中的 regsubsets 函数进行穷举搜索

对于最佳子集回归模型。然后比较为每个选择的调整后的r^2

子集大小。根据这个标准，哪个模型最好？你将不得不

查看summary.regsubsets 的组件。在那里你可以找到相关的。

在每种尺寸下选择的最佳模型的不同优化标准的值。

我有一个我调用的数据集 cigs

q=regsubsets(Sales~Age+HS+Income+Black+Female+Price, data=cigs, method="exhaustive")

所有这些都是正确的变量

summary(q)

Subset selection object
Call: regsubsets.formula(Sales ~ Age + HS + Income + Black + Female + 
    Price, data = cigs, method = "exhaustive")
6 Variables  (and intercept)
       Forced in Forced out
Age        FALSE      FALSE
HS         FALSE      FALSE
Income     FALSE      FALSE
Black      FALSE      FALSE
Female     FALSE      FALSE
Price      FALSE      FALSE
1 subsets of each size up to 6
Selection Algorithm: exhaustive
         Age HS  Income Black Female Price
1  ( 1 ) " " " " "*"    " "   " "    " "  
2  ( 1 ) " " " " "*"    " "   " "    "*"  
3  ( 1 ) "*" " " "*"    " "   " "    "*"  
4  ( 1 ) "*" " " "*"    "*"   " "    "*"  
5  ( 1 ) "*" " " "*"    "*"   "*"    "*"  
6  ( 1 ) "*" "*" "*"    "*"   "*"    "*"

知道为什么这没有给我任何关于r^2的信息吗？

【问题讨论】：

标签： r regression subset

【解决方案1】：

您看到的只是您的regsubsets 对象的打印摘要。如果您输入?regsubsets，您会看到summary.regsubsets 返回的组件比打印的组件多。例如，要访问r^2，只需键入：

summary_of_q <- summary(q)
summary_of_q$rsq

【讨论】：