data jdp;
 input life @@;
 cards;
1993 1993 1984 1997 1995 1991 1984 1996 1989 2001
;
run;
proc summary data=jdp;
 var life;
 output out=SummaryStats n=n mean=PointPredictor std=StdDev;
run;
data preds;
 set SummaryStats;
 df=n-1;
 L=.90;
 Q=(1+L)/2;
 T_Quantile=tinv(Q,df);
 StdErrPrediction=StdDev*sqrt(1+1/n);
 LowerPredLimit=PointPredictor-T_Quantile*StdErrPrediction;
 UpperPredLimit=PointPredictor+T_Quantile*StdErrPrediction;
 keep L df PointPredictor StdErrPrediction LowerPredLimit UpperPredLimit;
run;
proc print data=preds noobs;run;