我在Splus中有一个代码,但必须将其转换为R,这不是一件大事。但是我对这两个软件都很陌生。这是我正在努力的代码:
a <- leaps(x, y, wt = wt, method = "adjr2", nbest=nbest,
names = dimnames(x)[[2]], int=F)
x是一个独立变量矩阵,其中第一个长度(keep1)列对应于始终保存在BMA中的变量(贝叶斯模型求平均值 - 这不重要。本质上,x是矩阵),y也是一个矩阵。
当我选择我的数据并运行该功能时,我收到了以下错误:
leaps.exhaustive(a)中的错误:穷举搜索将是S L O W, 必须指定really.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,方法=&#34; adjr2&#34;,nbest = nbest,names = dimnames(x)[[2]] ,: unused arguments(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")
所以当我使用really.big时,R告诉我这是一个未使用过的论点,但是告诉我在我不使用时使用它。这里发生了什么?任何帮助表示赞赏!
答案 0 :(得分:0)
在此之前添加 T <- 60
a <- 跳跃(x, y, wt = wt, method = "adjr2", nbest=nbest, names = dimnames(x)[[2]], int=F,really.big=T)