求助：如何解决SAS问题呢？在线等，谢谢！-经管之家官网！

请问这个SAS程序出了什么问题，该如何解决呢？我把log贴出来，在线等。谢谢！
2573proc iml;
NOTE: IML Ready
2574
2575use Stag3data2;
2576read all into data;
NOTE: I/O required temporary file to be opened.
2577/*show names; */
2578/*print data;*/
2579start maxFunc(F) global(data);
2580sum=0;
2581Uf=J(1,3,0);
2582Um=J(1,3,0);
2583Ud=J(1,3,0);
2584Pf=J(1,3,0);
2585Pm=J(1,3,0);
2586Pd=J(1,3,0);
2587Vh=J(1,2,0);
2588G=J(1,3,0);
2589do i=1 to 1988;
2590prob=0;
2591do j=1 to 4;
2592
2593/*stag1*/
2594
2595/*Pf1 stands for fπ1*/
2596
2597z=(i-1)*4+j;
2598
2599/* change to z */
2600
2601Uf[1]=F[1]*data[z,5]+F[2]*data[z,6]+F[3]*data[z,7]+F[4]*data[z,8]
2602+F[5]*data[z,9]+F[6]*data[z,10]+F[7]*data[z,11]+F[8]*data[z,12]+F[9]*dat
a[z,13]+F[10]*data[z,14];
2603
2604Uf[2]=F[11]*data[z,15]+F[12]*data[z,16]+F[13]*data[z,17]+F[14]*data[z,18
]
2605+F[15]*data[z,19]+F[16]*data[z,20]+F[17]*data[z,21]+F[18]*data[z,22]+F[1
9]*data[z,23]+F[20]*data[z,24];
2606
2607Uf[3]=F[21];
2608
2609Pf[1]=exp(Uf[1])/(exp(Uf[1])+exp(Uf[2])+exp(Uf[3]));
2609! /* Pf[1] stand for f
π1 */
2610Pf[2]=exp(Uf[2])/(exp(Uf[1])+exp(Uf[2])+exp(Uf[3]));
2611Pf[3]=exp(Uf[3])/(exp(Uf[1])+exp(Uf[2])+exp(Uf[3]));
2612
2613/*Pm1 stand for mπ1*/
2614
2615Um[1]=F[22]*data[z,5]+F[23]*data[z,6]+F[24]*data[z,7]+F[25]*data[z,8]
2616+F[26]*data[z,9]+F[27]*data[z,10]+F[28]*data[z,11]+F[29]*data[z,12]+F[30
]*data[z,13]+F[31]*data[z,14];
2617
2618Um[2]=F[32]*data[z,15]+F[33]*data[z,16]+F[34]*data[z,17]+F[35]*data[z,18
]
2619+F[36]*data[z,19]+F[37]*data[z,20]+F[38]*data[z,21]+F[39]*data[z,22]+F[4
0]*data[z,23]+F[41]*data[z,24];
2620
2621Um[3]=F[42];
2622Pm[1]=exp(Um[1])/(exp(Um[1])+exp(Um[2])+exp(Um[3]));
2623Pm[2]=exp(Um[2])/(exp(Um[1])+exp(Um[2])+exp(Um[3]));
2624Pm[3]=exp(Um[3])/(exp(Um[1])+exp(Um[2])+exp(Um[3]));
2625
2626/*Pd1 stand for dπ1*/
2627
2628Ud[1]=F[43]*data[z,5]+F[44]*data[z,6]+F[45]*data[z,7]+F[46]*data[z,8]
2629+F[47]*data[z,9]+F[48]*data[z,10]+F[49]*data[z,11]+F[50]*data[z,12]+F[51
]*data[z,13]+F[52]*data[z,14];
2630
2631Ud[2]=F[53]*data[z,15]+F[54]*data[z,16]+F[55]*data[z,17]+F[56]*data[z,18
]
2632+F[57]*data[z,19]+F[58]*data[z,20]+F[59]*data[z,21]+F[60]*data[z,22]+F[6
1]*data[z,23]+F[62]*data[z,24];
2633
2634Ud[3]=F[63];
2635Pd[1]=exp(Ud[1])/(exp(Ud[1])+exp(Ud[2])+exp(Ud[3]));
2636Pd[2]=exp(Ud[2])/(exp(Ud[1])+exp(Ud[2])+exp(Ud[3]));
2637Pd[3]=exp(Ud[3])/(exp(Ud[1])+exp(Ud[2])+exp(Ud[3]));
2638
2639/*stag2*/
2640
2641Vf=
2642Vh[1]=F[64]*data[z,28]+F[65]*data[z,31]+F[66]*data[z,34];
2642! /*Vh stand for
house vote */
2643Vh[2]=F[64]*data[z,29]+F[65]*data[z,32]+F[66]*data[z,35];
2644
2645G[1]=(1-Pf[3]*Pm[3]*Pd[3])*(exp(Vh[1])/(exp(Vh[1])+exp(Vh[2])));
2645! /*G sta
nd for probability of household*/
2646G[2]=(1-Pf[3]*Pm[3]*Pd[3])*(exp(Vh[2])/(exp(Vh[1])+exp(Vh[2])));
2647G[3]=Pf[3]*Pm[3]*Pd[3];
2648
2649/*stag3*/
2650
2651Prof=Pf[1]*data[z,28]+Pf[2]*data[z,29]+Pf[3]*data[z,30];
2651! /*Pf means to c
hoose the available fπ1,2,3*/
2652Prom=Pm[1]*data[z,31]+Pm[2]*data[z,32]+Pm[3]*data[z,33];
2653Prod=Pd[1]*data[z,34]+Pd[2]*data[z,35]+Pd[3]*data[z,36];
2654
2655G[3]=G[1]*data[z,25]+G[2]*data[z,26];
2655! /* G[3] stand f
or group choice*/
2656fP11=1;
2657fP21=exp(Uf[1])/(exp(Uf[1])+exp(Uf[3]));
2657! /*fP21 :probabili
ty of switch from channel 2 to 1*/
2658fP23=exp(Uf[3])/(exp(Uf[1])+exp(Uf[3]));
2659fP12=exp(Uf[2])/(exp(Uf[2])+exp(Uf[3]));
2660fP22=1;
2661fP13=exp(Uf[3])/(exp(Uf[2])+exp(Uf[3]));
2662
2663mP11=1;
2664mP21=exp(Um[1])/(exp(Um[1])+exp(Um[3]));
2665mP23=exp(Um[3])/(exp(Um[1])+exp(Um[3]));
2666mP12=exp(Um[2])/(exp(Um[2])+exp(Um[3]));
2667mP22=1;
2668mP13=exp(Um[3])/(exp(Um[2])+exp(Um[3]));
2669
2670dP11=1;
2671dP21=exp(Ud[1])/(exp(Ud[1])+exp(Ud[3]));
2672dP23=exp(Ud[3])/(exp(Ud[1])+exp(Ud[3]));
2673dP12=exp(Ud[2])/(exp(Ud[2])+exp(Ud[3]));
2674dP22=1;
2675dP13=exp(Ud[3])/(exp(Ud[2])+exp(Ud[3]));
2676
2677fP=fP11*data[z,41]+fP21*data[z,42]+fP23*data[z,43]+fP12*data[z,44]+fP22*
data[z,45]+fP13*data[z,46];
2678mP=mP11*data[z,47]+mP21*data[z,48]+mP23*data[z,49]+mP12*data[z,50]+mP22*
data[z,51]+mP13*data[z,52];
2679dP=dP11*data[z,53]+dP21*data[z,54]+dP23*data[z,55]+dP12*data[z,56]+dP22*
data[z,57]+dP13*data[z,58];
2680
2681prob=prob + Prof*Prom*Prod*G[3]*fP*mP*dP;
2682end;
2683sum=sum*prob ;
2684end;
2685return(sum);
2686finish maxFunc;
NOTE: Module MAXFUNC defined.
2687
2688optn={1 2};
2689X=J(1,66,0);
2690call nlpcg(rc,xres,"maxFunc",X,optn);
ERROR: (execution) Invalid subscript or subscript out of range.
operation : [ at line 2601 column 16
operands: DATA, Z, *LIT1035
DATA 7948 rows 58 cols (numeric)
Z 1 row 1 col (numeric)
7949
*LIT1035 1 row 1 col (numeric)
5
statement : ASSIGN at line 2601 column 1
traceback : module MAXFUNC at line 2601 column 1
ERROR: (execution) Invalid subscript or subscript out of range.
operation : NLPCG at line 2690 column 1
operands: *LIT1332, X, OPTN
*LIT1332 1 row 1 col (character, size 7)
maxFunc
X 1 row 66 cols (numeric)
OPTN 1 row 2 cols (numeric)
1 2
statement : CALL at line 2690 column 1
2691
2692run;
NOTE: Module MAIN is undefined in IML; cannot be RUN.
2693quit;
NOTE: Exiting IML.
NOTE: 299 workspace compresses.
NOTE: SAS 系统由于错误而停止了该步的处理。
NOTE: “PROCEDURE IML”所用时间（总处理时间）:
实际时间 2.46 秒
CPU 时间 1.93 秒