我在 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 函数,所以还有另一个命令可以使用吗?帮助表示赞赏!