r中矩阵之间的计算

Question

我有一个与处理两个矩阵有关的问题。

假设第一个矩阵（P（i））是：

x <- rmultinom(2, size = 10, prob=c(0.5,0.5))

> x
      [,1] [,2]
[1,]    6    4
[2,]    4    6

..和第二个矩阵（Q（i））是：

y <- rmultinom(2, size = 10, prob=c(0.5,0.5))
> y
      [,1] [,2]
[1,]    6    7
[2,]    4    3

假设P（i）代表第一矩阵中的每个单元，Q（i）代表第二矩阵中的每个单元。我试图总结这个功能，（我找不到＆＃34; Sum＆＃34的符号）：

z = (Sum)(P(i)*log(P(i)/Q(i)))

非常感谢任何帮助。 谢谢！