不用等一天,看完我这个你可以去睡觉了。
data c;
input y;
cards;
12.15
12.12
12.01
12.28
12.09
12.03
12.01
12.11
12.06
12.14
;run;
proc means noprint data=Work.C;
var Y;
output out=work._ht1v n=_nobs_ mean=_mean_ std=_std_ var=_var_;
run;
data work._vstat_;
set work._ht1v end=_eof_ ;
title1 "实验结果";
keep _nobs_ _mean_ _var_ _std_ _chisq_ _df_ _prob_;
* --- calculate chi-square statistic and probability --- *;
_df_ = _nobs_ - 1;
* --- calculate confidence intervals --- *;
keep _lower_ _upper_;
_lower_ = _var_*_df_/cinv(1-(1-0.95)/2,_df_);
_upper_ = _var_*_df_/cinv( (1-0.95)/2,_df_);
file print;
* --- print results --- *;
put // @5 "关于样本的情况" //
@30 '方差' /
@30 _var_ best6. //;
* --- print confidence intervals --- *;
_ciprob_ = 0.95 * 100;
put // @5 _ciprob_
+(-1) "% 置信区间";
put / @20 '下限' @35 '上限' /
@20 11*'-' @35 11*'-';
put @23 _lower_ best6. @38 _upper_ best6. ;
run;
title1;
proc delete data=work._ht1v; run;