/* woolybear caterpillar cocoon example */
data woolybear;
   input day cocoontemp airtemp;
   difference=cocoontemp-airtemp;
cards;
1  15.1  10.4
2  14.6   9.2
3   6.8   2.2
4   6.8   2.6
5   8.0   4.1
6   8.7   3.7
7   3.6   1.7
8   5.3   2.0
9   7.0   3.0
10   7.1   3.5
11   9.6   4.5
12   9.5   4.4
;
proc print data=woolybear;
title 'woolybear cocoon temperature example';
/* ------------------------------------------ */
/* graphical summary */
/* ------------------------------------------ */
proc sgplot data=woolybear;
   title 'woolybear temperature difference distribution';
   histogram difference / binwidth=.5 datalabel=percent;
   density difference / type=kernel;
proc sgplot data=woolybear;
   title 'woolybear temperature difference distribution';
   histogram difference / binwidth=1 datalabel=percent;
   density difference / type=kernel;    
proc sgplot data=woolybear;
   title 'woolybear temperature difference distribution';
   hbox difference;
run;
title;
/* ------------------------------------------ */
/* numerical summary */
/* ------------------------------------------ */
proc means data=woolybear maxdec=4
   n min q1 median q3 max range qrange mean std;
/*p1 p5 p10 p20 p30 p40 p50 p60 p70 p80 p90 p95 p99;*/
   var difference;
proc means data=woolybear maxdec=4
/*   n min q1 median q3 max range qrange mean std; */
   p1 p5 p10 p20 p30 p40 p50 p60 p70 p80 p90 p95 p99;
   var difference;
proc univariate data=woolybear nextrval=5;
   var difference;
   ods select ExtremeValues BasicMeasures Quantiles;
   title 'woolybear temperature difference distribution';
run;
title;
/* -------------------------------------------------- */