最优子集回归程序里面都要改哪些变量啊

469482462

请问最优子集回归程序里面都要改哪些变量啊。我的因子有15个，一共54年数据，其中4年用来预测，我应该怎么修改去用程序呢？n,nn,ip各应该是多少？
c this is a program for calculating all possible subset regressions
c and selecting a optimal regression by couple score criterion
  program osr
parameter (n=50,nn=4,n1=n+nn)
        parameter (ip=14,iip=ip+1)                  
c      **********************************************************
c      * n:  sample size;                                       *
c      * nn: forecasting step number;                           *
c      * ip: predictor number;                                  *
c      * xy: array of predictors and predictand                 *
c      **********************************************************
  dimension x1(n1),x(n1)
common/aa/ex(30,30)
common/bb/coe(30,30)
common/cc/xy(100,30)
common/dd/a(30,30)
open(2,file=' ' ,status='old')
read(2,*) ((xy(i,j),i=1,n1),j=1,ip)
read(2,*) (xy(i,iip),i=1,n)
close(2)
open(4,file= ' ' ,status='new')
write(4,1)
    1 format(20x,'******OPTIMAL SUBSET REGRESSION******'/)
call sreg(x,x1,n,n1,ip,iip)
w=0.0
do 2 i=1,n
    2 w=w+(x(i)-x1(i))*(x(i)-x1(i))
v=sqrt(w/n)
write(4,3)v
    3 format(/20x,'THE FINAL RESULT OF FITTING AND FORECAST'/
     & (5x,'RMSE=',F15.2))
write(4,4)
    4 format(6x,'NO.',6x,'REG.',6X,'OBS.')
do 5 i=1,n1
write(4,6)i,x1(i),x(i)
    6 format(5x,i5,f10.2,f10.2)
    5 continue
stop
end
C------------------------------------END--------------------------------------

