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]
st: -nlcom- issues
Stas Kolenikov <firstname.lastname@example.org>
st: -nlcom- issues
Mon, 17 May 2010 12:33:45 -0500
I have a couple of issues with -nlcom- that seem odd. To be specific
(but I don't think that matters), I am estimating percentiles of an
income distribution in two years of data (by a custom eclass command),
and wrap them into a growth model (by -nlcom- -- I need a ratio of say
25th percentiles in two waves of data, pro-rated to be on the per
1. Time issue: each next -nlcom- takes longer to run than the previous
one, although at the face of it they should run at the same speed. I
think the issue came up a few months ago on the list, and was traced
down to how -testnl- (on which -nlcom- relies in computation of the
standard errors) handles memory. In the original post, it even crashed
Stata in some simulations when it used all available handles to
temporary objects (or something of that nature). I was under
impression that this oddity was fixed, but apparently it was not.
2. Numeric issue: the standard errors reported by -nlcom- tend to go
down with the number of parameters. That is, I estimate say
percentiles in numlist 10(10)90, and percentiles in numlist 2(2)98,
and I form my growth in the percentile estimates from these two sets
of estimates. While the estimates and their standard errors are
identical between the two lists when they intersect (i.e., for the
percentiles in the first list), the -nlcom- results match only in
point estimates. The standard errors for these identical point
estimates are smaller in the second case. I would say, disturbingly
smaller, by about 20-25%.
I'd be happy to share the code (and probably the data, too) with
tech-support, if needed.
Stas Kolenikov, also found at http://stas.kolenikov.name
Small print: I use this email account for mailing lists only.
* For searches and help try: