R语言怎么求置信区间

#调用程序包
library(sampling)
library(dplyr)
library(magrittr)
library(TeachingSampling)
library(lpSolve)
library(grid)
library(Matrix)
library(survival)
library(survey)

#读取数据
ly=read.table(result.txt,sep=" ",head=TRUE)
#简单随机抽样
N=nrow(ly)
n=20
s=srswor(n,N)
s
data.srswor=getdata(ly,s)
#估计
pw=rep(N/n,n)
fpc=rep(N,n)
agsrswor=as.data.frame(cbind(data.srswor,pw,fpc))
dsrswor<-svydesign(id = ~1,weights = ~pw,data = agsrswor,fpc = ~fpc)

#人均最终成绩，估计量的方差及其在95%的置信水平下的置信区间
svymean(Total,dsrswor,deff=TRUE)
#总最终成绩，估计量的方差及其在95%的置信水平下的置信区间
svytotal(Total,dsrswor,deff=TRUE)

#比率估计（用平时成绩做辅助变量）
#估计人均最终成绩，估计量的方差及其在95%的置信水平下的置信区间
ratio<-svyratio(Total,U_result,dsrswor)
lucym<-data.frame(U_result=mean(ly$U_result))
predict(ratio,lucym$U_result)
#估计总最终成绩，估计量的方差及其在95%的置信水平下的置信区间
ratio<-svyratio(Total,U_result,dsrswor)
lucyt<-data.frame(U_result=sum(ly$U_result))
predict(ratio,lucyt$U_result)

#回归估计（用平时成绩做辅助变量）
#估计人均最终成绩，估计量的方差及其在95%的置信水平下的置信区间
reg<-svyglm(TotalU_result,design=dsrswor)
lucym<-data.frame(U_result=mean(ly$U_result))
predict(reg,newdata = lucym)
#估计总最终成绩，估计量的方差及其在95%的置信水平下的置信区间
reg<-svyglm(TotalU_result,design=dsrswor)
lucyt<-data.frame(U_result=mean(ly$U_result))
predict(reg,newdata = lucyt,total = N)