topic Re: Multivariate analysis: Do I use Proc Corr or Proc GLM or Proc Reg? in SAS Programming

Multivariate analysis: Do I use Proc Corr or Proc GLM or Proc Reg?

hpatel3 — Thu, 18 Oct 2018 14:19:35 GMT

Hello!

I have a project at work where I'm being asked to analyze 7 variables with (5 of them being categorical(Yes/No) and 2 being numerical)

and their correlation with disease status (Variable ADV_HF): 1= they have the disease / 0 or blank = they don't have it.

I have not used multivariate analysis before and the different types are a little overwhelming.

Based on the SAS forums, I am under the impression that I shouldn't use Proc Reg since I have categorical variables, so should I use Proc GLM or Proc Corr? Will it make a huge difference?

What I've got so far is:

proc glm data= Hetal.ES_Regression;

Class Adv_HF;

model age_diag fam_Hx_ES fam_Hx_SD hx_sync LBBB EF_Reg AF_prior ;

1.I'm not sure how or if to use the contrast statement and manova statement.

2. How do I specify in the class statement that 1= disease state and 0- without disease?

3. Am I missing anything other key data step in this analysis?

Thank you!

Re: Multivariate analysis: Do I use Proc Corr or Proc GLM or Proc Reg?

PaigeMiller — Thu, 18 Oct 2018 14:42:35 GMT

I have a project at work where I'm being asked to analyze 7 variables with (5 of them being categorical(Yes/No) and 2 being numerical) and their correlation with disease status (Variable ADV_HF): 1= they have the disease / 0 or blank = they don't have it.

If the specific request you have is simply to analyze correlations, then you would use PROC CORR.

If the underlying reason is to fit a model, you should use PROC LOGISTIC (which is appropriate when your response variable is binary).

This is not (at least the way I use the word) a "multivariate" analysis, and no MANOVA would work here anyway. Multivariate would imply to me that you have multiple response variables, and if the multiple response variables are continuous, that is the only time when any MANOVA would work. So none of this applies to your situation, as I understand it.

Re: Multivariate analysis: Do I use Proc Corr or Proc GLM or Proc Reg?

Reeza — Thu, 18 Oct 2018 14:57:29 GMT

A really good tutorial is this one:
https://stats.idre.ucla.edu/sas/dae/logit-regression/

Note the use of the PARAM=REF option on the CLASS statement. You will want to do that. Additionally, check this example out:

https://documentation.sas.com/?docsetId=statug&docsetTarget=statug_logistic_examples02.htm&docsetVersion=14.3&locale=en

Re: Multivariate analysis: Do I use Proc Corr or Proc GLM or Proc Reg?

hpatel3 — Thu, 18 Oct 2018 15:58:38 GMT

Hi @PaigeMiller , thanks for responding.

I decided to use Proc Corr with the following code:

and it keeps giving me this error in my log:

Why does it keep telling me my variables do not match the type prescribed for this list? What am i doing wrong here

Do i need to denote which are categorical? Also do i need to use the "BY" statement to classify that I want these variables compared with Those who have the disease (1) vs those who dont (0)?

Re: Multivariate analysis: Do I use Proc Corr or Proc GLM or Proc Reg?

hpatel3 — Thu, 18 Oct 2018 16:00:58 GMT

Thanks for replying @Reeza!

I don't think Proc logistic works in this case because we're not looking for a specific question/outcome. Merely seeing if there is a correlation between the variables whether they have or don't have the disease.

My variables are potential risk factors and we want want to see if there's any correlation between these and disease status.

So I think I'd use proc Corr, yes?

Re: Multivariate analysis: Do I use Proc Corr or Proc GLM or Proc Reg?

Reeza — Thu, 18 Oct 2018 16:09:36 GMT

Nope, you have a binary outcome variable so using PROC CORR is not suitable here. You really are looking for logistic regression here and the odds ratios, whether you do it one variable at a time or a full model.

@hpatel3 wrote:

Thanks for replying @Reeza!

I don't think Proc logistic works in this case because we're not looking for a specific question/outcome. Merely seeing if there is a correlation between the variables whether they have or don't have the disease.

My variables are potential risk factors and we want want to see if there's any correlation between these and disease status.

So I think I'd use proc Corr, yes?

Re: Multivariate analysis: Do I use Proc Corr or Proc GLM or Proc Reg?

Reeza — Thu, 18 Oct 2018 16:10:45 GMT

