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