Bookmark and Share

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, is already up and running.

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

st: numeric accuracy

From   Stas Kolenikov <>
Subject   st: numeric accuracy
Date   Tue, 21 Sep 2010 12:15:54 -0500

I run a certain, rather extensive, microsimulation project on several
flavors of Stata (MP8 down to SE) and try to benchmark my results. On
some of the numeric outcomes I compare (percentiles of a
distribution), I get differences between different flavors as large as
300*c(epsdouble). The categorical outcomes always match, and some
other numeric outcomes (mean of a distribution), the difference is
~2*c(epsdouble). But the factor of 300 bothers me. I wonder what the
typical constants are in front of the c(epsdouble) that are used in
certification scripts. When I compared results of my -confa- with
-gllamm-, I was happy to find the difference to be of the order 1e-4.
But these are different estimation methods; here I am running exactly
the same code on exactly the same data, and I hoped to see a close
reproducibility from one run to another.

Stas Kolenikov, also found at
Small print: I use this email account for mailing lists only.
*   For searches and help try:

© Copyright 1996–2016 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   Site index