将似然方程转换为 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 个和号有关。我找不到一个示例来说明如何处理 2 个和号并遵循它们的结构。

我试图逐个调试它以检查我哪里出错了，每次我到达这个部分时，我得到错误：

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

任何建议将不胜感激。

Answer 1

尝试用 %*% 代替 * 进行矩阵乘法。

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

将似然方程转换为 R 代码

Converting likelihood equation into R code

equation

r

bayesian

log-likelihood