Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: st: nlcom and delta method
From
Beatrice Crozza <[email protected]>
To
[email protected]
Subject
Re: st: nlcom and delta method
Date
Tue, 3 Aug 2010 19:07:53 +0100
Dear Maarten,
thank you very much for your help.
It works perfectly!
The only thing is that when I calculate one non linear combination and
I save it with the post option and then I calculate another no linear
combination, after a different maximum likelihood, I cannot compare
the two values with the test command.
Do you have any idea of how to solve this problem, or there is no solution?
Because in order to calculate the test by myself I should also need
the covariance.
Thanks a lot.
Bea
2010/8/3 Maarten buis <[email protected]>:
> --- On Tue, 3/8/10, Beatrice Crozza wrote:
>> I did a maximum likelihood, using delta method, to get the
>> values.
>>
>> So, for example, one result of my maximum likelihood using
>> delta method is this one:
>>
>>
>> . nlcom (a: 1/(1+exp(-[ar]_b[_cons])))
>>
>> Now I wan to test a no linear combinations, using the
>> coefficient of a , however the usual command for the
>> coefficient: _b[a] doesn't work.
>
> if you want to access _b[a] then the command that created
> a must be an estimation command, and by default -nlcom-
> is not an estimation command. You can turn -nlcom- into
> an estimation command by specifying the -post- option.
> However now the results of your previous maximum
> likelihood command are overwritten. So this won't work if
> you want to use parameters from your maximum likelihood
> command. In that case it is probably best to do the entire
> opperation in one -nlcom- command. In that case it is
> often convenient to store a string containing the command
> used to create a (n.b. not the value but the command) in
> a local. Say that your second non-linear combination is
> a / [ar]_b[x], then I would do that as follows:
>
> local a "( 1/(1+exp(-[ar]_b[_cons])) )"
> nlcom (b: `a'/[ar]_b[x])
>
> Hope this helps,
> Maarten
>
> --------------------------
> Maarten L. Buis
> Institut fuer Soziologie
> Universitaet Tuebingen
> Wilhelmstrasse 36
> 72074 Tuebingen
> Germany
>
> http://www.maartenbuis.nl
> --------------------------
>
>
>
>
>
> *
> * 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/
>
*
* 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/