r - 複数のパラメータの下で尤度計算をベクトル化する方法は?

Question

私はベルヌーイ混合を実装しようとしていますが、ループせずに計算を正しくベクトル化する方法を考えていました。

さまざまなバージョンの apply を試しましたが、目的の出力が得られません (dim = c(5,4,2)。コンポーネントのパラメーターはマトリックスではなくリストにする必要がありますか?

set.seed(123)

#Data
X <- matrix(sample(c(0,1), 20, replace = TRUE, prob = c(.6, .4)),
               nrow = 5, ncol = 4)

#Params
parameters <-  matrix(runif(8), nrow = 2, ncol = 4)

#Would like to vectorize this
dbinom(X, 1, parameters[1,], log = TRUE)
dbinom(X, 1, parameters[2,], log = TRUE)