R语言 总体均数的95%置信区间、几何均数（GMT）及其置信区间的计算

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单侧
---------------------------------------------------------------------------------------<