如何修改BAS模型函数(R)中的向量

如何修改BAS模型函数(R)中的向量,r,vector,lm,R,Vector,Lm,我试图运行一个bas.lm模型,我想指定一个空模型,但我一点也不懂代码 我几乎建立了整个模型,但我希望第一个模型只包含预测因子1、2、3,第二个模型包含所有3+新预测因子。为此,我认为使用了“bestmodel”参数,但我不明白它是如何工作的。文档中说它是:“可选的二进制向量,表示初始化采样的模型。如果空采样从空模型开始”,空采样的源代码是 bestmodel <- c(1, rep(0, p - 1)) bestmodel Mean_Accuracy Predictor1 Predic

我试图运行一个bas.lm模型,我想指定一个空模型,但我一点也不懂代码

我几乎建立了整个模型,但我希望第一个模型只包含预测因子1、2、3,第二个模型包含所有3+新预测因子。为此,我认为使用了“bestmodel”参数,但我不明白它是如何工作的。文档中说它是:“可选的二进制向量,表示初始化采样的模型。如果空采样从空模型开始”,空采样的源代码是

bestmodel <- c(1, rep(0, p - 1))
bestmodel
Mean_Accuracy Predictor1 Predictor2 Predictor3   NewPredictor
0.73             3          7           6             8
0.50             4          5           1             5
0.86             1          3           4             1

Accuracy <- bas.lm(Mean_Accuracy ~ .,
                     data = df,
                     method = "BAS",
                     bestmodel = NULL,
                     prior = "JZS",
                     include.always = Mean_Accuracy ~ Predictor1 + Predictor2 + Predictor3, 
                     #bestmodel = c(?,?),
                     modelprior = uniform())

summary(Accuracy)
           P(B != 0 | Y)   model 1       model 2
Intercept              1    1.0000  1.000000e+00
NewPredictor           1    1.0000  0.000000e+00
Predictor1             1    1.0000  1.000000e+00
Predictor2             1    1.0000  1.000000e+00
Predictor3             1    1.0000  1.000000e+00
BF                    NA    1.0000 4.642017e-139
PostProbs             NA    1.0000  0.000000e+00
R2                    NA    0.3506  2.710000e-01
dim                   NA    7.0000  6.000000e+00
logmarg               NA 1179.9872  8.614630e+02