Ariel
unread,Apr 21, 2012, 2:06:45 PM4/21/12Sign in to reply to author
Sign in to forward
You do not have permission to delete messages in this group
Sign in to report message
Either email addresses are anonymous for this group or you need the view member email addresses permission to view the original message
to meds...@googlegroups.com
Dear All,
I am trying to build a model to explain an ordinal response variable y with 4 levels: y0, y1, y2 and y3. The independent variable in this model is v.
v is a categorical variable with three categories va, vb and vc. The assumption of proportional odds fails for v.
Therefore I am using the partial proportional odds (PPO) model. The model gives me 3 ORs for y1, y2 and y3
that correspond to the odds for events {y>=y1}, {y>=y2}, and {y>=y3}
I would like to perform a test for trend for each of these 3 events such that each test gives a p-value.
For example, the test for trend for y1 would explain if there is a trend between the change in level of v and in the odds for {y>=y1}.
Is there a way to use the partial PO model directly to compute a trend test for the 3 events?
Here’s a solution I am thinking of:
I first convert the categorical variable, v into a continuous variable with values 1,2, and 3 and build a PPO model using this continuous variable - let’s call this model PPOC.
For each outcome level i, I divide the samples in the data into 2 sets: {y>=yi} and {y<yi}, and compute the likelihood of this split using PPO models with
a) only an intercept term, and b) both intercept and exposure (i.e. PPOC).
The difference in likelihoods of a) and b) gives me a test for this level.
Is this the right approach? Will this test be equivalent to a test for trend?
Are there any other approaches that I can adopt to determine if a trend exists between v and y for each of the 3 levels of y?
Any pointers would be very helpful!