topic Re: PHREG multilevel response in Statistical Procedures

PHREG multilevel response

lillymaginta1 — Sun, 27 Dec 2020 14:02:02 GMT

Hi Everyone,

How can I modify PHREG statement to run for a multilevel response variable (stm= 0, 1, 2). The following statement yield one HR, I need to get two HR for stm(1 vs. 0 and 2 vs. 0)

proc phreg data=s ;

   class drug (ref='1') stm (ref='0') / param=ref;

   model period1*stm(0)= drug  year/ risklimits;

run;

any help is appreciated

Re: PHREG multilevel response

SASJedi — Mon, 04 Jan 2021 15:58:38 GMT

Moved this post to Statistical Procedures, where it should get better exposure.

Re: PHREG multilevel response

SteveDenham — Tue, 05 Jan 2021 18:54:43 GMT

Try adding a HAZARDRATIO statement, such as:

hazardratio stm/diff=ref;

SteveDenham

Re: PHREG multilevel response

lillymaginta1 — Wed, 06 Jan 2021 14:51:48 GMT

Thank you for the response! I tried the statement but I received the following error, the variable specified in the hazardratio statement must be in the model.

Re: PHREG multilevel response

StatDave — Wed, 06 Jan 2021 16:12:07 GMT

If the response that you want to model is multinomial (either ordinal or nominal), then you should be using PROC LOGISTIC, not PROC PHREG. You can find examples of both in the Examples section of the LOGISTIC documentation. From your description, it sounds like you want to treat the response levels as nominal. The comparisons you mention will be obtained if you use the LINK=GLOGIT option in the MODEL statement and the DESCENDING response variable option:

model stm(descending) = drug year / link=glogit;