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!
SAS Innovate 2025 is scheduled for May 6-9 in Orlando, FL. Sign up to be first to learn about the agenda and registration!
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.
Ready to level-up your skills? Choose your own adventure.