我在Splus中有一个代码,但必须将其转换为R,这不是一件大事。但是我对这两个软件都很陌生。这是我正在努力的代码:
a <- regsubsets(x, y, wt = wt, method = "adjr2", nbest=nbest,
names = dimnames(x)[[2]], int=F)
x是一个独立变量矩阵,其中第一个长度(keep1)列对应于始终保存在BMA中的变量(贝叶斯模型平均值 - 这不重要。本质上,x是矩阵),y是a矩阵也是。
最初,这是使用“跳跃”功能,但矩阵太大,所以R让我使用regsubsets。现在,它说“adjr2”是一个“模糊或无法识别的方法名称”。我很确定这是飞跃功能中公认的方法吗?是否有另一个命令要使用,因为我使用的是regsubset函数,而不是leaps函数?帮助赞赏!
答案 0 :(得分:0)
regsubsets()
与leaps()
的参数不同。对于regsubsets()
,“方法”指的是模型选择的类型。来自https://cran.r-project.org/web/packages/leaps/leaps.pdf文档的第5页:
方法:使用穷举搜索,前进选择,后退选择或 顺序替换搜索。
所以adjr2
不是regsubsets()
的选项。如果要从regsubsets()
的模型搜索输出中获取adjr2,可以这样调用:
> summary(regsubsets_object)$adjr2