程序如下：
data npk;/*建立数据集npk*/
input n p k w;/*建立变量n、p、k、w*/
nn=n*n;/* 建立变量nn其值为n乘n*/
pp=p*p;kk=k*k;np=n*p;nk=n*k;pk=p*k;
cards;/*以下是数据体*/
0 196 372 15.18
34 196 372 21.36
67 196 372 25.72
101 196 372 32.29
135 196 372 34.03
202 196 372 39.45
259 196 372 43.15
336 196 372 43.36
404 196 372 40.83
471 196 372 30.75
259 0 372 33.46
259 24 372 32.47
259 49 372 36.06
259 73 372 37.96
259 98 372 41.04
259 147 372 40.09
259 196 372 41.26
259 245 372 42.17
259 294 372 40.36
259 342 372 42.73
259 196 0 18.98
259 196 47 27.35
259 196 93 34.86
259 196 140 38.52
259 196 186 38.44
259 196 279 37.73
259 196 372 38.43
259 196 465 43.87
259 196 558 42.77
259 196 651 46.22
;/*数据体结束*/
proc reg data=npk;/*调用reg过程用数据集npk中数据*/
model w=n p k nn pp kk np nk pk/selection=stepwise;/*回归模型以w为因变量，以n p k nn pp kk np nk pk为自变量；逐步选入预报因子*/
run