yydy441.sas

来自「高老师的经典多元统计分析教材源代码非常经典使用」· SAS 代码 · 共 48 行

/* YYDYFX441.sas  */


data d431;
  input year x1-x5 y1 y2;
  cards;
1949   0.9      0.8     0.14     6.63     0.24      1.47      7.31
1950   1.0      2.1     0.15     7.07     0.46      1.25      7.42
1951   2.9      6.3     0.33     7.60     1.02      2.05     11.13
1952   5.0      4.4     0.78    12.88     1.61      2.49     16.08
1953   8.2     13.3     1.18    15.86     1.63      3.16     22.86
1954  13.1     16.8     1.56    18.79     1.93      3.87     29.52
1955  23.8     17.8     2.11    14.63     2.31      4.50     34.54
1956  34.8     27.8     3.09    19.79     3.32      6.09     41.22
1957  35.4     22.1     3.58    16.50     4.44      6.78     47.54
1958  47.0     32.2     7.31    26.22     7.18     10.73     60.00
1959  62.6     33.2     9.61    28.00     8.77     17.65     78.00
1960  68.0     55.6    12.85    27.56     9.89     26.84     96.20
1961  35.3     24.4     6.76    10.95     5.58     24.20     52.37
1962  31.3     17.9     5.08    10.15     6.03     20.08     37.77
1963  35.2     24.8     5.54    14.23     7.18     19.28     40.07
1964  45.3     37.8     7.14    20.38     8.80     22.89     50.36
1965  49.5     78.8    11.20    26.56    10.45     28.94     65.33
1966  59.7    101.6    15.89    33.18    12.51     39.05     83.64
1967  47.8     74.9    10.86    23.90    11.42     39.09     68.16
1968  17.7     40.2     5.10    17.56     9.03     26.81     41.64
1969  36.0     73.3    13.14    27.20     8.05     37.19     67.30
1970  62.0    138.6    25.54    36.28    10.30     54.09    103.57
1971  97.0    247.0    31.31    41.53    14.18     77.39    135.80
1972  95.2    270.0    28.79    40.24    15.19     84.02    118.10
1973 118.4    233.5    28.03    38.20    15.77     88.39    119.62
1974  99.9    205.0    26.50    31.54    12.29     86.32    112.39
1975 151.0    288.0    38.61    46.87    17.36    107.94    144.41
1976 108.0    262.2    31.46    38.62    15.10    102.76    130.66
1977 162.5    358.6    46.21    52.48    20.48    118.84    175.10
1978 238.2    454.8    55.86    55.96    26.40    139.30    214.44
 ;
 proc print;
 run;
 proc reg data=d431;
    model y1 y2=x1-x5/selection=stepwise
                      sle=0.05 sls=0.05;
run;
 proc reg data=d431;
    model y1 y2=x3-x5;
run;
quit;