topic Re: Can I get stratum-specific estimates by using strata statement in proc phreg? in Statistical Procedures

Can I get stratum-specific estimates by using strata statement in proc phreg?

Christina96 — Tue, 31 Jan 2023 15:57:55 GMT

Hi! I am doing stratified cox proportional hazard regression analyses. I used by statement with proc sort stated ahead.

My colleague told me to use strata statement. And I tried but I did not get stratum-specific results. Instead, I only got one HR (assuming there were 3 strata). I wondered how to do this with strata statement. Or, I should stick with the by statement with proc sort.

Thank you in advance!

Here's the code:

proc phreg data=mydata namelen=32;
class A(ref="1") B(ref="0") C(ref="aa");
model personyrs*Y(0)=X
					 A B C/rl;      
strata D;
hazardratio X;
run;

Re: Can I get stratum-specific estimates by using strata statement in proc phreg?

OsoGris — Wed, 01 Feb 2023 14:01:05 GMT

You do not get stratum specific estimates using the STRATA statement in PROC PHREG. PHREG multiplies the partial likelihoods from each strata and then optimizes (determines the coefficients) from the resulting product. It sounds like maybe you want to use a BY statement instead to get estimates specific to a particular group or groups.

Re: Can I get stratum-specific estimates by using strata statement in proc phreg?

Christina96 — Wed, 01 Feb 2023 15:36:01 GMT

Thank you very much! Others also suggested that WHERE statement would work. I wondered if there's any difference. Based on my understanding, I can get three stratum-specific estimates by running PHREG once. But for WHERE, I need to run it 3 times if I have 3 strata. Please correct me if I am wrong. I did not find any detail explanation about WHERE statement under PHREG in SAS help center.

Thank you!

Re: Can I get stratum-specific estimates by using strata statement in proc phreg?

OsoGris — Wed, 01 Feb 2023 16:39:05 GMT

I don't think you want a WHERE statement which would mean you would have to run PHREG 3 times to get 3 sets of estimates. Whereas a BY statement can handle this with one run of PHREG.

Re: Can I get stratum-specific estimates by using strata statement in proc phreg?

Christina96 — Wed, 01 Feb 2023 18:11:55 GMT

Thank you very much! I also think BY would be the best fit but I was not sure about it. And, thanks for the explanation on the STRATA statement in PHREG.