求助：一个非常奇怪的sas程序问题？-经管之家官网！

程序没有报错，当是不能运行出结果，有提示“NOTE: Module MAIN is undefined in IML; cannot be RUN.”而我在之前的用run来跑程序时事没有问题的，现在非常confused，不知道如何解决好？大家有什么好想法在此分享，小弟不胜感激。在网上查了原因是"If a RUN statement is submitted without specifying the name of the module that you wish to execute, the RUN statement attempts to run the module named MAIN. If you have defined a module without a name, or with a name of MAIN, then RUN; is a valid statement that will run the module whose name defaulted to MAIN. "
2429proc iml;
NOTE: IML Ready
2430
2431use Stag3data2;
2432read all into data;
NOTE: I/O required temporary file to be opened.
2433/*show names; */
2434/*print data;*/
2435start maxFunc(F) global(data);
2436sum=0;
2437Uf=J(1,3,0);
2438Um=J(1,3,0);
2439Ud=J(1,3,0);
2440Pf=J(1,3,0);
2441Pm=J(1,3,0);
2442Pd=J(1,3,0);
2443Vh=J(1,2,0);
2444G=J(1,3,0);
2445/*do i=1 to 1988;*/
2446do i=1 to 994;
2447prob=0;
2448do j=1 to 4;
2449
2450/*stag1*/
2451
2452/*Pf1 stands for fπ1*/
2453
2454z=(i-1)*4+j;
2455
2456/* change to z */
2457
2458Uf[1]=F[1]*data[z,5]+F[2]*data[z,6]+F[3]*data[z,7]+F[4]*data[z,8]
2459+F[5]*data[z,9]+F[6]*data[z,10]+F[7]*data[z,11]+F[8]*data[z,12]+F[9]*data[z,13]+F[10]*data[z,14];
2460
2461Uf[2]=F[11]*data[z,15]+F[12]*data[z,16]+F[13]*data[z,17]+F[14]*data[z,18]
2462+F[15]*data[z,19]+F[16]*data[z,20]+F[17]*data[z,21]+F[18]*data[z,22]+F[19]*data[z,23]+F[20]*data[z,24];
2463
2464Uf[3]=F[21];
2465
2466Pf[1]=exp(Uf[1])/(exp(Uf[1])+exp(Uf[2])+exp(Uf[3]));
2466! /* Pf[1] stand for fπ1 */
2467Pf[2]=exp(Uf[2])/(exp(Uf[1])+exp(Uf[2])+exp(Uf[3]));
2468Pf[3]=exp(Uf[3])/(exp(Uf[1])+exp(Uf[2])+exp(Uf[3]));
2469
2470/*Pm1 stand for mπ1*/
2471
2472Um[1]=F[22]*data[z,5]+F[23]*data[z,6]+F[24]*data[z,7]+F[25]*data[z,8]
2473+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];
2474
2475Um[2]=F[32]*data[z,15]+F[33]*data[z,16]+F[34]*data[z,17]+F[35]*data[z,18]
2476+F[36]*data[z,19]+F[37]*data[z,20]+F[38]*data[z,21]+F[39]*data[z,22]+F[40]*data[z,23]+F[41]*data[z,24];
2477
2478Um[3]=F[42];
2479Pm[1]=exp(Um[1])/(exp(Um[1])+exp(Um[2])+exp(Um[3]));
2480Pm[2]=exp(Um[2])/(exp(Um[1])+exp(Um[2])+exp(Um[3]));
2481Pm[3]=exp(Um[3])/(exp(Um[1])+exp(Um[2])+exp(Um[3]));
2482
2483/*Pd1 stand for dπ1*/
2484
2485Ud[1]=F[43]*data[z,5]+F[44]*data[z,6]+F[45]*data[z,7]+F[46]*data[z,8]
2486+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];
2487
2488Ud[2]=F[53]*data[z,15]+F[54]*data[z,16]+F[55]*data[z,17]+F[56]*data[z,18]
2489+F[57]*data[z,19]+F[58]*data[z,20]+F[59]*data[z,21]+F[60]*data[z,22]+F[61]*data[z,23]+F[62]*data[z,24];
2490
2491Ud[3]=F[63];
2492Pd[1]=exp(Ud[1])/(exp(Ud[1])+exp(Ud[2])+exp(Ud[3]));
2493Pd[2]=exp(Ud[2])/(exp(Ud[1])+exp(Ud[2])+exp(Ud[3]));
2494Pd[3]=exp(Ud[3])/(exp(Ud[1])+exp(Ud[2])+exp(Ud[3]));
2495
2496/*stag2*/
2497
2498Vf=
2499Vh[1]=F[64]*data[z,28]+F[65]*data[z,31]+F[66]*data[z,34];
2499! /*Vh stand for house vote */
2500Vh[2]=F[64]*data[z,29]+F[65]*data[z,32]+F[66]*data[z,35];
2501
2502G[1]=(1-Pf[3]*Pm[3]*Pd[3])*(exp(Vh[1])/(exp(Vh[1])+exp(Vh[2])));
2502! /*G stand for probability of household*/
2503G[2]=(1-Pf[3]*Pm[3]*Pd[3])*(exp(Vh[2])/(exp(Vh[1])+exp(Vh[2])));
2504G[3]=Pf[3]*Pm[3]*Pd[3];
2505
2506/*stag3*/
2507
2508Prof=Pf[1]*data[z,28]+Pf[2]*data[z,29]+Pf[3]*data[z,30];
2508! /*Pf means to choose the available fπ1,2,3*/
2509Prom=Pm[1]*data[z,31]+Pm[2]*data[z,32]+Pm[3]*data[z,33];
2510Prod=Pd[1]*data[z,34]+Pd[2]*data[z,35]+Pd[3]*data[z,36];
2511
2512G[3]=G[1]*data[z,25]+G[2]*data[z,26];
2512! /* G[3] stand for group choice*/
2513fP11=1;
2514fP21=exp(Uf[1])/(exp(Uf[1])+exp(Uf[3]));
2514! /*fP21 :probability of switch from channel 2 to 1*/
2515fP23=exp(Uf[3])/(exp(Uf[1])+exp(Uf[3]));
2516fP12=exp(Uf[2])/(exp(Uf[2])+exp(Uf[3]));
2517fP22=1;
2518fP13=exp(Uf[3])/(exp(Uf[2])+exp(Uf[3]));
2519
2520mP11=1;
2521mP21=exp(Um[1])/(exp(Um[1])+exp(Um[3]));
2522mP23=exp(Um[3])/(exp(Um[1])+exp(Um[3]));
2523mP12=exp(Um[2])/(exp(Um[2])+exp(Um[3]));
2524mP22=1;
2525mP13=exp(Um[3])/(exp(Um[2])+exp(Um[3]));
2526
2527dP11=1;
2528dP21=exp(Ud[1])/(exp(Ud[1])+exp(Ud[3]));
2529dP23=exp(Ud[3])/(exp(Ud[1])+exp(Ud[3]));
2530dP12=exp(Ud[2])/(exp(Ud[2])+exp(Ud[3]));
2531dP22=1;
2532dP13=exp(Ud[3])/(exp(Ud[2])+exp(Ud[3]));
2533
2534fP=fP11*data[z,41]+fP21*data[z,42]+fP23*data[z,43]+fP12*data[z,44]+fP22*data[z,45]+fP13*data[z,46];
2535mP=mP11*data[z,47]+mP21*data[z,48]+mP23*data[z,49]+mP12*data[z,50]+mP22*data[z,51]+mP13*data[z,52];
2536dP=dP11*data[z,53]+dP21*data[z,54]+dP23*data[z,55]+dP12*data[z,56]+dP22*data[z,57]+dP13*data[z,58];
2537
2538prob=prob + Prof*Prom*Prod*G[3]*fP*mP*dP;
2539end;
2540sum=sum*prob ;
2541end;
2542return(sum);
2543finish maxFunc;
NOTE: Module MAXFUNC defined.
2544
2545optn={1 2};
2546X=J(1,66,0);
2547X[1,1] = 0.7865;
2548X[1,2] = -1.173888;
2549X[1,3] = -0.901912;
2550X[1,4] = 2.52522;
2551X[1,5] = 1.934948;
2552X[1,6] = 5.106454;
2553X[1,7] = -2.823555;
2554X[1,8] = 3.528326;
2555X[1,9] = -2.357012;
2556X[1,10] = -0.525734;
2557X[1,11] = -12.653302;
2558X[1,12] = -3.099612;
2559X[1,13] = 4.539191;
2560X[1,14] = 2.001213;
2561X[1,15] = 0.153087;
2562X[1,16] = 2.228049;
2563X[1,17] = -1.640786;
2564X[1,18] = -0.220159;
2565X[1,19] = -1.042585;
2566X[1,20] = 5.483434;
2567X[1,21] = -1.852272 ;
2568X[1,64]=3;
2569X[1,65]=7;
2570X[1,66]=10;
2571call nlpcg(rc,xres,"maxFunc",X,optn);
NOTE: ABSGCONV convergence criterion satisfied.
2572run;
NOTE: Module MAIN is undefined in IML; cannot be RUN.
2573quit;
NOTE: Exiting IML.
NOTE: 30078 workspace compresses.
NOTE: “PROCEDURE IML”所用时间（总处理时间）:
实际时间 3:02.79
CPU 时间 3:02.34