topic Re: adjusting proc glm for ANOVA in SAS Procedures

adjusting proc glm for ANOVA

Manije72 — Mon, 01 Apr 2024 17:28:24 GMT

I would appreciate it if you could answer my question quickly.

I want to see the dietary iron among case (patient) and control (healthy). The variable for disease is risk factor. I want to adjust sex, alcohol (yes/no), drug (yes/no), and BMI. Is this the correct code?

proc glm data=SS;

class risk factor sex alcohol drugs;

model Iron= risk factor sex alcohol drugs BMI_A / solution;

lsmeans risk factor/ stderr pdiff cov out=adjmeans;

run;

Re: adjusting proc glm for ANOVA

PaigeMiller — Mon, 01 Apr 2024 20:04:29 GMT

@Manije72 wrote:

Hi

I would appreciate it if you could answer my question quickly.

I want to see the dietary iron among case (patient) and control (healthy). The variable for disease is risk factor. I want to adjust sex, alcohol (yes/no), drug (yes/no), and BMI. Is this the correct code?

proc glm data=SS;

   class risk factor sex alcohol drugs;

   model Iron= risk factor sex alcohol drugs BMI_A / solution;

   lsmeans risk factor/ stderr pdiff cov out=adjmeans;

run;

If you want to compare the dietary iron between case and control, this is a categorical variable which I do not see in your model.

Also "risk factor" is not a valid variable name in SAS, so it can't be used like that (although risk_factor is a valid variable name).

Re: adjusting proc glm for ANOVA

Manije72 — Mon, 01 Apr 2024 20:57:10 GMT

Thank you for your response.

I named "patients/healthy) as "risk factor" in my data. It has 0(healthy) and 1(patients) code.

Re: adjusting proc glm for ANOVA

Manije72 — Mon, 01 Apr 2024 21:01:22 GMT

Hi,

Thank you for your response.

I named case/control group "risk_factor" in my data. It has 0(control) and 1(case) codes.

I used "risk_factor" in SAS.

Re: adjusting proc glm for ANOVA

PaigeMiller — Mon, 01 Apr 2024 21:49:14 GMT

When you have a variable name in your data such as risk_factor, please use that variable name in your text as well, possibly with other descriptions. Do not make us guess what variable you are referring to.

With that understanding, your code looks reasonable to me.