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!
Build your skills. Make connections. Enjoy creative freedom. Maybe change the world. Registration is now open through August 30th. Visit the SAS Hackathon homepage.
Register today!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.
Select SAS Training centers are offering in-person courses. View upcoming courses for: