Notice: On March 31, it was **announced** that Statalist is moving from an email list to a **forum**. The old list will shut down on April 23, and its replacement, **statalist.org** is already up and running.

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

From |
Dirk Enzmann <dirk.enzmann@uni-hamburg.de> |

To |
statalist-digest@hsphsun2.harvard.edu |

Subject |
st: Mata finds nonreal although I can't find a no nonreal |

Date |
Tue, 27 Dec 2011 04:14:46 +0100 |

I just encountered the problem that Mata gives me the error message invsym(): 3253 nonreal found where real required <istmt>: - function returned error although I believe that there is no nonreal value in the matrix. Below is the code that produces the error in the line W1 = A1*invsym(A1'*A1)

Any suggestions for a solution of this problem? * -------- (start) ---------- mata P1 = ( .43506, .44920 \ .68702, .06891 \ .77649, -.09335 \ -.60949, .27331 \ .02898, .65232 \ -.16600, .43717 \ .82794, -.04080 \ .79320, -.11331 \ -.04388, .62171 \ .03556, .71095 \ .08361, .75362 \ -.23253, .58472 \ .71299, -.20058 \ -.19725, .55294 \ .51660, .03276 ) RF1 = ( 1.00, .00 \ .00, 1.00 ) S1 = P1*RF1 S1 V1 = NULL eigv1 = NULL eigensystem(RF1,V1,eigv1) A1 = S1*(V1:/sqrt(eigv1))*V1' A1 rssqA1 = diagonal(A1*A1') rssqA1 cssqA1 = diagonal(A1'*A1) cssqA1 W1 = A1*invsym(A1'*A1) W1 end * -------- (end) ------------ Dirk ======================================== Dr. Dirk Enzmann Institute of Criminal Sciences Dept. of Criminology Rothenbaumchaussee 33 D-20148 Hamburg Germany phone: +49-(0)40-42838.7498 (office) +49-(0)40-42838.4591 (Mrs Billon) fax: +49-(0)40-42838.2344 email: dirk.enzmann@uni-hamburg.de http://www2.jura.uni-hamburg.de/instkrim/kriminologie/Mitarbeiter/Enzmann/Enzmann.html ======================================== * * For searches and help try: * http://www.stata.com/help.cgi?search * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

**Follow-Ups**:**st: Mata problem with nonreals [was: Re: Mata finds nonreal although I can't find a no nonreal]***From:*Dirk Enzmann <dirk.enzmann@uni-hamburg.de>

- Prev by Date:
**st: XTIVREG2: Are my instruments proper?** - Next by Date:
**RE: st: XTIVREG2: Are my instruments proper?** - Previous by thread:
**st: XTIVREG2: Are my instruments proper?** - Next by thread:
**st: Mata problem with nonreals [was: Re: Mata finds nonreal although I can't find a no nonreal]** - Index(es):