BookmarkSubscribeRSS Feed
🔒 This topic is solved and locked. Need further help from the community? Please sign in and ask a new question.
ammarhm
Lapis Lazuli | Level 10

Hi

I know this forum is mainly for help with the SAS code itself and I have received some fantastic answers here, I would really appreciate if anyone could please help me with this.

I am doing a regression analysis, looking at the effect of 2 drugs (A or B) in a population where we stratified individuals by ethnic group:

Estimates for European Population

proc logistic data=Studydatabase;

class sex drug study;

model Efficacy (event='Good')=   drug age sex weight egfr;

where population='European';

run;

EffectOddsRatioEstLowerCLUpperCLProbChiSq
Age1.0050.9731.0370.7729
Drug A vs B1.8621.2292.8210.0033
Sex Female vs Male1.4690.9252.3340.1033
Weight1.0151.0011.0280.0346
eGFR1.0060.9951.0160.2967

________________________________________

Estimates for Asian Population

proc logistic data=Studydatabase;

class sex drug study;

model Efficacy (event='Good')=   drug age sex weight egfr;

where population='Asian';

run;


ObsEffectOddsRatioEstLowerCLUpperCLProbChiSq
1Age1.0470.9941.1030.0839
2Drug A vs B0.7190.2172.3820.5892
3Sex Female vs Male1.9320.4458.3840.3795
4Weight0.9870.9541.0200.4357
5eGFR1.0140.9911.0360.2300


As you can notice, the influence of  drug is significant in the European, but not in the Asian population. One could interpret this as the effect of the drug being modified according to the study subpopulation, and hence expect the p value for interaction to be significant, but it is not at all!


proc logistic data=Studydatabase;

class sex drug study;

model effectnum (event='Good')=  study drug age sex weight egfr population*drug;

run;

Interactin term P value

ProbChiSq

0.1874

I know this forum is mainly for questions regarding the SAS code itself, but I would appreciate any help from all the SAS and statistics experts here

Best wishes

Am

Any explanation please? Am I missing something

1 ACCEPTED SOLUTION

Accepted Solutions
SteveDenham
Jade | Level 19

You had the following code:

proc logistic data=Studydatabase;

class sex drug study;

model effectnum (event='Good')=  study drug age sex weight egfr population*drug;

run;

Now one thing I would do before continuing is to include 'population' in the class statement, which would give:

proc logistic data=Studydatabase;

class sex drug study population;

model effectnum (event='Good')=  study drug age sex weight egfr population*drug;

run;

Now in both of these model statements, there is a variable called 'weight', right there before egfr.  It is not in the class statement, so I assume that it is a continuous covariate, and the Odds Ratio estimate of 0.987 says that the odds of a good event decrease slightly for each one unit increase in weight.

I suggested adding the term weight*population to the model, as I believe the event occurrence relationship to weight may differ in the two populations, and this needs to be removed as a possible confounder of the drug effect.  This would result in the following code:

proc logistic data=Studydatabase;

class sex drug study population;

model effectnum (event='Good')=  study drug age sex weight egfr population*weight population*drug;

run;

Steve Denham

View solution in original post

4 REPLIES 4
SteveDenham
Jade | Level 19

Looking at the overlap in the confidence bounds on the two odds ratios, I am not surprised that the test for the point interaction is not significant.

European (1.229, 2.821)

Asian       (0.219, 2.382)

Note also, that there is a significant weight effect in the European dataset that is not found in the Asian dataset.  This may be a confounded variable.  What happens if you add population*weight to the model?

Steve Denham

ammarhm
Lapis Lazuli | Level 10

Dear Steve

Thank you very much, I have also noticed the overlap in the confidence intervals.

Could you please how to add the *weight term in the code? I mean if I do this I would get an error as I don't have a variable (weight):

proc logistic data=Studydatabase;

class sex drug study;

model effectnum (event='Good')=  study drug age sex weight egfr population*weight;

run;

Thanks again

SteveDenham
Jade | Level 19

You had the following code:

proc logistic data=Studydatabase;

class sex drug study;

model effectnum (event='Good')=  study drug age sex weight egfr population*drug;

run;

Now one thing I would do before continuing is to include 'population' in the class statement, which would give:

proc logistic data=Studydatabase;

class sex drug study population;

model effectnum (event='Good')=  study drug age sex weight egfr population*drug;

run;

Now in both of these model statements, there is a variable called 'weight', right there before egfr.  It is not in the class statement, so I assume that it is a continuous covariate, and the Odds Ratio estimate of 0.987 says that the odds of a good event decrease slightly for each one unit increase in weight.

I suggested adding the term weight*population to the model, as I believe the event occurrence relationship to weight may differ in the two populations, and this needs to be removed as a possible confounder of the drug effect.  This would result in the following code:

proc logistic data=Studydatabase;

class sex drug study population;

model effectnum (event='Good')=  study drug age sex weight egfr population*weight population*drug;

run;

Steve Denham

ammarhm
Lapis Lazuli | Level 10

Dear Steve

Thank you very much, I have also noticed the overlap in the confidence intervals.

Could you please how to add the *weight term in the code? I mean if I do this I would get an error as I don't have a variable (weight):

proc logistic data=Studydatabase;

class sex drug study;

model effectnum (event='Good')=  study drug age sex weight egfr population*weight;

run;

Thanks again

sas-innovate-2024.png

Join us for SAS Innovate April 16-19 at the Aria in Las Vegas. Bring the team and save big with our group pricing for a limited time only.

Pre-conference courses and tutorials are filling up fast and are always a sellout. Register today to reserve your seat.

 

Register now!

What is Bayesian Analysis?

Learn the difference between classical and Bayesian statistical approaches and see a few PROC examples to perform Bayesian analysis in this video.

Find more tutorials on the SAS Users YouTube channel.

Click image to register for webinarClick image to register for webinar

Classroom Training Available!

Select SAS Training centers are offering in-person courses. View upcoming courses for:

View all other training opportunities.

Discussion stats
  • 4 replies
  • 875 views
  • 4 likes
  • 2 in conversation