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