最近用TSP分析几个数据的时候有的时候出现warning,而有的时候却不出现,方法是三段阶最小二乘法。
程序命令是:
OPTIONS MEMORY=20;
FREQ N;
SMPL 1 5468;
PANEL(ID=NO,MEAN) L C PR PL DIR DP IP BL PPR;
GENR LL=LOG(L);
GENR LPL=LOG(PL);
GENR LDP=LOG(DP);
GENR LBL=LOG(BL);
GENR DDDP=PR-PPR>0;
GENR DDSP=PR-PPR<=0;
GENR DDP=(PR-PPR)*DDDP;
GENR DSP=(PR-PPR)*DDSP;
? PRINT DDDP DDP DDSP DSP;
PARAM A0 A1 A2 A3 B0 B1 B2 B3 B4 B5 M1 M2;
FRML EQD LL=A0+A1*PR+A2*IP+A3*LPL-(M1/(1-M1))*(B1-A1)*DDP;
FRML EQS LL=B0+B1*(PR-DIR)+B2*LDP+B3*LBL+B4*LPL+B5*IP-(M2/(1-M2))*(B1-A1)*DSP;
LSQ (INST=(C IP LPL DIR LDP LL DDP DSP)) EQD;
LSQ (INST=(C IP LPL DIR LDP LL DDP DSP)) EQS;
LSQ (INST=(C IP LPL DIR LDP LL DDP DSP)
,MAXITW=0,WNAME=OWN) EQD EQS;
STOP; END;
noprint;
load NO YEAR L PR PL DIR DP IP BL PPR;
结果里面出现错误的地方是(一部分):
Standard
Parameter Estimate Error t-statistic P-value
A0 -23.0698 8.26965 -2.78969 [.005]
A1 -3.33639 .221767 -15.0446 [.000]
A2 -.425091E-02 .012209 -.348174 [.728]
A3 3.37706 .593448 5.69058 [.000]
M1 -.177639E+10 0. 0. [1.00]
B1 .036742 .061883 .593728 [.553]
B0 2.09118 .855461 2.44450 [.015]
B2 .015868 .011979 1.32460 [.185]
B3 .996501 .858339E-02 116.096 [.000]
B4 -.181335 .077785 -2.33123 [.020]
B5 .312239E-03 .181204E-02 .172313 [.863]
M2 .082014 .032605 2.51537 [.012]
*** WARNING in command 18 Procedure LSQ: At least one coefficient in the
table above could not be estimated due to singularity of the data or
derivatives.
先谢谢各位大仙了。


雷达卡




京公网安备 11010802022788号







