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

Question

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

> 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). 

%*% это матричное произведение и * поэлементный продукт.