各位大神，遇到一个很大的难题，有3组数据，想分别用NLIN模型运算，放在一张图中，但是程序错误不可用。研究了好多天，图片只能显示3组数据的散点图和其中1组数据的NLIN拟合图，如图所示，之前的程序也在下面。我最终的结果是想把3组数据的散点图，以及3组数据的NLIN模型拟合图均放在一张图片中，请大神们赐教，拜托了，谢谢

data expd;
input x y z o @@;
datalines;

23  4497  4331  3206
23  4472  4283  3143
23  4158  4004  2962
25  4241  4083  3031
25  3985  3845  2866
26  4240  4092  3062
26  4915  4742  3530
27  3271  3147  2339
27  3133  3013  2234
27  4850  4681  3488
27  4964  4796  3588
28  3760  3624  2685
28  5055  4883  3695
28  4522  4354  3270
29  5458  5262  3907
29  4575  4423  3254
30  4933  4749  3574
30  5464  5286  3968

;

proc nlin data = expd best = 10 method = gauss;
parms a=6000 to 12000 by 100 b=-3 to -5 by 1 ;
model y=a*(1-exp(-exp(b)*(x)));
output out = expout p = ygs;
run;
proc nlin data = expd best = 10 method = gauss;
parms c=6000 to 12000 by 100 d=-3 to -5 by 1 ;
model z=c*(1-exp(-exp(d)*(x)));
output out = expout p = ygsa;
run;
proc nlin data = expd best = 10 method = gauss;
parms e=6000 to 12000 by 100 f=-3 to -5 by 1 ;
model o=e*(1-exp(-exp(f)*(x)));
output out = expout p = ygsb;
run;
goptions reset = global gunit = pct cback = white noborder
         htitle = 3 htext = 3 ftext = calibri colors = (black);
proc gplot data = expout;
plot y*x ygs*x z*x o*x ygsa*x ygsb*x /haxis=axis1 vaxis=axis2  overlay noframe legend hminor=0;
symbol1 i=none v=dot cv=red h=2.5 w=2;
symbol2 i=spline v=none l=1 h=2.5 w=2.5;
symbol3 i=none v=star cv=blue h=2.5 w=2;
symbol4 i=none v=plus cv=green h=2.5 w=2;
symbol5 i=spline v=none l=3 h=2.5 w=2.5;
symbol6 i=spline v=none l=2 h=2.5 w=2.5;
axis1 minor=(number=1)  order=(20 to 120 by 20)
label=(h=13pt f="calibri" "Body Weight (kg)") ;
axis2 minor=(number=1) order=(0 to 12000 by 2000)
label=(h=13pt a=90 r=0 f=calibri "NE Intake (kcal/day)");
run;