首先,论坛里面提到了下面这个函数
geometric_mean=geomean(of var1-var5);
但是这个函数是给若干个变量求几何均数的,并不是给一个变量的多个观测求几何均数,无法满足需求。
其次,论坛里大家分享了手工计算,过程太麻烦了,我把程序也贴上,大家看看就行了,不必使用这个。
data test;
set myData;
LogV1 = log(V1);
run;
proc means data=test alpha=0.1;
var LogV1;
output out=myStats
mean=meanLogV1
lclm=lclmLogV1
uclm=uclmLogV1;
run;
proc sql;
select exp(meanLogV1) as geometricMean,
exp(lclmLogV1) as lclmGeoMean,
exp(uclmLogV1) as uclmGeoMean
from myStats;
quit;
【方法一】使用ttest,没错就是做T检验的那个过程,没想到了吧
- proc ttest data=sashelp.class dist=lognormal;
- var age;
- ods select ConfLimits;
- run;
- proc surveymeans data=sashelp.class geomean;
- var age;
- run;