1

我在 Splus 中有一个代码,但必须将其转换为 R,这不是什么大事。但是我对这两种软件都很陌生。这是我正在努力的代码:

a <- regsubsets(x, y, wt = wt, method = "adjr2", nbest=nbest, 
     names = dimnames(x)[[2]], int=F)

x 是一个自变量矩阵,其中第一个 length(keep1) 列对应于始终保存在 BMA 中的变量(贝叶斯模型平均——这并不重要。本质上,x 是一个矩阵),并且 y 也是一个矩阵。

最初,这是使用“跳跃”函数,但矩阵太大,所以 R 让我使用 regsubsets。现在,它说“adjr2”是一个“模棱两可或无法识别的方法名称”。我很确定这是跳跃功能中公认的方法?因为我使用的是 regsubset 函数,而不是 jumps 函数,所以还有另一个命令可以使用吗?帮助表示赞赏!

4

1 回答 1

0

regsubsets()没有与 相同的参数leaps()。对于regsubsets()“方法”是指模型选择的类型。从https://cran.r-project.org/web/packages/leaps/leaps.pdf文档的第 5 页:

方法:使用穷举搜索、前向选择、后向选择或顺序替换进行搜索。

所以adjr2不是regsubsets(). 如果要从 的模型搜索输出中获取 adjr2 regsubsets(),可以这样调用:

> summary(regsubsets_object)$adjr2
于 2015-09-24T20:36:24.057 回答