Why when using priors, greta doesnt recognize matrix multiplication? is there a workaround?
R
(X <- matrix(1:10, ncol=1))
[,1]
[1,] 1
[2,] 2
[3,] 3
[4,] 4
[5,] 5
[6,] 6
[7,] 7
[8,] 8
[9,] 9
[10,] 10
(B <- matrix(1:3, nrow=1))
[,1] [,2] [,3]
[1,] 1 2 3
X%*%B
[,1] [,2] [,3]
[1,] 1 2 3
[2,] 2 4 6
[3,] 3 6 9
[4,] 4 8 12
[5,] 5 10 15
[6,] 6 12 18
[7,] 7 14 21
[8,] 8 16 24
[9,] 9 18 27
[10,] 10 20 30greta
(X <- as_data(X))
greta array (data)
[,1]
[1,] 1
[2,] 2
[3,] 3
[4,] 4
[5,] 5
[6,] 6
[7,] 7
[8,] 8
[9,] 9
[10,] 10
(B <- as_data(B))
greta array (data)
[,1] [,2] [,3]
[1,] 1 2 3
X%*%B
[,1] [,2] [,3]
[1,] 1 2 3
[2,] 2 4 6
[3,] 3 6 9
[4,] 4 8 12
[5,] 5 10 15
[6,] 6 12 18
[7,] 7 14 21
[8,] 8 16 24
[9,] 9 18 27
[10,] 10 20 30B is prior
(B <- dirichlet(alpha=ones(1,3)))
greta array (variable following a dirichlet distribution)
[,1] [,2] [,3]
[1,] ? ? ?
X%*%B
[,1] [,2] [,3]
[1,] NA NA NA
[2,] NA NA NA
[3,] NA NA NA
[4,] NA NA NA
[5,] NA NA NA
[6,] NA NA NA
[7,] NA NA NA
[8,] NA NA NA
[9,] NA NA NA
[10,] NA NA NA