使用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)有关?指导非常感谢!