0

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

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

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

当我选择我的数据并运行这是其中一部分的函数时,我收到以下错误:

jumps.exhaustive(a) 中的错误:穷举搜索会很慢,必须指定 real.big=T

因此,我将初始命令更改为:

 a <- leaps(x, y, wt = wt, method = "adjr2", nbest=nbest, 
 names = dimnames(x)[[2]], int=F, really.big=T)

当我再次运行该函数时,出现此错误:

跳跃错误(x,y,wt = wt,method = "adjr2",nbest = nbest,names = dimnames(x)[[2]],:未使用的参数(really.big = T)

该函数以这种方式调用:

out.dkt<-dkt.bma(bigx, bigy, bigz.sys, Tper=3,  hendog, name.x, name.z, 
emethod="LS", info.crit="BIC", h, keep1,  leaps.genR2="F",
nbest=1000, hier=1, q=0.50, OW=0,OR=1000,
filo="U:\\Programanddatafiles\\dkt-jae-ls-bic-hier1-sys-nokeep-nbest-check.txt")

所以当我使用真的.big 时,R 告诉我这是一个未使用的参数,但告诉我在我不使用时使用它。这里发生了什么?任何帮助表示赞赏!

4

1 回答 1

1

在此之前添加 T <- 60

a <- 跳跃(x,y,wt = wt,method = "adjr2",nbest=nbest,names = dimnames(x)[[2]],int=F,really.big=T)

于 2021-03-16T21:56:38.460 回答