用r求置信区间,r语言置信区间代码

setwd(‘D:/fc’);x=read.table(‘2.6.1.txt’);x=x$V1;#读数据
n=length(x);alpha=0.5;#置信水平
p0=0.5;#分位数
q0=1200;#要检验的数值
sn=sum(x<q0);#统计量

#符号检验，双边
binom.test(sn,n,p0,alt=c(‘two.sided’))
#或
pvalue=2*min(pbinom(sn,n,p0),1-pbinom(sn-1,n,p0));pvalue

#符号检验，左边（原假设大于等于）
quantile(x,p0);
binom.test(sn,n,p=p0,alt=‘greater’)
#或
pvalue=1-pbinom(sn-1, n,p0);pvalue

#符号检验，右边（原假设小于等于）
quantile(x,p0);
binom.test(sn,n,p=p0,alt=‘less’)
#或
pvalue=pbinom(sn, n,p0);pvalue

#符号检验分位数的置信区间
k1=qbinom(alpha/2,n,p0);if (pbinom(k1,n,p0)==alpha/2) k1=k1+1;
k2=qbinom(1-alpha/2,n,p0)-1;
sx=sort(x);L=sx[k1];U=sx[k2+1];list(lower=L,upper=U,conf.level=1-alpha)