[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]
SV: RE: Re: st: Linear Trend Tests of ORs
Young Hee Rho wrote:
I have encountered many "trend tests" of linearity concerning odds ratios (OR) of a categorical variable.
For example, I am modeling a logistic model Y=b1x1 + b2x2 + b3x3 +b4. x2 is a 5-level categorical variable, for example the level of drinking (while Y is the presence/absence of hyperuricemia). When the results are displayed, the ORs of the 5 levels are shown and the linear trend is shown as a single p value. The individual ORs may not have significance, however the overall trend does.
To do a test for linear trend, you may use the log likelihood ration test!
First you run the logistic regression with the categorical variable expanded by 'xi' (getting 4 estimates relative to the reference category of b2) and store the log likelihood in 'A':
xi:logistic y b1 i.b2 b3 b4
estimate store A
then you repeat the regression without the 'xi' expansion of the categorical variable (Now you only one estimate of b2, which is the linear effect of b2).
xi:logistic y b1 b2 b3 b4
(Note: the 'i.' in front of b2 is removed).
You then simply need to se if the reduced model (model 2) is as good as your previous model (Model 1). You do that using the likelihood ration test:
To conclude that you have a linear trend the p-value of the lrtest needs to be insignificant (Model 2 is not significantly worse than Model 1) AND the estimate for b2 (the linear effect per category in model 2) must be significant!
Kim Lyngby Mikkelsen
Fra: firstname.lastname@example.org [mailto:email@example.com] På vegne af Rho YH
Sendt: 23. maj 2006 05:23
Emne: RE: RE: Re: st: Linear Trend Tests of ORs
I have just found out that the tabodds command may meet what I wanted - linear trend of ORs,
however making multivariate adjustments is not easy (I tried and it gave no results after adjusting
for >2 or 3 variables.) Is there any "immediate command" by just inputing the OR (and CI, if needed)
and the independent variable category and produces a p value?
> ---- Original Message ----
> From : Rho YH [firstname.lastname@example.org]
> To : email@example.com
> Date : 2006년 5월 23일(화) 09:43:23
> Subject : RE: Re: st: Linear Trend Tests of ORs
>Hmm.. It looks like the aformentioned Cochrane-Armitage Test, however I'll check it out.
>> ---- Original Message ----
>> From : Suzy [firstname.lastname@example.org]
>> To : email@example.com
>> Date : 2006년 5월 22일(월) 21:24:22
>> Subject : Re: st: Linear Trend Tests of ORs
>>Perhaps Szklo and Nieto's book can help: Epidemiology. Beyond the
>>Basics, discusses test for trend (dose reponse) in Appendix B (pp 459-462).
>>Formula is from Mantel:
>>Mantel N. Chi square tests with one degree of freedom: etensions of the
>>Manetel-Haenszel procedure. J Am Stat Assoc. 1963;58: 690-700.
>>Hope this helps.
>>Young Hee Rho wrote:
>>>I have encountered many "trend tests" of linearity concerning odds ratios (OR) of a
>>>For example, I am modeling a logistic model Y=b1x1 + b2x2 + b3x3 +b4. x2 is a 5-level
>>>categorical variable, for example the level of drinking (while Y is the presence/absence of
>>>hyperuricemia). When the results are displayed, the ORs of the 5 levels are shown and
>>>the linear trend is shown as a single p value. The individual ORs may not have significance,
>>>however the overall trend does. It is said that it was tested through regressing the median of
>>>the levels on the ORs. Otherwise in other cases, there are many trend tests of linearity
>>>expresed in many papers, however, the actual method is not explained in detail. (It does not
>>>apear to come from polynomial contrasts of ANOVA nor from categorical trend tests
>>>(Cochrane-Armitage) since the arformentioned test is from values coming from
>>>one categorical variable having several estimates. How is this done and how much methods
>>>exsist on this topic? Are there any useful references?
>>>** For those who got twice this article, I sent this article again since it did not seem to register on
>>>Statalist. Many apologies if there was a duplicate delivery.
>>* For searches and help try:
* For searches and help try: