Documentation Example 2 for PROC TRANSREG
/****************************************************************/
/* S A S S A M P L E L I B R A R Y */
/* */
/* NAME: TREGEX2 */
/* TITLE: Documentation Example 2 for PROC TRANSREG */
/* PRODUCT: STAT */
/* SYSTEM: ALL */
/* KEYS: Box-Cox */
/* PROCS: TRANSREG */
/* DATA: */
/* */
/* SUPPORT: saswfk UPDATE: July 25, 2010 */
/* REF: PROC TRANSREG, EXAMPLE 2 */
/* MISC: */
/****************************************************************/
title 'Yarn Strength';
proc format;
value a -1 = 8 0 = 9 1 = 10;
value l -1 = 250 0 = 300 1 = 350;
value o -1 = 40 0 = 45 1 = 50;
run;
data yarn;
input Fail Amplitude Length Load @@;
format amplitude a. length l. load o.;
label fail = 'Time in Cycles until Failure';
datalines;
674 -1 -1 -1 370 -1 -1 0 292 -1 -1 1 338 0 -1 -1
266 0 -1 0 210 0 -1 1 170 1 -1 -1 118 1 -1 0
90 1 -1 1 1414 -1 0 -1 1198 -1 0 0 634 -1 0 1
1022 0 0 -1 620 0 0 0 438 0 0 1 442 1 0 -1
332 1 0 0 220 1 0 1 3636 -1 1 -1 3184 -1 1 0
2000 -1 1 1 1568 0 1 -1 1070 0 1 0 566 0 1 1
1140 1 1 -1 884 1 1 0 360 1 1 1
;
ods graphics on;
proc transreg details data=yarn ss2
plots=(transformation(dependent) obp);
model BoxCox(fail / convenient lambda=-2 to 2 by 0.05) =
qpoint(length amplitude load);
run;
title 'Atmospheric Pressure Changes Between'
' Easter Island & Darwin, Australia';
ods graphics on;
proc transreg data=sashelp.enso;
model identity(pressure) = pbspline(year);
run;
proc transreg data=sashelp.enso;
model identity(pressure) = pbspline(year / sbc);
run;
proc transreg data=sashelp.enso;
model identity(pressure) = pbspline(year / sbc lambda=2 10000 range);
run;
proc transreg data=sashelp.enso;
model identity(pressure) = pbspline(year / sbc lambda=.1 .5 1 5
10 50 100 500 to 2500 by 500);
run;