Thanks for the response, it was really helpful! Based on what you've told me and stuff I've read elsewhere, I'm planning on reporting the geometric mean (gm) and the standard deviation (sd) for gm. I derived these two by exponentiating the mean and sd for my log transformed variable, "ln_vit." Essentially, exponentiating both mean & sd of "ln_vit" = gm & sd(of gm) for "vit", proc means data=work.run1;
var vit ln_vit;
class disease;
run; My last question I would appreciate help with is doing ANOVA for my skewed variable, "vit," by disease severity (i.e. three categories - mild, moderate, & severe). My guess is that I do ANOVA based on "ln_vit" (the log transformed values for "vit"). Correct me if I'm wrong - the p value I get for this ANOVA would also be the p value if I compared gm of vit, by disease severity? PROC GLM DATA=work.run1 PLOTS=none;
CLASS disease_severity;
MODEL ln_vit=disease_severity;
MEANS disease_severity;
RUN; Side note - Here are the results you suggested I post, regarding my initial questions:
... View more