Я работаю над кодированием модели GP в R и вычислением ее функции правдоподобия.
У меня проблема с преобразованием приведенного ниже уравнения в код 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)
Я думаю, что это связано с двумя знаками суммы. Я не смог найти пример, показывающий, как работать с двумя знаками сумм и следовать их структуре.
Я пытался отладить его по частям, чтобы проверить, где я ошибся, и каждый раз, когда я добираюсь до этой части, я получаю сообщение об ошибке:
(t(zlt -olt)) * (zlt- olt)
Любой совет будет высоко ценится.