Несоответствующие массивы |Ошибка расчета Т в квадрате |р

Question

Я пытаюсь вычислить T в квадрате.У меня есть следующие параметры:

> invS #inverse variance covariance matrix          
          x1        x2
x1  0.005536320 -0.001167908
x2 -0.001167908  0.002635186

> n # number of rows
[1] 11

> d_mean
   x1        x2 
-9.363636 13.272727 

Когда я пытаюсь вычислить T в квадрате:

> Tsq <- n* d_mean*invS*t(d_mean)

... Я получаю эту ошибку:

Error in n* d_mean*invS*t(d_mean) : non-conformable arrays

Что я делаю неправильно?

Answer 1

Попробуйте:

Tsq <- n * d_mean %*% invS %*% t(d_mean). 

%*% - матричное произведение, а * - поэлементное произведение.