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

Community Tip - You can subscribe to a forum, label or individual post and receive email notifications when someone posts a new topic or reply. Learn more! X

Errors working with a largematrix

ptc-5412730
1-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

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?

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

Top Tags