R/qci.R

In bokov/powertrip: Analysis of longevity data.

`qci` <-
function(x,p=.5,cn=.95,showqs=F){
	x<-sort(x);
	cn<-1.01+cn; n<-length(x); np<-n*p; er<-cn*sqrt(np*(1-p));
	lb<-ceiling(np-er); ub<-ceiling(np+er);
	lb[lb<1]<-1; ub[ub>n]<-n;
	#print(np); print(er); print(lb); print(ub);
        if(showqs) qs<-quantile(x,p) else qs<-NULL
	cbind(x[lb],qs,x[ub]);
}