sas如何用数据算出公示里的未知的数

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

现在有个公式是y=a(1-e^(-bx)) 已经有data关于y和x  我就想知道怎么用sas来算出公式里的a和b 希望大神能给点建议

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：如何用 Data 如何

谢谢楼主的分享啊~

我没分享什么阿

2个未知数，一个方程，可以解出来？你是说怎么用图表示出这个公式吗？matlab或者r可以表示公式。

I had done a project about Chemical Kinetics first order increase. The formula was very similar to yours.
I use the code like

1. dert.y = b *(a- y);
   
2.         fit y initial =(y = y0) / time =x dynamic;
   
3.         run;

复制代码

in MODEL proc. and initialize y0 as 0. Of course you can fit it directly. NLIN works as well.
Jingju

longgb246 发表于 2014-5-25 09:10
2个未知数，一个方程，可以解出来？你是说怎么用图表示出这个公式吗？matlab或者r可以表示公式。

是两个未知数 但我现在有X 和Y 的数据 想求出a 和b

jingju11 发表于 2014-5-25 11:18
I had done a project about Chemical Kinetics first order increase. The formula was very similar to y ...

i saw the proc nlin example but in i don't know how to set the initial parms

Here is a trick you can play with nlin. See  below,

1) Suppose your data is generated by
a=0.5123; b=0.72135;
y=a*(1-exp(-b*x));

2) a and b can be recovered with proc nlin.

data t1;
a=0.5123; b=0.72135;
x=0.3;
y=a*(1-exp(-b*x));
output;
x=0.731;
y=a*(1-exp(-b*x));
output;

run;

proc nlin data=t1;
parms a=1 b=1;
model y=a*(1-exp(-b*x));
run;

The regression error is zero as expected.


                                             Sum of        Mean               Approx
           Source                    DF     Squares      Square    F Value    Pr > F

           Model                      2      0.0540      0.0270        .       .
           Error                      0    1.11E-16           .
           Corrected Total            1     0.00608


                                                Approx
                  Parameter      Estimate    Std Error    Approximate 95% Confidence Limits

                  a                0.5123            .           .           .
                  b                0.7214            .           .           .

bobguy 发表于 2014-5-26 07:33
Here is a trick you can play with nlin. See  below,

1) Suppose your data is generated by

万分感谢