Thank you very much Mr Lucchetti. Mr Cottrell also replied and suggested another way, using<br><br>scalar meq = (m1=m2) with values 1 if right, NA if wrong.<br><br>I&#39;ll be checking with both, then. <br><br>Thank you very much!<br>
<br><br><div class="gmail_quote">2009/5/11 Riccardo (Jack) Lucchetti <span dir="ltr">&lt;<a href="mailto:r.lucchetti@univpm.it">r.lucchetti@univpm.it</a>&gt;</span><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div class="im">On Mon, 11 May 2009, Sebastiano Putoto wrote:<br>
<br>
</div><div><div></div><div class="h5"><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Dear Gretl users,<br>
<br>
I know my question may seem (it actually is) very banal, but I can&#39;t figure out a way to check if two matrices are equal.<br>
<br>
The point is verifying idempotence and symmetry for an orthogonal projector. The mathematical proof isn&#39;t required here: the point is showing that the actual calculated<br>
matrices are orthogonal projectors and therefore have the specified properties; I got to show the numbers are the same. And apart from arbitrary extraction of submatrices to<br>
prove they get out same, I can&#39;t think of anything &#39;scientific&#39;.<br>
<br>
Hope you will be able to help me with the matter.<br>
<br>
Regards,<br>
<br>
</blockquote>
<br></div></div>
How about<br>
<br>
   check = sumc(sumr(abs(mat1-mat2))) &lt; 1.0e-07<br>
<br>
or something similar?<br>
<br>
<br>
Riccardo (Jack) Lucchetti<br>
Dipartimento di Economia<br>
Universitą Politecnica delle Marche<br>
<br>
<a href="mailto:r.lucchetti@univpm.it" target="_blank">r.lucchetti@univpm.it</a><br>
<a href="http://www.econ.univpm.it/lucchetti" target="_blank">http://www.econ.univpm.it/lucchetti</a><br>_______________________________________________<br>
Gretl-users mailing list<br>
<a href="mailto:Gretl-users@lists.wfu.edu">Gretl-users@lists.wfu.edu</a><br>
<a href="http://lists.wfu.edu/mailman/listinfo/gretl-users" target="_blank">http://lists.wfu.edu/mailman/listinfo/gretl-users</a><br></blockquote></div><br><br clear="all"><br>-- <br>Cordiali saluti/Salutations/Best Regards<br>
<br>Sebastiano Putoto<br>