| 查看: 1411 | 回复: 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 ] |
回复此楼» 猜你喜欢
综述论文作为代表作会不会影响评审专家的印象分？
已经有12人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有4人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有29人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有9人回复
-
是这周出结果还是下周出结果？
已经有3人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有9人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有4人回复
-
帮忙看看fileCode
已经有13人回复
-
为什么网上很多人说本周 12号出结果
已经有7人回复
-
读博的好处
已经有3人回复
找到一些相关的精华帖子，希望有用哦~
求助一个函数的解法编程的源程序（有效期5天）
已经有3人回复
- <求助>学习SPSS SAS 还是R语言。
已经有21人回复
- 求助 matlab编程,对于你来说可能很简单
已经有4人回复
- 求助SAS 多元非线性回归拟合编程
已经有3人回复
- 【求助】同样一个问题用matlab和fortran编程计算，结果相差比较大，不知道为什么
已经有4人回复
- 【求助】有关“SW0PC_FXGP编程软件”的介绍
已经有1人回复
- 【求助】SAS编程出了错，看不懂，请教高手
已经有8人回复
- 【求助】蒙卡编程问题
已经有15人回复
- 【求助】小町算编程【已完成】
已经有15人回复
- 【求助】求助：vb编程中用牛顿迭代解三次方程为什么只得到一个根？【已完成】
已经有15人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