And your issue with CORR is that it needs numeric variables so you should fix your data first and ensure that numeric variables are numbers. Categorical variables can be either numeric or categorical.

Re: Multivariate analysis: Do I use Proc Corr or Proc GLM or Proc Reg?

PaigeMiller — Thu, 18 Oct 2018 16:17:53 GMT

PROC CORR can produce a correlation-like number called Kendall's Tau for categorical variables. However, they probably do need to be converted to category numbers (a numeric variable) in order for PROC CORR to process them.

Re: Multivariate analysis: Do I use Proc Corr or Proc GLM or Proc Reg?

hpatel3 — Thu, 18 Oct 2018 18:42:55 GMT

It worked! I got results! (I Think)

My only question left is about status. Does multivariate analysis take into account that for developing the disease (event status=1) vs not developing the disease (event status=0) need to be specified anywhere? Or does SAS automatically assume 0=no event and 1=event?? My code and results are as followed:

Re: Multivariate analysis: Do I use Proc Corr or Proc GLM or Proc Reg?

Reeza — Thu, 18 Oct 2018 19:07:04 GMT

For a univariate analysis with a single binary outcome, I would recommend ANOVA\T-Test for the continuous variables and chi square for categorical variables.

For your output, you have Pearson correlation coefficients and SAS makes no assumptions regarding the value of 0/1 being a particular record type. You probably want Tau-b or Tau-c instead.

https://en.wikipedia.org/wiki/Kendall_rank_correlation_coefficient

Re: Multivariate analysis: Do I use Proc Corr or Proc GLM or Proc Reg?

hpatel3 — Fri, 19 Oct 2018 17:33:37 GMT

@PaigeMiller , I converted to categorical variables and I do see Kendall's tau, but I figured pearson's would be the number to look at. Why do you suggest Tau instead?

Should I be using pearson's on the continuous variables and Tau on the categorical ones?

Re: Multivariate analysis: Do I use Proc Corr or Proc GLM or Proc Reg?

PaigeMiller — Fri, 19 Oct 2018 18:29:48 GMT

Yes, as @Reeza pointed out, and I omitted, you want Kendall's Tau-B (or maybe it's a lowercase b) which is appropriate for computing a correlation-like measure when you have two ordinal variables. You use Pearson when you have two continuous variables.

I still have issues with your requirement that you want a measure like correlation when you have predictor (X) variables and response (Y) variables, correlation is not meant for that case; some measure of how well the X predict the Y is the appropriate statistic.

I can't agree with this statement from Reeza

For a univariate analysis with a single binary outcome, I would recommend ANOVA\T-Test for the continuous variables

This isn't correlation, which seems to be what you are asking for, although I don't understand why; and it also seems to reverse the role of X and Y. You don't do ANOVA or t-tests with binary Y, you do it for binary X. For continuous X variables, and binary Y, logistic regression is still what I would use, and the measure you want is the odds ratio or the slope of the logistic regression.

Re: Multivariate analysis: Do I use Proc Corr or Proc GLM or Proc Reg?

Reeza — Fri, 19 Oct 2018 19:26:31 GMT

I would treat the 0/1 outcome as groups and checking for difference between the groups. So is the mean for the continuous variables different for those who have disease X and those who don’t. If that’s not appropriate for a t-test/ANOVA, what test would you use?

Re: Multivariate analysis: Do I use Proc Corr or Proc GLM or Proc Reg?

PaigeMiller — Fri, 19 Oct 2018 20:04:19 GMT

As I stated, you seem to have reversed X and Y. Logistic regression is what you use with binary Y.

But, the whole issue remains unclear as to what the original poster really wants to achieve, and so I think until he clarifies the situation, I'm going to pause here. On the one hand he wants correlation but on the other hand he was talking about fitting a model with PROC REG or PROC GLM.

Re: Multivariate analysis: Do I use Proc Corr or Proc GLM or Proc Reg?

hpatel3 — Mon, 22 Oct 2018 14:34:35 GMT

@PaigeMiller @Reeza

After going back and reviewing, I don't think Proc reg or corr would give me the best results. I stated originally that I want correlation between variables and I still do. I thought Proc Corr/GLM/REG would give me that. That's why I was trying to use them. But someone on another thread stated that Proc Logistic would work for both continuous and categorical variables, so if I can just use that, that would work I believe? I don't need to fit the model or anything I believe.

So TO SEE IF THERE IS A CORRELATION between my risk factors and whether they develop the disease, I tried this:

