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; ... View more