SAS中文论坛

标题: 新手求助！方程组求解的问题，内详进来看看！ [打印本页]

作者: shiyiming 时间: 2008-7-25 08:50
标题: 新手求助！方程组求解的问题，内详进来看看！
我想把var数据集里面的观测值一一带入model过程中求出 A和Delta_A，然后利用这两个值求出每条观测值对应的DD和EDF，为什么输出的内容显示只对第一条观测值求解了，而第二条没呢？

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;

欢迎光临 SAS中文论坛 (https://mysas.net/forum/)