SAS 多个NLIN过程放一张图中，做PLOT图求教

各位大神，遇到一个很大的难题，有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;

crazebird 发表于 2019-3-9 22:42
你好，splot过程怎么写，你能帮改写一下吗，SAS我基本是小白，这里面的程序我都是照葫芦画瓢的，万分感谢 ...

帮你改好了，把数据填充完整就可以得到你给出的结果了

1. data expd;
   
2.         input x y z o @@;
   
3.         datalines;
   
4.         23  4497  4331  3206
   
5.         23  4472  4283  3143
   
6.         23  4158  4004  2962
   
7.         25  4241  4083  3031
   
8.         25  3985  3845  2866
   
9.         26  4240  4092  3062
   
10.         26  4915  4742  3530
    
11.         27  3271  3147  2339
    
12.         27  3133  3013  2234
    
13.         27  4850  4681  3488
    
14.         27  4964  4796  3588
    
15.         28  3760  3624  2685
    
16.         28  5055  4883  3695
    
17.         28  4522  4354  3270
    
18.         29  5458  5262  3907
    
19.         29  4575  4423  3254
    
20.         30  4933  4749  3574
    
21.         30  5464  5286  3968
    
22. ;
    
23. 
    
24. proc nlin data = expd best = 10 method = gauss;
    
25.         parms a=6000 to 12000 by 100 b=-3 to -5 by -1 ;
    
26.         model y=a*(1-exp(-exp(b)*(x)));
    
27.         output out = expout1 p = ygs;
    
28. run;
    
29. 
    
30. proc nlin data = expd best = 10 method = gauss;
    
31.         parms c=6000 to 12000 by 100 d=-3 to -5 by -1 ;
    
32.         model z=c*(1-exp(-exp(d)*(x)));
    
33.         output out = expout2 p = ygsa;
    
34. run;
    
35. 
    
36. proc nlin data = expd best = 10 method = gauss;
    
37.         parms e=6000 to 12000 by 100 f=-3 to -5 by -1 ;
    
38.         model o=e*(1-exp(-exp(f)*(x)));
    
39.         output out = expout3 p = ygsb;
    
40. run;
    
41. 
    
42. data expout;
    
43.         merge expout1-expout3;
    
44. run;
    
45. 
    
46. goptions reset = global gunit = pct cback = white noborder htitle = 3 htext = 3 colors = (black);
    
47. proc gplot data = expout;
    
48.         plot y*x ygs*x z*x o*x ygsa*x ygsb*x /haxis=axis1 vaxis=axis2 overlay noframe legend hminor=0;
    
49. 
    
50.         symbol1 i=none v=dot cv=red h=2.5 w=2;
    
51.         symbol2 i=spline v=none l=1 h=2.5 w=2.5;
    
52.         symbol3 i=none v=star cv=blue h=2.5 w=2;
    
53.         symbol4 i=none v=plus cv=green h=2.5 w=2;
    
54.         symbol5 i=spline v=none l=3 h=2.5 w=2.5;
    
55.         symbol6 i=spline v=none l=2 h=2.5 w=2.5;
    
56. 
    
57.         axis1 minor=(number=1) order=(20 to 120 by 20) label=(h=13pt f="calibri" "Body Weight (kg)") ;
    
58.         axis2 minor=(number=1) order=(0 to 12000 by 2000) label=(h=13pt a=90 r=0 f=calibri "NE Intake (kcal/day)");
    
59. run;
    
60. quit;

复制代码

这里能写多个plot语句的话就可以，否则建议使用sgplot过程

whymath 发表于 2019-3-9 18:44
这里能写多个plot语句的话就可以，否则建议使用sgplot过程

你好，splot过程怎么写，你能帮改写一下吗，SAS我基本是小白，这里面的程序我都是照葫芦画瓢的，万分感谢了

这个如果用其他软件来解决的话，so easy

nkunku 发表于 2019-3-10 08:11
这个如果用其他软件来解决的话，so easy

你好，用什么软件啊，能分享一下吗，谢谢

crazebird 发表于 2019-3-10 09:33
你好，用什么软件啊，能分享一下吗，谢谢

Origin等

nkunku 发表于 2019-3-10 14:19
Origin等

你好，不会用阿

crazebird 发表于 2019-3-8 23:11
各位大神，遇到一个很大的难题，有3组数据，想分别用NLIN模型运算，放在一张图中，但是程序错误不可用。 ...

自己来评论一下，各位大神们，有会SAS编程的吗，帮我编一下啊，真的真的是万分感谢了

whymath 发表于 2019-3-8 23:11
帮你改好了，把数据填充完整就可以得到你给出的结果了

大神，真的是太感谢你了，太厉害了，现在还有一个问题，图例 legend，能否让其只显示想要展示的图例，比如一共6组数据，只想显示其中的4个或者5个图例，原代码和图如下所示

proc gplot data = expout;
plot y*x ygs*x z*x o*x ygsa*x ygsb*x /haxis=axis1 vaxis=axis2 overlay noframe legend=legend1 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;
legend1 down=3 position=(top right inside) frame
label=(f=calibri h=1.5 "") value=(f=calibri "" "LO" "DO" "WO" "DO" "EE");
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;
quit;