BookmarkSubscribeRSS Feed
tatami
Fluorite | Level 6

I am running MMRM on log10 transformed outcome. There is some medical reasoning for log10 transformation instead of natural log. I would like to model using SAS PROC MIXED. To get geometric mean and geometric ratio, I would need to reverse the transformation. I understand that direct reverse transformation gives median instead of the expected value (Solved: Non-normal data; PROC GLIMMIX - SAS Support Communities). How should I go about this? I would like to get GMR of treatment versus control at a specific time point.

 

proc mixed data=final2;
class id trt avisit;
model log10aval=trt avisit log10base trtp*avisit / S;
repeated avisit / subject=id type= UN R;
LSMEANS trt*avisit / PDIFF CL;
estimate "Visit 9 - Week 52" trt 1 -1 trt*avisit 0 0 0 0 1 0 0 0 0 -1 / CL;
run;

 

 

1 REPLY 1
SteveDenham
Jade | Level 19

If there is a reason for using a log10 transform for this field, then the backtransformed (10**estimate, 10**stderr) values are what you should probably report to be consistent with the literature.  If you truly need the expected value and the standard deviation of the expected value (standard error), there are two methods that we have used.  These are the delta method and the omega method.  The latter is spelled out in the PROC GLIMMIX documentation, referring to a lognormal distribution.  You can use the formulas there, after converting the log10 values to natural logs (multiply by a constant).  For the delta method, it turns out that Wikipedia has good formulas for E(X) and Var(X).

 

SteveDenham

sas-innovate-2024.png

Available on demand!

Missed SAS Innovate Las Vegas? Watch all the action for free! View the keynotes, general sessions and 22 breakouts on demand.

 

Register now!

What is ANOVA?

ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.

Find more tutorials on the SAS Users YouTube channel.

Discussion stats
  • 1 reply
  • 1023 views
  • 0 likes
  • 2 in conversation