libname Hetal "\\tuftsmc\home\hpatel3\SAS Datasets";
run; 



proc logistic data=hetal.es_regre;
class age_diag EF_reg Hx_Sync FHX_ES FHX_SD LBBB AF_prior * Hx_Sync FHX_ES FHX_SD LBBB AF_prior;
model adv_hf (event="1")= age_diag EF_REG Hx_sync FHX_ES FHX_SD LBBB AF_prior;
run;

However got this error:

I'm not sure exactly what should be done based on this error. Is asking me to clarify numeric or character variables??? I'm not sure how to write this out in my code.

My class statement is: all variables * categorical variables <-- That is what I was supposed to do, correct?

And my model statement is: disease status variable = all variables <--is this correct as well?

Thanks!

Re: Multivariate analysis: Do I use Proc Corr or Proc GLM or Proc Reg?

Reeza — Mon, 22 Oct 2018 14:40:19 GMT

and their correlation with disease status (Variable ADV_HF): 1= they have the disease / 0 or blank = they don't have it.

I interpreted this as the outcome, or Y variable, which is binary to me. I believe you're correct with the OP not understanding what they want to do though.

@PaigeMiller wrote:

As I stated, you seem to have reversed X and Y. Logistic regression is what you use with binary Y.

But, the whole issue remains unclear as to what the original poster really wants to achieve, and so I think until he clarifies the situation, I'm going to pause here. On the one hand he wants correlation but on the other hand he was talking about fitting a model with PROC REG or PROC GLM.

Re: Multivariate analysis: Do I use Proc Corr or Proc GLM or Proc Reg?

PaigeMiller — Mon, 22 Oct 2018 14:59:45 GMT

You can't use an asterisk in a CLASS statement.

In the future, please click on the {i} icon and paste your SASLOG into the window that appears.

I stated originally that I want correlation between variables and I still do. I thought Proc Corr/GLM/REG would give me that. That's why I was trying to use them. But someone on another thread stated that Proc Logistic would work for both continuous and categorical variables, so if I can just use that, that would work I believe? I don't need to fit the model or anything I believe.

Unfortunately, I still believe that this is not a good approach, in my opinion you DO need a model, you don't want correlations. Why? Because a model treats your Y variable (presence or absence of the disease) as the response. Correlations do not treat the presence or absence of the disease as the response.

But ... whether to use a model or not ... whether to use correlations ... that's your decision, not mine.

Finally, let's clear things up once and for all, PROC LOGISTIC works with class Y variables, and continuous or class X variables. So, it would seem to fit your situation, except you still state you don't want a model, and PROC LOGISTIC gives you a model. And it has to have a Y variable and one or more X variables.

Re: Multivariate analysis: Do I use Proc Corr or Proc GLM or Proc Reg?

Reeza — Mon, 22 Oct 2018 14:49:26 GMT

If you don't need a model, why are you running a model?

I really really feel you need to consult a statistician here because you're not sure what you're doing and you're definitely not quite understanding the appropriate analysis techniques. Or ask a teacher, colleague or supervisor for some assistance. Or read a journal article that does this and see how they do the analysis and present the data. The typical approach is formulate a hypothesis, do univariate analysis, then do multivariate analysis.

@hpatel3 wrote:

@PaigeMiller @Reeza

After going back and reviewing, I don't think Proc reg or corr would give me the best results. I stated originally that I want correlation between variables and I still do. I thought Proc Corr/GLM/REG would give me that. That's why I was trying to use them. But someone on another thread stated that Proc Logistic would work for both continuous and categorical variables, so if I can just use that, that would work I believe? I don't need to fit the model or anything I believe.

So TO SEE IF THERE IS A CORRELATION between my risk factors and whether they develop the disease, I tried this:
libname Hetal "\\tuftsmc\home\hpatel3\SAS Datasets";
run; 



proc logistic data=hetal.es_regre;
class age_diag EF_reg Hx_Sync FHX_ES FHX_SD LBBB AF_prior * Hx_Sync FHX_ES FHX_SD LBBB AF_prior;
model adv_hf (event="1")= age_diag EF_REG Hx_sync FHX_ES FHX_SD LBBB AF_prior;
run;
However got this error:

I'm not sure exactly what should be done based on this error. Is asking me to clarify numeric or character variables??? I'm not sure how to write this out in my code.

My class statement is: all variables * categorical variables <-- That is what I was supposed to do, correct?

And my model statement is: disease status variable = all variables <--is this correct as well?

Thanks!