数据如下：
   145.141   129.149   139.691   135.841   138.256   144.270   153.183   138.790   126.827   144.205
   142.198   136.652   133.028   132.915   136.238   156.167   124.063   163.581   159.460   157.769
   119.615   123.198   172.549   125.433   145.005   156.127   173.551   158.686   138.491   136.976
   138.172   164.701   159.409   171.185   147.891   150.412   157.834   129.385   143.262   144.372
   118.150   150.587   137.935   146.142   144.729   151.477   158.385   149.471   149.151   146.577
   132.571   153.984   151.319   155.348
   118.094   107.336   114.885   118.798   117.939   115.742   107.373   125.505   114.759   125.052
   115.352   117.968   125.824   122.708   120.868   119.296   121.147   129.299   130.997   130.508
   119.785   134.947   129.207   128.807   114.050   124.219   132.980   134.471   128.729   130.545
   128.730   128.154   134.123   128.575   134.968   134.952   126.003   129.385   123.466   124.858
   114.907   132.991   125.532   129.362   122.047   120.759   123.972   114.292   119.647   119.129
   101.129   116.727   111.509   118.303
   143.471   172.836   126.557    94.742   155.469   148.404   198.879   130.765   132.889   146.304
   111.809   123.903    95.433   103.282   163.277    96.396   180.058   116.041    99.991    91.359
   171.568    79.615    97.516   108.760   113.353   176.919   125.925   115.968    85.351   128.506
   104.848    99.140    86.181   151.777   134.817   104.188   128.537   106.470   111.349   146.102
   177.490   145.237   139.789    57.748   167.401   202.635    86.238   114.302   114.593   183.517
   171.403   144.350   172.773   137.064
   123.911   108.881   120.305   115.212   124.924   127.808   125.153   117.197   106.823   125.478
   119.018   124.408   123.262   112.798   125.051   128.536   123.164   133.889   131.954   124.030
   135.518   133.835   142.851   119.878   128.634   129.113   141.152   124.403   116.140   126.236
   127.219   147.249   133.641   122.328   136.414   123.304   128.465   128.137   113.085   118.602
   109.634   128.441   136.613   125.206   122.685   121.164   136.815   117.448   132.176   131.258
   116.676   117.442   125.698   134.751
   187.532   179.260   110.413   171.855   151.179   147.334   185.264   196.508   135.986   133.731
   164.898   155.859   197.550   203.223   196.797   177.995   150.936   180.247   119.702   151.415
   209.930   175.767   237.866   212.911   142.384   211.770   155.124   186.183   216.592   210.713
   195.675   183.740   220.632   193.905   222.869   110.112   143.661   170.177   176.089   184.211
   184.234   201.868   197.426   177.713   207.434   168.423   208.379   204.617   174.418   153.695
   152.401   208.099   202.347   169.077
   178.581   154.109   165.315   160.532   147.686   152.903   162.047   138.196   149.025   149.965
   169.545   155.103   151.449   166.409   164.836   162.762   169.878   144.512   166.008   137.947
   149.799   150.691   147.752   160.186   161.137   147.643   142.878   150.608   168.095   162.545
   147.965   135.755   137.545   150.080   136.768   146.663   167.928   163.920   181.676   172.085
   160.506   168.975   135.178   163.928   142.569   153.636   170.995   172.706   160.876   148.736
   170.135   169.300   170.612   149.764
    24.974    25.051    25.010    25.466    24.921    24.654    24.843    25.201    25.133    24.745
    24.979    24.789    24.501    25.084    25.009    24.836    24.675    24.609    24.557    24.702
    24.835    24.773    24.537    24.976    24.735    24.741    24.409    24.827    25.096    24.791
    25.108    24.798    24.360    24.416    24.983    25.090    24.431    24.777    25.419    25.525
    25.379    25.455    25.039    25.121    25.077    25.454    25.018    25.437    25.368    25.131
    25.429    25.480    25.044    25.211
     9.202     9.750     9.474     9.191     9.493     9.735     9.611     9.813     9.173     9.520
     9.938     9.249     9.779     9.293     9.694     9.049     9.211     9.033     9.990     9.735
     9.305     9.889     9.650     9.715     9.658    10.364    10.754     9.755     9.446    10.379
    10.024    10.651    10.005    10.857     9.889     9.945    10.150    10.181     9.976     9.288
     9.968     9.472    10.464    10.262    10.208    10.040    10.193     9.599     9.662     9.535
     9.321     9.620     9.451    10.884
    23.477    23.631    23.592    23.712    23.132    23.294    23.346    23.356    23.324    23.264
    23.472    23.452    23.261    23.505    23.643    23.436    23.378    23.198    23.284    23.140
    23.441    23.369    23.157    23.287    23.246    23.229    23.349    23.592    23.669    23.271
    23.491    23.289    23.023    23.241    23.741    23.447    23.217    23.686    24.375    24.118
    24.138    23.991    23.762    23.812    23.756    23.626    23.869    24.106    23.846    23.801
    24.067    24.061    23.692    23.550
    26.980    26.424    26.258    26.770    26.323    26.407    27.140    26.687    27.432    26.995
    26.282    25.997    26.193    26.158    26.224    26.092    26.146    26.800    26.586    26.923
    26.800    26.821    26.487    26.752    26.278    26.771    26.723    26.881    26.649    26.803
    27.268    26.614    26.814    26.788    26.477    27.180    26.725    26.460    26.113    25.976
    26.282    26.497    26.597    26.945    26.742    26.691    26.740    26.437    26.372    27.032
    26.417    26.218    26.573    26.798
     4.088     4.673     4.586     4.978     4.662     4.388     4.897     4.536     4.820     4.527
     5.126     4.924     4.353     4.968     4.509     4.597     4.268     4.605     4.510     4.712
     4.475     4.823     4.727     4.884     5.036     4.858     4.929     4.380     4.960     4.852
     5.810     4.873     4.472     4.929     5.413     5.720     5.629     5.080     4.740     4.408
     4.439     4.690     4.710     4.802     4.600     5.553     5.265     4.468     4.962     4.670
     5.305     5.012     5.331     5.601
    24.206    24.073    23.954    24.093    23.371    23.713    24.037    23.798    23.668    23.858
    23.916    23.502    23.569    23.773    24.063    23.830    23.672    23.451    23.672    23.154
    23.666    23.774    23.187    23.740    23.612    23.787    23.829    24.080    24.546    23.935
    23.734    23.697    23.288    23.475    23.969    23.841    23.505    23.831    24.448    24.411
    24.561    24.253    24.113    24.232    23.778    24.384    24.188    24.571    24.031    24.004
    24.327    24.352    24.150    23.887
    22.304    22.235    21.785    22.190    21.847    22.014    22.309    21.778    22.169    21.918
    21.990    21.881    21.947    22.093    22.037    21.872    21.857    21.743    22.025    21.402
    21.664    21.735    21.510    21.926    21.663    21.744    21.738    22.267    22.373    22.169
    21.688    21.875    21.666    21.913    21.922    21.771    21.993    22.006    22.858    22.689
    22.739    22.554    22.111    22.359    22.119    22.327    22.337    22.666    22.252    22.039
    22.297    22.475    22.349    22.082
    26.037    25.912    26.521    26.142    25.942    26.473    26.046    26.619    27.044    26.404
    25.545    26.065    26.338    25.657    25.634    25.477    26.033    26.373    26.041    26.979
    26.317    26.606    26.319    26.222    26.096    26.704    26.662    26.653    25.549    26.693
    26.974    26.581    26.981    26.994    26.753    26.533    26.537    26.486    25.512    25.604
    25.900    26.425    26.756    26.405    26.719    25.963    26.165    25.400    25.997    26.612
    25.795    25.717    25.952    26.259
    12.159    12.527    13.031    12.285    12.147    11.518    12.019    12.164    11.427    12.084
    11.724    11.678    12.107    11.706    11.750    12.023    12.685    12.195    12.039    12.406
    12.410    11.983    12.004    12.388    12.290    12.453    12.224    11.734    11.776    12.026
    12.550    12.555    12.365    12.543    12.228    12.817    12.789    12.154    11.547    11.668
    11.675    11.512    12.402    12.065    12.795    12.231    11.752    11.559    12.106    12.042
    11.757    12.086    11.614    13.163

飞飞

求大神们指点，我也很想知道

Clover

楼主，请问这个问题解决了吗？我最近也在用这个程序，也有相同的疑问。

迷糊加油

关键部分call sreg(x,x1,n,n1,ip,iip)没有贴出来。。。{:eb303:}