SAS中,logistic方程的参数已经拟合出来了,如何在散点图中添加趋势线或者是拟合曲线。想要得到图1的效果,但是我用以下程序做的图如图2,和图一完全不一样,请大神们帮忙看看哪里错了。急求用SAS绘制如图1的编程,谢谢。
data tt;
input t y;
cards;
0 0
3 1.09
4 8.19
5 94.09
6 211.41
7 328.76
8 492.9
9 306.48
11 231.13
15 155.09
;
proc nlin data=tt best=10 method=marquardt;
parms a=1000 to 6000 by 50 kg=1 to 1.5 by 0.01 tcg=5 to 7 by 0.1
kd=3 to 5 by 1 ;
model y=(a/(1+exp(kg*(tcg-t)))-a/(1+exp(kg*tcg)))*(exp(-t/kd));
run;quit;
data yc;
set tt;
yp=(2087.8/(1+exp(1.4983*(6.5275-t)))-2087.8/(1+exp(1.4983*6.5275)))*(exp(-t/5.2313));
proc gplot data=yc;
symbol1 i=spline v=dot;
symbol2 i=join v=none co=red;
plot y*t=1 yp*t=2 /overlay;
run;

图1
图2