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

Community Tip - Your Friends List is a way to easily have access to the community members that you interact with the most! 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