楼主: ohmypig
1817 1

[编程问题求助] How to solve a system of nonlinear equations with cross-equation correlations? [推广有奖]

  • 0关注
  • 0粉丝

本科生

79%

还不是VIP/贵宾

-

威望
0
论坛币
288 个
通用积分
28.5459
学术水平
3 点
热心指数
4 点
信用等级
2 点
经验
829 点
帖子
79
精华
0
在线时间
106 小时
注册时间
2010-7-15
最后登录
2018-1-31

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
Hi all,
I try to solve a system of nonlinear equations with STATA.
My data looks like this:
id        s                   x              e
1  .0064259    -1.008586        ?
1  .2882765        2.701381        ?
1  .1368153        2.010069        ?
where id=1,2,3,...J,


Formula: s_j=exp(x_j+e_j)/A, and A=(1+exp(x_1+e_1)+...+exp(x_j+e_j)+...+exp(x_J+e_J))
I have exactly J equations with J unknowns. I want to know e_js.

I found http://www.stata.com/support/faqs/programming/system-of-nonlinear-equations/ and http://www.stata.com/statalist/archive/2011-06/msg00674.html
but I don't know how to add "A".
I also found "nlsur", which allows me to add "A". Here is what I did.
[1] I reformatted the data as the following:
id             s1              x1                    s2                     x2              s3           x3                     s11            x11             s22                  x22               s33           x33
1        .0064259        -1.008586        .2882765        2.701381        .1368153        2.010069        .0064259        -1.008586                               
1        .0064259        -1.008586        .2882765        2.701381        .1368153        2.010069                                         .2882765        2.701381               
1        .0064259        -1.008586        .2882765        2.701381        .1368153        2.010069                                                                                .1368153        2.010069       

...

code is as the following:
nlsur (s11=exp(x11+{e1})/(1+exp(x1+{e1})+exp(x2+{e2})+...exp(xJ+{eJ}))) /*
*/     (s22=exp(x22+{e2})/(1+exp(x1+{e1})+exp(x2+{e2})+...exp(xJ+{eJ}))) /*
*/ ... (sJJ=exp(xJJ+{eJ})/(1+exp(x1+{e1})+exp(x2+{e2})+...exp(xJ+{eJ})))


STATA returns "cannot have fewer observations than parameters".

Can anyone give me some hint to proceed?

Thanks, Amy

二维码

扫码加我 拉你入群

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

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

关键词:correlations correlation Equations Nonlinear relations matrix

沙发
ohmypig 发表于 2013-9-20 02:47:39 |只看作者 |坛友微信交流群
Unfortunately, "nlsur" doesn't do the job.

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注jltj
拉您入交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-11-9 23:23