将似然方程转换为R码

Question

我正在用R编码GP模型并计算其似然函数。

我在将下面的公式转换为R代码时遇到问题：

我使用了下面的代码，但不断收到错误消息：

Error in (t(zlt - olt)) * (zlt - olt) : non-conformable arrays

我的代码是：

ConditionalLikelihood <- (-(N/2)) * (log(sig2e)) -(1/(2*(sig2e)))*(t(zlt -olt)) * (zlt- olt) -(sum(T*r)/2) *(log (sig2eta * SEta)) (-1/(2*(sig2e))) * (t(zlt- xb)) * 1/SEta *(zlt- xb)

我认为这与两个和号有关。我找不到一个示例，说明如何处理2个和号并遵循其结构。

我试图逐部分调试它，以检查哪里出错了，并且每次进入此部分时，都会出现错误：

(t(zlt -olt)) * (zlt- olt)

任何建议将不胜感激。

Answer 1

尝试％*％进行矩阵乘法，而不是*。

(t(zlt -olt)) %*% (zlt- olt)