***************************************************************; ** PROGRAM: stability_testing.sas **; ** AUTHOR: Christopher R. Bilder **; ** Department of Statistics **; ** Oklahoma State University **; ** bilder@okstate.edu **; ** www.chrisbilder.com **; ** DATE: 7-25-00 **; ** PURPOSE: Find the estimated regression model for the **; ** stability testing activity. Create a scatter **; ** plot with the estimated regression line, **; ** confidence interval bands, and prediction **; ** interval bands. **; ** NOTES: **; ** 1) This free program is not supported in any way. **; ** 2) Do not distribute copies of this program without **; ** the written permission of the author **; ** 3) Copyright 2000 Christopher R. Bilder **; ** **; ***************************************************************; dm 'log;clear;output;clear;'; options ps=90 ls=65 pageno=1; options compress=no; options mprint symbolgen mlogic; title1 'Stability Testing'; data set1; input time potency; datalines; 3 0.99 6 1.02 9 0.98 12 0.985 15 0.975 18 0.96 21 0.965 24 0.961 30 0.952 36 0.950 48 0.940 60 0.951 run; title2 'The data set'; proc print data=set1; run; title2 'Proc reg output'; proc reg data=set1; model potency = time / clm cli p influence; run; goptions reset=global border ftext=swiss gunit=cm htext=0.4; goptions display noprompt; proc gplot data=set1; title2 'Potency vs. Time'; *plot potency*time twice so that CI and PI bands can be drawn; plot potency*time potency*time / overlay frame haxis=axis1 vaxis=axis2 vref=0.95 vref=1.05 lvref=3; * rlclm95 for CIs; symbol1 v=dot i=rlclm95 height=0.2 ci=red cv=green; * rlcli95 for PIs; symbol2 v=dot i=rlcli95 height=0.2 ci=blue cv=green; axis1 label = ('Time' ) length = 12; axis2 label = (a=90 'Potency' ) length = 7 order = (0.90 to 1.10 by 0.05); run; quit;