source: main/adopters/nj/trunk/src/main/ibisq/qModules3/birth/AvgGestAge.def @ 26928

f label Average&Gestational&Age
f type special
#f xml_out_map_file XMLAverage.map
#f out_detail lbl_not_used__see_xml_out_map_file 
# data_where binfdead^=2
f data_where bgestage<48
#f variable bgestage

--------BoNdArY--------
1 script
OPTIONS MPRINT MLOGIC SYMBOLGEN;
OPTION SPOOL;

proc summary data=tmp;
  var bgestage;
  class %cross1% 
  ?cross2? %cross2%
  ;
  format mean 5.2 count 9.0 sum 9.0 stderr 6.2;
  output out=tmp (drop=_TYPE_ _FREQ_) 
    mean=mean n=count sum=sum stderr=stderr;
 run;
 
 *****************************************************************;
 * Student's t distribution. df=count-1. The "Z" value is really  ;
 * the critical value of "t", two tailed, with an alpha of 0.025. ;
 *****************************************************************;
 data tmp;
        set tmp;
       if count=1  then Z=  .;
  else if count=2  then Z= 12.706;
  else if count=3  then Z=  4.303;
  else if count=4  then Z=  3.182;
  else if count=5  then Z=  2.776;
  else if count=6  then Z=  2.571;
  else if count=7  then Z=  2.447;
  else if count=8  then Z=  2.365;
  else if count=9  then Z=  2.306;
  else if count=10 then Z=  2.262;
  else if count=11 then Z=  2.228;
  else if count=12 then Z=  2.201;
  else if count=13 then Z=  2.179;
  else if count=14 then Z=  2.160;
  else if count=15 then Z=  2.145;
  else if count=16 then Z=  2.131;
  else if count=17 then Z=  2.120;
  else if count=18 then Z=  2.110;
  else if count=19 then Z=  2.101;
  else if count=20 then Z=  2.093;
  else if count=21 then Z=  2.086;
  else if count=22 then Z=  2.074;
  else if count=23 then Z=  2.069;
  else if count=24 then Z=  2.064;
  else if count=25 then Z=  2.060;
  else if count=26 then Z=  2.056;
  else if count=27 then Z=  2.052;
  else if count=28 then Z=  2.048;
  else if count=29 then Z=  2.045;
  else if count=30 then Z=  2.042;
  else if count>30 then Z=  1.96;

  t1=mean-(Z*stderr);
  if (t1<0) then t1=0;
  if count in (0 1) then t1=.;
  LL=put(t1, 8.2);
  UL=put((mean+(z*stderr)), 8.2);
  LL=compress(LL);
  UL=compress(UL);
  n=count;   *ibis-q needs a count variable named 'n';
  run;
  
/* data tmp;
  set tmp;
  if (0<n<20) then do;
   mean = .A;
   LL = put('**', 8.0);
   UL = put('**', 8.0);
  redflag = put('Suppressed', $12.);
   end;
   run; */
 proc print data=tmp noobs;  title2 '6. TMP - final dataset to pass to IBIS View app'; 
 run;
 
--------BoNdArY--------
f out_variable average
#########################################
# definition for output file
f xml_out_map_file XMLAverage.map
 --------BoNdArY--------
f out_detail lbl_not_used__see_xml_out_map_file
  mean 15.2
  LL 15.2
  UL 15.2
  sum 15.0
  n 15.2
--------BoNdArY--------

--------BoNdArY--------
1 suppressed_variables 1
mean
--------BoNdArY--------