topic Re: Assistance with adjusted geometric mean in Statistical Procedures

Assistance with adjusted geometric mean

leaning2sas — Tue, 31 Jan 2023 17:05:36 GMT

Hello!

I am a basic user and have been trying to get a geometric mean of the variable of interest (concentration of a chemical) in the urine stratified by different years and age. Since this chemical is present in urine, it must be adjusted for creatinine. The variable of interest has been converted into its log, as it is not normally distributed. I used proc glimmix as this was nonnormal data and adjusted the model for creatinine, to obtain the predicted mean, and then used those to get the geometric means. I am not sure if my code is right or even if I am using the right approach. I would appreciate it if you could guide me regarding the same.

proc glimmix data = NH.demo0118_phth_fv2;
class timephth;
model mbp = timephth creatinine_phth/solution; * URXUCR=urinary_creatinine;
output out=NH.ph_ad pred=p_mbp_mean;
run;

proc surveymeans data = NH.ph_ad allgeo;
domain age timephth;
var p_mbp_mean;
title "Geometric mean of XX years of MBP";
ODS output domaingeomeans = NH.mbp_gm;
run;

Re: Assistance with adjusted geometric mean

SteveDenham — Wed, 01 Feb 2023 16:30:33 GMT

If you want an adjusted geometric mean, I would suggest the following code:

proc glimmix data=NH.demo0118_phth_fv2;
ods output lsmeans=geomeans;
class timephth;
model value=timephth creatinine_phth/dist=lognormal;
lsmeans timephth/cl;
run;

data geomeans2;
set geomeans;
mean=exp(estimate);
run;

This will give you the geometric mean at the ARITHMETIC mean of creatinine_phth. You can also get confidence bounds by exponentiating the variables lower and upper in the geomeans dataset.

SteveDenham

Re: Assistance with adjusted geometric mean

Rick_SAS — Thu, 02 Feb 2023 10:43:29 GMT

In case you wonder why Steve's suggestion works, recall that the geometric mean for a set of positive numbers is equivalent to computing the arithmetic mean of the log-transformed data and then using the exponential function to back-transform the result. See https://blogs.sas.com/content/iml/2019/09/30/what-is-a-geometric-mean.html

Re: Assistance with adjusted geometric mean

leaning2sas — Fri, 03 Feb 2023 04:02:02 GMT

Thank you very much! I appreciate your assistance with the code. I was confused while writing the code, but now I understand.

Re: Assistance with adjusted geometric mean

leaning2sas — Fri, 03 Feb 2023 04:03:44 GMT

I appreciate the explanation. The link was very helpful in understanding geometric mean since it is not used often in data analysis.

Re: Assistance with adjusted geometric mean

leaning2sas — Fri, 03 Feb 2023 05:24:17 GMT

Just a follow-up question, since I'm using a survey dataset, can the same codes be used? Also, how do I stratify by age? For e.g., I want creatinine adjusted geometric mean for age groups 6-12 and 13-18. I really appreciate any help you can provide. Thanks in advance!

Re: Assistance with adjusted geometric mean

SteveDenham — Thu, 09 Feb 2023 18:29:53 GMT

If the data come from a survey, you may want to change to one of the appropriate survey procedures such as SURVEYMEANS or SURVEYREG. SURVEYMEANS can provide geometric means directly, while SURVEYREG would require transformation/back-transformation. Shifting to the SURVEY procedures requires you to have knowledge of the complex survey design used to collect the data.

There are several ways to convert continuous variables such as age to categorical variables in a DATA step, but beware of loss of power/increase in standard error when doing so. There is a lot of discussion in the literature about the dangers of categorization of continuous variables. From the short example you present, you would have to convince me that the difference between a 12 year old and a 13 year old (different categories) was markedly greater than between a 12 year old and a 6 year old (same category).

SteveDenham