限制R的BMA包可以处理的解释变量的数量?

时间:2013-12-29 08:05:24

标签: r linear-regression bayesian

使用R的BMA(贝叶斯模型平均)包,我想运行以下代码:

result = bic.glm(x,y,prior.param = c(1,1,1,1,0.5,1,0.5,0.5,0.5,1,1,1,1,1,0.5,1,
1,1,1,1,1,1,1,1,1,1,1,1,0.5,1), glm.family="gaussian",factor.type=TRUE)

当我的x矩阵有29个变量(并且我的先前参数向量也长度为29)时,bic.glm函数有效。但是,当我向x添加第30个变量并向我的prior.param向量添加第30个值时,我收到以下错误:

Error in dropcols(leaps.x, y, glm.family, wt, maxCol) : dropped == 0

这是否意味着此函数可以处理的解释变量的最大数量?或者错误是否与完全不同的东西(可能缺少值/ NA)有关?指导非常感谢!

0 个答案:

没有答案