1.求总体均数的95%置信区间,可根据自身需要更该alpha值，计算相应的置信区间

--------------------------------------------------------------------------------------
##总体均数95%置信区间的计算alpha<-0.05                        #置信区间 n <- length(data)                  #计算变量长度mean<- mean(data)                  #计算均值Z1 <- qt(1-alpha/2,n-1)          #计算对应的Z值（双侧）Z2 <- qt(1-alpha,n-1)            #计算对应的Z值（单侧）C1 <- round(mean-Z1*sqrt(sum((data-mean)^2)/(n-1))/sqrt(n),2)   #双侧下限，取两位小数C2 <- round(mean+Z1*sqrt(sum((data-mean)^2)/(n-1))/sqrt(n),2)   #双侧上限，取两位小数C3 <- round(mean-Z2*sqrt(sum((data-mean)^2)/(n-1))/sqrt(n),2)   #单侧下限，取两位小数C4 <- round(mean+Z2*sqrt(sum((data-mean)^2)/(n-1))/sqrt(n),2)   #单侧上限，取两位小数CI1 <- c(C1,C2)          #95%CI双侧CI2 <- c(C3,C4)          #95%CI单侧
---------------------------------------------------------------------------------------

2.几何均数（GMT）及其置信区间

几何均数95%置信区间的计算类似于总体均属95%置信区间的计算，只是在计算几何均属95%置信区间时，须提前将数据进行对数转化（采用log()函数进行对数转换），经上述步骤计算经对数转换后的95%置信区间，再接着取反对数，如下述代码，即可求得结果。

--------------------------------------------------------------------------------------
几何均数（GMT）的计算
#方法1：
library("psych")
geo_mean<-geometric.mean(data[,x])#方法2：
geo_mean <- exp(mean(log(data[,x])))--------------------------------------------------------------------------------------
几何均数95%置信区间的计算geometry.CI1 <-exp(C1)
geometry.CI2 <-exp(C2)结果即为：geo_mean(geometry.CI1 ,geometry.CI2 )