我正在尝试将矩阵数据输入到brm()函数中以运行信号回归。brm来自 brms 包，它提供了一个使用 Stan 拟合贝叶斯模型的接口。信号回归是指在更大的模型中使用另一个协变量对一个协变量进行建模，并使用by如下参数：model <- brm(response ~ s(matrix1, by = matrix2) + ..., data = Data)。问题是，我无法使用“data”参数输入矩阵，因为它只允许data.frame输入一个对象。
这是我的代码以及我试图绕过该约束而获得的错误......
首先，我的可重现代码导致模型构建：
library(brms)
#100 rows, 4 columns. Each cell contains a number between 1 and 10
Data <- data.frame(runif(100,1,10),runif(100,1,10),runif(100,1,10),runif(100,1,10))
#Assign names to the columns
names(Data) <- c("d0_10","d0_100","d0_1000","d0_10000")
Data$Density <- as.matrix(Data)%*%c(-1,10,5,1)
#the coefficients we are modelling
d <- c(-1,10,5,1)
#Made a matrix with 4 columns with values 10, 100, 1000, 10000 which are evaluation …