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: multilevel mixed effects piecewise model with 1 unknown knot location
From
Jordan Silberman <[email protected]>
To
[email protected]
Subject
Re: st: multilevel mixed effects piecewise model with 1 unknown knot location
Date
Wed, 23 Jan 2013 16:25:47 -0500
Many thanks for your advice. Jordan
On Wed, Jan 23, 2013 at 8:51 AM, Maarten Buis <[email protected]> wrote:
> On Wed, Jan 23, 2013 at 1:12 PM, Jordan Silberman wrote:
>> Does anyone know of a way to estimate in Stata a multilevel model that
>> is piecewise, with 1 knot, at an unknown location?
>
> That will require quite some custom programming. The problem is that
> multilevel models are estimated using some form of maximum likelihood,
> and the likelihood of such a hockey stick model is not smooth; at the
> maximum you won't be able to compute the first and second derivative.
> The tools in Stata for maximizing the likelihood require that the
> first and second derivatives exist, so not only will you have to write
> your own likelihood evaluator, but you also have to rewrite/adapt the
> tools for maximizing that likelihood. I would say that is way too much
> work to be realistic.
>
> What you could do is a grid search, i.e. estimate many different
> models with the knot at different locations and pick the one that fits
> best
>
> -- Maarten
>
> ---------------------------------
> Maarten L. Buis
> WZB
> Reichpietschufer 50
> 10785 Berlin
> Germany
>
> http://www.maartenbuis.nl
> ---------------------------------
> *
> * For searches and help try:
> * http://www.stata.com/help.cgi?search
> * http://www.stata.com/support/faqs/resources/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/faqs/resources/statalist-faq/
* http://www.ats.ucla.edu/stat/stata/