我想使用dlm
用一个以上的解释变量设置一个DLM。
包。
我已经使用dlmGibbsDIG
函数通过MCMC方法从观察方程和系统方程的方差进行了估算。
输出给出了我的dV
值和对角矩阵dW
。现在,我想使用dlm class
函数用我的信息构建一个dlmModReg
:
mod <- dlmModReg(X = test_data, addInt = T, dV = dV, dW = dW)
test_data
包含我的设计矩阵,其中包含20个变量。 dV
包括我的观察方程式(单值)的方差,dW
包括我的数据集中每个自变量的方差(21x21,带截距)。
运行代码行会导致以下错误消息:
Error in dlmModReg(X = test_data, addInt = T, dV = dV, dW = dW) :
Inconsistent dimensions of arguments
如果我在不使用dW
矩阵的情况下运行代码,则dlm class
是使用21x21 W
矩阵创建的,每个条目均为零(因此)。因此,我不理解我的dW
对角矩阵引起的错误消息。
但是如果没有这些信息,我将无法执行包括过滤和平滑在内的以下步骤来估计自变量的回归系数。
感谢您的帮助和解答。