在sas中计算中位数的95%置信区间:
ods exclude all;
/*ods select Quantiles;*/
ods output Quantiles=Quantiles_out BasicMeasures=BasicMeasures_out;
proc univariate data=adexsum(where=(PARAMCD="CI_TIMES")) CIPCTLDF /*CIPCTLNORMAl*/ ciquantdf(alpha=0.05);
var AVAL;
/* output out=median_ci pctlpre=Median pctlpts=50 pctlname=Median_CI range=range min=min max=max;*/
run;
ods exclude close;
通过ods将计算结果输出,Quantiles=Quantiles_out,Quantiles_out指定输出数据集的名称,CIPCTLDF和CIPCTLNORMAl是两种不同的计算方法,若为正态分布则使用CIPCTLNORMAl,否则使用CIPCTLDF,输出结果中以下几列分别对应了正态分布和自由分布的95%CI的上下限。