| 查看: 1408 | 回复: 0 | ||
arc2360金虫 (小有名气)
|
[求助]
求助检查一个简单SAS编程（将两条曲线放一个图形中）
|
|
求大师帮忙看下面程序，前面是两条非线性曲线的拟合，后面想把两条曲线放一个图中，但程序运行后没有显示，求助求助，谢谢！！ Data L1; input t y @@; cards; 0 3.804 24 3.877 48 3.930 72 5.01 96 5.959 120 5.813 144 6.55 168 7.094 192 7.587 216 6.8 240 6.82 264 7.084 288 7.065 ; Proc nlin method=Marquardt; Parms Y0=3 to 4 by 0.05 Ymax=6.5 to 7.5 by 0.1 u=0.01 to 0.05 by 0.005 M=100 to 150 by 5; Temp=exp(-u*(t-M)); model y=Y0+(Ymax-Y0) * exp(-exp(-u * (t-M))); der.Y0=1-exp(-temp); der.Ymax=exp(-temp); der.u=(Ymax-Y0) * exp(-temp) * (t-M) * temp; der.M=(Ymax-Y0) * exp(-temp) * u * (-temp); output out=expout p=py1; run; Data L1; input t y @@; cards; 0 3.804 24 3.877 48 3.930 72 5.01 96 5.959 120 5.813 144 6.55 168 7.094 192 7.587 216 6.8 240 6.82 264 7.084 288 7.065 ; Proc nlin method= Marquardt; Parms Y0=3 to 4 by 0.05 Ymax=6.5 to 7.5 by 0.1 u=0.01 to 0.05 by 0.005 M=55 to 85 by 5; Temp=exp(-u*(t-M)); model y=Y0+(Ymax-Y0)/(1+exp(-u * (t-M))); der.Y0=1-1/(1+temp); der.Ymax=1/(1+temp); der.u=(Ymax-Y0) * temp * (t-M)/((1+temp)*(1+temp)); der.M=-(Ymax-Y0) * u * temp/((1+temp)*(1+temp)); output out=expout p=py2; run; goptions reset=global gunit=pct cback=white border htitle=6 htext=3 ftext=swissb colors=(back); proc gplot data=expout; plot y*t py1*t py2*t /haxis=axis1 vaxis=axis2 overlay; symbol1 i=none v=plus cv=red h=2.5 w=2; symbol2 i=join v=none l=1 h=2.5 w=2; axis1 order=0 to 300 by 20; axis2 order=3 to 8 by 0.5; title1 'y=Y0+(Ymax-Y0)*exp(-exp(-u*(t-M)))' ; title2 ' y=Y0+(Ymax-Y0)/(1+exp(-u * (t-M)))' ; title3 'proc nlin method= Levenberg-Marquardt '; run; [ Last edited by arc2360 on 2012-3-21 at 21:51 ] |
回复此楼» 猜你喜欢
听说今天filecode变了
已经有34人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
基金中了
已经有6人回复
-
filecode
已经有14人回复
-
固定端突然变了，今天
已经有7人回复
-
filecode
已经有11人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有6人回复
-
一个有机合成实验室都需要哪些设备？
已经有12人回复
-
求各位大神看下
已经有33人回复
-
咨询面上基金
已经有7人回复
找到一些相关的精华帖子，希望有用哦~
求助一个函数的解法编程的源程序（有效期5天）
已经有3人回复
- <求助>学习SPSS SAS 还是R语言。
已经有21人回复
- 求助 matlab编程,对于你来说可能很简单
已经有4人回复
- 求助SAS 多元非线性回归拟合编程
已经有3人回复
- 【求助】同样一个问题用matlab和fortran编程计算，结果相差比较大，不知道为什么
已经有4人回复
- 【求助】有关“SW0PC_FXGP编程软件”的介绍
已经有1人回复
- 【求助】SAS编程出了错，看不懂，请教高手
已经有8人回复
- 【求助】蒙卡编程问题
已经有15人回复
- 【求助】小町算编程【已完成】
已经有15人回复
- 【求助】求助：vb编程中用牛顿迭代解三次方程为什么只得到一个根？【已完成】
已经有15人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
50