Greta中的元素乘法

Question

我正在用Greta编写一个模型，遇到了一个小小的绊脚石:使用Greta数组的元素乘法。

可重现性误差：

library(greta)
library(MCMCpack)

以下工作(此处不涉及Greta数组)：

pi <- rdirichlet(1, alpha = rep(1, 5))
tpi <- t(array(pi, dim = c(5, 10)))
tpi*2

以下内容不起作用(此处涉及的Greta数组)：

pi <- dirichlet(alpha = rep(1, 5), dim = 1)
tpi <- t(array(pi, dim = c(5, 10)))
tpi*2
Error in tpi * 2 : non-numeric argument to binary operator

以下内容也有效(尽管pi不是数组，而tpi是数组)：

2*pi

但是我需要在我的模型中执行数组扩展操作(步骤2)。

谢谢您的建议。

Answer 1

这是一个在Greta中重塑数组的公开发行。这个问题应该在下一个版本中解决，但同时，Greta @goldingn (参见这里的评论)提出了一个备选方案：

而不是：

tpi <- t(array(pi, dim = c(5, 10)))

使用：

tpi <- do.call(rbind, replicate(10, pi))

下面的元素乘法将不再产生错误。