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: Re: gllamm and u-shaped cumulative logits over time

From   "Joseph Coveney" <>
To   <>
Subject   st: Re: gllamm and u-shaped cumulative logits over time
Date   Wed, 20 Mar 2013 11:14:24 +0900

Thomas Herold wrote:

I am working on a longitudinal dataset dealing with depression, measured 
as ordinal variable (1: none, 2: mild, 3:moderate, 4: severe). There are 
5 measurement points (t0-t4) and two different groups (T0=conventional 
treatment GP, T1=new treatment). I want to find out whether the 
treatment has an influence on the development of the depression score. I 
would normally use the -gllamm- command with the model (without a random 
slope) looking something like this:

gllamm depr week treatment interact, i(id) link(ologit) adapt eform

, where week is the time tx-t0 and interact=week*treatment.

But here is my problem: The development over time (and over time by 
group) is not linear but clearly (!) u-shaped. There is a short-term 
effect in both groups, i.e. the depression scores decline considerably 
from t0 to t1. However, in the long run (t1-t4), the scores in both 
groups increase again and often even exceed the initial value at t0. For 
obvious reasons, the above model fits the data very poorly.

Any help would be much appreciated.


How about indicator variables?

xi: i.treatment*i.week
gllamm depr _I*, i(id) family(binomial) link(ologit) adapt eform

As an alternative, you could try polynomials in time (along with the interaction
terms) or fractional polynomials (-help fracgen-, or even -fracpoly-)

Or splines--there are the official Stata -mkspline- and the user-written
-bspline- suite of commands. 

Joseph Coveney

*   For searches and help try:

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