今天遇到了一个计算几何均数的问题,在论坛里面搜了一下发现没有很好的解决办法,我就去研究了一下,发现了两个简便的计算方法,与大家分享一下。
首先,论坛里面提到了下面这个函数
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检验的那个过程,没想到了吧
【方法二】使用surveymeans过程
结果如下图所示