我正在用 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)
任何建议将不胜感激。