Examples of the data that you currently have would be helpful.
Geometric means normally involve the natural log. Are you sure that Log10 is needed?
An example:
data heart;
set sashelp.heart;
lChol = log(Cholesterol);
label lChol="Log(Cholesterol)";
run;
proc genmod data=heart;
class smoking_status;
model lChol = smoking_status;
lsmeans smoking_status / diff exp cl;
run;
I used Proc genmod on some dummy sas-data preparing for real data to calculate GMT, and GMT ratios. I want to compare groups 1-3 separately vs group 4 as regards the ratios and the results are below but I cannot find a way to output the lsmeans gr / diff exp cl results into a dataset for further processing (add to table outputs). Any suggestions? Also the procedure used for the would be appreciated (guess I have to investigate the SAS documentation):
proc genmod data=gmt;
class gr;
model lChol = gr;
lsmeans gr / diff exp cl;
*output out=out_gmt ;
run;
Add the statement
ods trace on;
before the procedure call to get a list of available output datasets in the log. Then use a statement such as
ods output LSMeans=myLsmeans Diffs=myDiffs;
among the procedure statements to output the datasets.
@HF_clin wrote:
... guess I have to investigate the SAS documentation...
definitely!
Are you ready for the spotlight? We're accepting content ideas for SAS Innovate 2025 to be held May 6-9 in Orlando, FL. The call is open until September 25. Read more here about why you should contribute and what is in it for you!
Learn the difference between classical and Bayesian statistical approaches and see a few PROC examples to perform Bayesian analysis in this video.
Find more tutorials on the SAS Users YouTube channel.