cancel
Showing results for 
Search instead for 
Did you mean: 
cancel
Showing results for 
Search instead for 
Did you mean: 

Errors working with a largematrix

Newbie

Errors working with a largematrix

Z is a 51 row X 3 column matrix. I need to obtain the 51 X 51 matrix H = Z(Z'Z)^(-1)Z', where ' denotes transpose and ^(-1) means inverse. H is idempotent, so as a check, either H'H or HH' should yield an identity matrix. Straightforward evaluation of H evidently yields an incorrect result, as H'H is not an identity matrix.. Can MathCad handle matrices this large? Or are there special procedures for dealing with such matrices?

2 REPLIES 2

Re: Errors working with a largematrix

Matrices of that size should be no problem for Mathcad or for Prime and there are no special treatments.

Why don't you tell us which version of Mathcad or Prime you are using and why not attach your worksheet showing what you think is wrong?

Re: Errors working with a largematrix

Sorry, I just read your post more carefully.

H is idempotent, so as a check, either H'H or HH' should yield an identity matrix.

Why? Doesn't idempotent simply mean that H^2 = H?

This is true for H:

idempotence.png

Announcements