我想把var数据集里面的观测值一一带入model过程中求出 A和Delta_A,为什么输出的内容显示只解了一次?
data VAR;
input stkcd $ E L DP Delta_E R;
cards;
000408 366316464.1 789769262 715198830.4 0.43 0.0225
000498 785495875.1 245927678 678942212 0.41 0.0225
;
run;
PROC MODEL DATA=VAR;
FIT START=(A=1000000000 Delta_A=0.2)/OUTEST=ESTIMATE;
D1=(LOG(A/L)+(R+Delta_A*Delta_A/2))/Delta_A;
D2=D1-Delta_A;
EQ.F1=(E+L*EXP(-R)*PROBNORM(D2))/PROBNORM(D1)-A;
EQ.F2=Delta_E*E/(PROBNORM(D1)*Delta_A)-A;
SOLVE A Delta_A/SOLVEPRINT;
ID E L Delta_E R;
RUN;
DATA EDF(keep=stkcd A Delta_A DD EDF);
MERGE VAR ESTIMATE;
DD=(A-DP)/(A*Delta_A);
EDF=(1-PROBNORM(DD))*100;
RUN;
proc print data=EDF;run;
[此贴子已经被作者于2008-7-5 14:31:08编辑过]