Take the 2-minute tour ×
Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. It's 100% free, no registration required.

I want to implement the following formula (taken from Kaiser, 1970) in R where $R$ is square matrix of correlations:

$$S = (\textrm{diag } R^{-1})^{-1/2}$$

I understand the diagonal and inverse operations, but I am unclear on the meaning of raising a square matrix to a negative half power. Thus, my questions

  1. What does it mean to raise a square matrix to a negative half power?
  2. What general ideas of linear algebra does this assume?
  3. (if this is in scope) How would this be implemented in R?

References

  • Kaiser, H. F. (1970). A second generation little jiffy. Psychometrika, 35(4), 401-415.
share|improve this question

2 Answers 2

up vote 4 down vote accepted

If you look at Powers of diagonal matrices, it would be the reciprocal of the square root of each term in the diagonal.

Lets do an example for a diagonal matrix:

$$A=\begin{bmatrix}2 & 0 \\ 0 & 3\end{bmatrix}$$

$$A^{-1/2} = \begin{bmatrix} \frac{1}{\sqrt{2}} & 0 \\ 0 & \frac{1}{\sqrt{3}}\end{bmatrix} = \frac{1}{6}\begin{bmatrix}3 \sqrt{2} & 0 \\ 0 & 2 \sqrt{3}\end{bmatrix}$$

Clear?

share|improve this answer
    
Many thanks. I should have realised that. –  Jeromy Anglim Mar 25 '13 at 4:45
    
@JeromyAnglim: no problem - I am blind to most things until I work them! Regards –  Amzoti Mar 25 '13 at 4:48
1  
Clear to me! ;-) –  amWhy Apr 13 '13 at 0:27
    
Same here...and it can be embarrassing when I miss something in answers here on math.se! So it it does keep oneself humble. –  amWhy Apr 13 '13 at 2:26

The general idea (I think so) is Matrix function and Taylor expansion. Let $R$ is square matrix, we may write it in the form $R = I+X$ ($I$ - identity matrix), then using Taylor expansion: $$ R^{-1/2} = (I+X)^{-1/2} = I -\frac{1}{2}X + \frac{3}{8}X^2 + \cdots $$ I would be nice if $X$ is small.

share|improve this answer
1  
I see, that I've answered bit another question, not actually yours. The question I've answered is: what is $A^{-1/2}$ when $A$ is a square matrix. ) –  nikita2 Mar 25 '13 at 5:06
1  
It is worth noting that the inverse square root of a square matrix is not unique, and that the one obtained using this particular taylor expansion is just one possible such matrix. –  Glen O Mar 25 '13 at 5:31
    
Yes, even $\sqrt{1}$ is not unique. –  nikita2 Mar 25 '13 at 6:27
    
Thanks for the answer. The power of a diagonal diagonal matrix was my specific problem, but I was interested in reading about the more general problem involving any square matrix. –  Jeromy Anglim Mar 25 '13 at 6:57

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.