Statalist


[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: st: r(198) again


From   jjc.li@utoronto.ca
To   statalist@hsphsun2.harvard.edu
Subject   Re: st: r(198) again
Date   Mon, 16 Mar 2009 12:55:48 -0400

Is there any way to make stata keep going?

I have to press "more" or "go" so many times.

Jingjing


Quoting Eva Poen <eva.poen@gmail.com>:

I have no idea, really. -set trace on- before the call to -nlsur-, and
see what you can learn from the output roundabout where the error
occurs. It is strange that it happens after 34 successful iterations.

However, by having run

.       nlsurwellbehav lnc sl sk sm lnpl lnpk lnpm lnpe lnq t d1 d2 d3 d4 d5
d6, at (a)

you probably have now overwritten your original left hand side
variables. Don't do this. If you call -nlsurwellbehav- directly, you
_must_ use copies of your left hand side variables.

Therefore, restore your data to their original state, -set trace on-,
call -nlsur wellbehav @ lnc sl sk sm ... - and see where the error
occurs.

Oh, and one more thing: in your parameters() option, a0 should come
last, not first, as it is the last parameter (at position 62) in your
program.

Eva


2009/3/16  <jjc.li@utoronto.ca>:
Hi Eva,

I corrected space before if, then run -nlsur wellbehav @ lnc sl sk sm lnpl
lnpk lnpm lnpe lnq t d1 d2 d3 d4 d5 d6, parameters(a0 al ak am aq aqq blq
bkq bmq bqt bq1 bq2 bq3 bq4 bq5 bq6 ai aii blt bkt bmt ad1 ad2 ad3 ad4 ad5
ad6 bt1 bt2 bt3 bt4 bt5 bt6 bl1 bl2 bl3 bl4 bl5 bl6 bk1 bk2 bk3 bk4 bk5 bk6
 bm1 bm2 bm3 bm4 bm5 bm6 be1 be2 be3 be4 be5 dll dkl dkk dml dmk dmm) neq(4)
ifgnls-

It keeps going and seems OK. However,  when calculating FGNLS estimates,
there's error again:


Iteration 32:  Scaled RSS =  452.8862
Iteration 33:  Scaled RSS =  452.8512
Iteration 34:  Scaled RSS =  452.7397
Iteration 35:  殳 invalid name
r(198);


So I did what you told again:


.       matrix a = J(1,62,0)

.       set trace on

.       nlsurwellbehav lnc sl sk sm lnpl lnpk lnpm lnpe lnq t d1 d2 d3 d4 d5
d6, at (a)
       begin nlsurwellbehav ---
-       version 10.1
-       syntax varlist(min=16 max=16) [if], at(name)
-       tokenize `varlist'
=       tokenize lnc sl sk sm lnpl lnpk lnpm lnpe lnq t d1 d2 d3 d4 d5 d6
-       args lnc sl sk sm lnpl lnpk lnpm lnpe lnq t d1 d2 d3 d4 d5 d6
-       tempname al ak am ae
-       scalar `al'=`at'[1,1]
=       scalar __000000=a[1,1]
-       scalar `ak'=`at'[1,2]
=       scalar __000001=a[1,2]
-       scalar `am'=`at'[1,3]
=       scalar __000002=a[1,3]
-       scalar `ae'=1-`al'-`ak'-`am'
=       scalar __000003=1-__000000-__000001-__000002
-       tempname aq aqq
-       scalar `aq'=`at'[1,4]
=       scalar __000004=a[1,4]
-       scalar `aqq'=`at'[1,5]
=       scalar __000005=a[1,5]
-       tempname blq bkq bmq beq
-       scalar `blq'=`at'[1,6]
=       scalar __000006=a[1,6]
-       scalar `bkq'=`at'[1,7]
=       scalar __000007=a[1,7]
-       scalar `bmq'=`at'[1,8]
=       scalar __000008=a[1,8]
-       scalar `beq'=-`blq'-`bkq'-`bmq'
=       scalar __000009=-__000006-__000007-__000008
-       tempname bqt bq1 bq2 bq3 bq4 bq5 bq6
-       scalar `bqt'=`at'[1,9]
=       scalar __00000A=a[1,9]
-       scalar `bq1'=`at'[1,10]
=       scalar __00000B=a[1,10]
-       scalar `bq2'=`at'[1,11]
=       scalar __00000C=a[1,11]
-       scalar `bq3'=`at'[1,12]
=       scalar __00000D=a[1,12]
-       scalar `bq4'=`at'[1,13]
=       scalar __00000E=a[1,13]
-       scalar `bq5'=`at'[1,14]
=       scalar __00000F=a[1,14]
-       scalar `bq6'=`at'[1,15]
=       scalar __00000G=a[1,15]
-       tempname ai aii blt bkt bmt bet
-       scalar `ai'=`at'[1,16]
=       scalar __00000H=a[1,16]
-       scalar `aii'=`at'[1,17]
=       scalar __00000I=a[1,17]
-       scalar `blt'=`at'[1,18]
=       scalar __00000J=a[1,18]
-       scalar `bkt'=`at'[1,19]
=       scalar __00000K=a[1,19]
-       scalar `bmt'=`at'[1,20]
=       scalar __00000L=a[1,20]
-       scalar `bet'=-`blt'-`bkt'-`bmt'
=       scalar __00000M=-__00000J-__00000K-__00000L
-       tempname ad1 ad2 ad3 ad4 ad5 ad6 bt1 bt2 bt3 bt4 bt5 bt6
-       scalar `ad1'=`at'[1,21]
=       scalar __00000N=a[1,21]
-       scalar `ad2'=`at'[1,22]
=       scalar __00000O=a[1,22]
-       scalar `ad3'=`at'[1,23]
=       scalar __00000P=a[1,23]
-       scalar `ad4'=`at'[1,24]
=       scalar __00000Q=a[1,24]
-       scalar `ad5'=`at'[1,25]
=       scalar __00000R=a[1,25]
-       scalar `ad6'=`at'[1,26]
=       scalar __00000S=a[1,26]
-       scalar `bt1'=`at'[1,27]
=       scalar __00000T=a[1,27]
-       scalar `bt2'=`at'[1,28]
=       scalar __00000U=a[1,28]
-       scalar `bt3'=`at'[1,29]
=       scalar __00000V=a[1,29]
-       scalar `bt4'=`at'[1,30]
=       scalar __00000W=a[1,30]
-       scalar `bt5'=`at'[1,31]
=       scalar __00000X=a[1,31]
-       scalar `bt6'=`at'[1,32]
=       scalar __00000Y=a[1,32]
-       tempname bl1 bl2 bl3 bl4 bl5 bl6 bk1 bk2 bk3 bk4 bk5 bk6 bm1 bm2 bm3
bm4 bm5 bm6 be1 be

       2 be3 be4 be5 be6

-       scalar `bl1'=`at'[1,33]
=       scalar __00000Z=a[1,33]
-       scalar `bl2'=`at'[1,34]
=       scalar __000010=a[1,34]
-       scalar `bl3'=`at'[1,35]
=       scalar __000011=a[1,35]
-       scalar `bl4'=`at'[1,36]
=       scalar __000012=a[1,36]
-       scalar `bl5'=`at'[1,37]
=       scalar __000013=a[1,37]
-       scalar `bl6'=`at'[1,38]
=       scalar __000014=a[1,38]
-       scalar `bk1'=`at'[1,39]
=       scalar __000015=a[1,39]
-       scalar `bk2'=`at'[1,40]
=       scalar __000016=a[1,40]
-       scalar `bk3'=`at'[1,41]
=       scalar __000017=a[1,41]
-       scalar `bk4'=`at'[1,42]
=       scalar __000018=a[1,42]
-       scalar `bk5'=`at'[1,43]
=       scalar __000019=a[1,43]
-       scalar `bk6'=`at'[1,44]
=       scalar __00001A=a[1,44]
-       scalar `bm1'=`at'[1,45]
=       scalar __00001B=a[1,45]
-       scalar `bm2'=`at'[1,46]
=       scalar __00001C=a[1,46]
-       scalar `bm3'=`at'[1,47]
=       scalar __00001D=a[1,47]
-       scalar `bm4'=`at'[1,48]
=       scalar __00001E=a[1,48]
-       scalar `bm5'=`at'[1,49]
=       scalar __00001F=a[1,49]
-       scalar `bm6'=`at'[1,50]
=       scalar __00001G=a[1,50]
-       scalar `be1'=`at'[1,51]
=       scalar __00001H=a[1,51]
-       scalar `be2'=`at'[1,52]
=       scalar __00001I=a[1,52]
-       scalar `be3'=`at'[1,53]
=       scalar __00001J=a[1,53]
-       scalar `be4'=`at'[1,54]
=       scalar __00001K=a[1,54]
-       scalar `be5'=`at'[1,55]
=       scalar __00001L=a[1,55]
-       scalar
`be6'=-`bl1'-`bl2'-`bl3'-`bl4'-`bl5'-`bl6'-`bm1'-`bm2'-`bm3'-`bm4'-`bm5'-`bm6'-`

       bk1'-`bk2'-`bk3'-`bk4'-`bk5'-`bk6'-`be1'-`be2'-`be3'-`be4'-`be5'

=       scalar
__00001M=-__00000Z-__000010-__000011-__000012-__000013-__000014-__00001B-__00001


C-__00001D-__00001E-__00001F-__00001G-__000015-__000016-__000017-__000018-__000019-__00
       001A-__00001H-__00001I-__00001J-__00001K-__00001L

-       tempname dll dkl dkk dml dmk dmm a0
-       scalar `dll'=`at'[1,56]
=       scalar __00001N=a[1,56]
-       scalar `dkl'=`at'[1,57]
=       scalar __00001O=a[1,57]
-       scalar `dkk'=`at'[1,58]
=       scalar __00001P=a[1,58]
-       scalar `dml'=`at'[1,59]
=       scalar __00001Q=a[1,59]
-       scalar `dmk'=`at'[1,60]
=       scalar __00001R=a[1,60]
-       scalar `dmm'=`at'[1,61]
=       scalar __00001S=a[1,61]
-       scalar `a0'=`at'[1,62]
=       scalar __00001T=a[1,62]
-       quietly {
-       replace `lnc'
=`a0'+`aq'*`lnq'+`al'*`lnpl'+`ak'*`lnpk'+`am'*`lnpm'+(1-`al'-`ak'-`am')*`


lnpe'+`blq'*`lnpl'*`lnq'+`bkq'*`lnpk'*`lnq'+`bmq'*`lnpm'*`lnq'+(-`blq'-`bkq'-`bmq')*`ln

pe'*`lnq'+`aqq'*0.5*`lnq'*`lnq'+`ai'*`t'+`aii'*0.5*`t'*`t'+`blt'*`lnpl'*`t'+`bkt'*`lnpk

'*`t'+`bmt'*`lnpm'*`t'+(-`blt'-`bkt'-`bmt')*`lnpe'*`t'+`bqt'*`lnq'*`t'+`ad1'*`d1'+`ad2'

*`d2'+`ad3'*`d3'+`ad4'*`d4'+`ad5'*`d5'+`ad6'*`d6'+`bl1'*`lnpl'*`d1'+`bl2'*`lnpl'*`d2'+`

bl3'*`lnpl'*`d3'+`bl4'*`lnpl'*`d4'+`bl5'*`lnpl'*`d5'+`bl6'*`lnpl'*`d6'+`bk1'*`lnpk'*`d1

'+`bk2'*`lnpk'*`d2'+`bk3'*`lnpk'*`d3'+`bk4'*`lnpk'*`d4'+`bk5'*`lnpk'*`d5'+`bk6'*`lnpk'*

`d6'+`bm1'*`lnpm'*`d1'+`bm2'*`lnpm'*`d2'+`bm3'*`lnpm'*`d3'+`bm4'*`lnpm'*`d4'+`bm5'*`lnp

m'*`d5'+`bm6'*`lnpm'*`d6'+`be1'*`lnpe'*`d1'+`be2'*`lnpe'*`d2'+`be3'*`lnpe'*`d3'+`be4'*`

lnpe'*`d4'+`be5'*`lnpe'*`d5'+(-`bl1'-`bl2'-`bl3'-`bl4'-`bl5'-`bl6'-`bm1'-`bm2'-`bm3'-`b

m4'-`bm5'-`bm6'-`bk1'-`bk2'-`bk3'-`bk4'-`bk5'-`bk6'-`be1'-`be2'-`be3'-`be4'-`be5')*`lnp

e'*`d6'+`bq1'*`lnq'*`d1'+`bq2'*`lnq'*`d2'+`bq3'*`lnq'*`d3'+`bq4'*`lnq'*`d4'+`bq5'*`lnq'

*`d5'+`bq6'*`lnq'*`d6'+`bt1'*`t'*`d1'+`bt2'*`t'*`d2'+`bt3'*`t'*`d3'+`bt4'*`t'*`d4'+`bt5

'*`t'*`d5'+`bt6'*`t'*`d6'+(-`dll'*`dll'+`al'-`al'*`al')*0.5*`lnpl'*`lnpl'+(-`dll'*`dkl'

-`al'*`ak')*`lnpl'*`lnpk'+(-`dll'*`dml'-`al'*`am')*`lnpl'*`lnpm'+(`dll'*`dll'-`al'+`al'

*`al'+`dll'*`dkl'+`al'*`ak'+`dll'*`dml'+`al'*`am')*`lnpl'*`lnpe'+(-`dkl'*`dkl'-`dkk'*`d

kk'+`ak'-`ak'*`ak')*0.5*`lnpk'*`lnpk'+(-`dkl'*`dml'-`dkk'*`dmk'-`ak'*`am')*`lnpk'*`lnpm

'+(`dll'*`dkl'+`al'*`ak'+`dkl'*`dkl'+`dkk'*`dkk'-`ak'+`ak'*`ak'+`dkl'*`dml'+`dkk'*`dmk'

+`ak'*`am')*`lnpk'*`lnpe'+(-`dml'*`dml'-`dmk'*`dmk'-`dmm'*`dmm'+`am'-`am'*`am')*0.5*`ln

pm'*`lnpm'+(`dll'*`dml'+`al'*`am'+`dkl'*`dml'+`dkk'*`dmk'+`ak'*`am'+`dml'*`dml'+`dmk'*`

dmk'+`dmm'*`dmm'-`am'+`am'*`am')*`lnpm'*`lnpe'+(-`dll'*`dll'+`al'-`al'*`al'-`dll'*`dkl'

-`al'*`ak'-`dll'*`dml'-`al'*`am'-`dll'*`dml'-`al'*`am'-`dkl'*`dml'-`dkk'*`dmk'-`ak'*`am

'-`dml'*`dml'-`dmk'*`dmk'-`dmm'*`dmm'+`am'-`am'*`am'-`dll'*`dkl'-`al'*`ak'-`dkl'*`dkl'-

`dkk'*`dkk'+`ak'-`ak'*`ak'-`dkl'*`dml'-`dkk'*`dmk'-`ak'*`am')*0.5*`lnpe'*`lnpe'
`if'

=       replace lnc
=__00001T+__000004*lnq+__000000*lnpl+__000001*lnpk+__000002*lnpm+(1-__00000


0-__000001-__000002)*lnpe+__000006*lnpl*lnq+__000007*lnpk*lnq+__000008*lnpm*lnq+(-__000

006-__000007-__000008)*lnpe*lnq+__000005*0.5*lnq*lnq+__00000H*t+__00000I*0.5*t*t+__0000

0J*lnpl*t+__00000K*lnpk*t+__00000L*lnpm*t+(-__00000J-__00000K-__00000L)*lnpe*t+__00000A

*lnq*t+__00000N*d1+__00000O*d2+__00000P*d3+__00000Q*d4+__00000R*d5+__00000S*d6+__00000Z

*lnpl*d1+__000010*lnpl*d2+__000011*lnpl*d3+__000012*lnpl*d4+__000013*lnpl*d5+__000014*l

npl*d6+__000015*lnpk*d1+__000016*lnpk*d2+__000017*lnpk*d3+__000018*lnpk*d4+__000019*lnp

k*d5+__00001A*lnpk*d6+__00001B*lnpm*d1+__00001C*lnpm*d2+__00001D*lnpm*d3+__00001E*lnpm*

d4+__00001F*lnpm*d5+__00001G*lnpm*d6+__00001H*lnpe*d1+__00001I*lnpe*d2+__00001J*lnpe*d3

+__00001K*lnpe*d4+__00001L*lnpe*d5+(-__00000Z-__000010-__000011-__000012-__000013-__000

014-__00001B-__00001C-__00001D-__00001E-__00001F-__00001G-__000015-__000016-__000017-__

000018-__000019-__00001A-__00001H-__00001I-__00001J-__00001K-__00001L)*lnpe*d6+__00000B

*lnq*d1+__00000C*lnq*d2+__00000D*lnq*d3+__00000E*lnq*d4+__00000F*lnq*d5+__00000G*lnq*d6

+__00000T*t*d1+__00000U*t*d2+__00000V*t*d3+__00000W*t*d4+__00000X*t*d5+__00000Y*t*d6+(-

__00001N*__00001N+__000000-__000000*__000000)*0.5*lnpl*lnpl+(-__00001N*__00001O-__00000

0*__000001)*lnpl*lnpk+(-__00001N*__00001Q-__000000*__000002)*lnpl*lnpm+(__00001N*__0000

1N-__000000+__000000*__000000+__00001N*__00001O+__000000*__000001+__00001N*__00001Q+__0

00000*__000002)*lnpl*lnpe+(-__00001O*__00001O-__00001P*__00001P+__000001-__000001*__000

001)*0.5*lnpk*lnpk+(-__00001O*__00001Q-__00001P*__00001R-__000001*__000002)*lnpk*lnpm+(

__00001N*__00001O+__000000*__000001+__00001O*__00001O+__00001P*__00001P-__000001+__0000

01*__000001+__00001O*__00001Q+__00001P*__00001R+__000001*__000002)*lnpk*lnpe+(-__00001Q

*__00001Q-__00001R*__00001R-__00001S*__00001S+__000002-__000002*__000002)*0.5*lnpm*lnpm

+(__00001N*__00001Q+__000000*__000002+__00001O*__00001Q+__00001P*__00001R+__000001*__00

0002+__00001Q*__00001Q+__00001R*__00001R+__00001S*__00001S-__000002+__000002*__000002)*

lnpm*lnpe+(-__00001N*__00001N+__000000-__000000*__000000-__00001N*__00001O-__000000*__0

00001-__00001N*__00001Q-__000000*__000002-__00001N*__00001Q-__000000*__000002-__00001O*

__00001Q-__00001P*__00001R-__000001*__000002-__00001Q*__00001Q-__00001R*__00001R-__0000

1S*__00001S+__000002-__000002*__000002-__00001N*__00001O-__000000*__000001-__00001O*__0

0001O-__00001P*__00001P+__000001-__000001*__000001-__00001O*__00001Q-__00001P*__00001R-
       __000001*__000002)*0.5*lnpe*lnpe

-       replace `sl' =
`al'+`blq'*`lnq'+(-`dll'+`al'-`al'*`al')*`lnpl'+(-`dll'*`dkl'-`al'*`ak')


*`lnpk'+(-`dll'*`dml'-`al'*`am')*`lnpm'+(`dll'-`al'+`al'*`al'+`dll'*`dkl'+`al'*`ak'+`dl

l'*`dml'+`al'*`am')*`lnpe'+`blt'*`t'+`bl1'*`d1'+`bl2'*`d2'+`bl3'*`d3'+`bl4'*`d4'+`bl5'*
       `d5'+`bl6'*`d6' `if'

=       replace sl =
__000000+__000006*lnq+(-__00001N+__000000-__000000*__000000)*lnpl+(-__0000


1N*__00001O-__000000*__000001)*lnpk+(-__00001N*__00001Q-__000000*__000002)*lnpm+(__0000

1N-__000000+__000000*__000000+__00001N*__00001O+__000000*__000001+__00001N*__00001Q+__0

00000*__000002)*lnpe+__00000J*t+__00000Z*d1+__000010*d2+__000011*d3+__000012*d4+__00001
       3*d5+__000014*d6

-       replace `sk' =
`ak'+`bkq'*`lnq'+(-`dll'*`dkl'-`al'*`ak')*`lnpl'+(-`dkl'*`dkl'-`dkk'*`dk


k'+`ak'-`ak'*`ak')*`lnpk'+(-`dkl'*`dml'-`dkk'*`dmk'-`ak'*`am')*`lnpm'+(`dll'*`dkl'+`al'

*`ak'+`dkl'*`dkl'+`dkk'*`dkk'-`ak'+`ak'*`ak'+`dkl'*`dml'+`dkk'*`dmk'+`ak'*`am')*`lnpe'+

`bkt'*`t'+`bk1'*`d1'+`bk2'*`d2'+`bk3'*`d3'+`bk4'*`d4'+`bk5'*`d5'+`bk6'*`d6'
`if'

=       replace sk =
__000001+__000007*lnq+(-__00001N*__00001O-__000000*__000001)*lnpl+(-__0000


1O*__00001O-__00001P*__00001P+__000001-__000001*__000001)*lnpk+(-__00001O*__00001Q-__00

001P*__00001R-__000001*__000002)*lnpm+(__00001N*__00001O+__000000*__000001+__00001O*__0

0001O+__00001P*__00001P-__000001+__000001*__000001+__00001O*__00001Q+__00001P*__00001R+

__000001*__000002)*lnpe+__00000K*t+__000015*d1+__000016*d2+__000017*d3+__000018*d4+__00
       0019*d5+__00001A*d6

-       replace `sm' =
`am'+`bmq'*`lnq'+(-`dll'*`dml'-`al'*`am')*`lnpl'+(-`dkl'*`dml'-`dkk'*`dm


k'-`ak'*`am')*`lnpk'+(-`dml'*`dml'-`dmk'*`dmk'-`dmm'*`dmm'+`am'-`am'*`am')*`lnpm'+(`dll

'*`dml'+`al'*`am'+`dkl'*`dml'+`dkk'*`dmk'+`ak'*`am'+`dml'*`dml'+`dmk'*`dmk'+`dmm'*`dmm'

-`am'+`am'*`am')*`lnpe'+`bmt'*`t'+`bm1'*`d1'+`bm2'*`d2'+`bm3'*`d3'+`bm4'*`d4'+`bm5'*`d5
       '+`bm6'*`d6' `if'

=       replace sm =
__000002+__000008*lnq+(-__00001N*__00001Q-__000000*__000002)*lnpl+(-__0000


1O*__00001Q-__00001P*__00001R-__000001*__000002)*lnpk+(-__00001Q*__00001Q-__00001R*__00

001R-__00001S*__00001S+__000002-__000002*__000002)*lnpm+(__00001N*__00001Q+__000000*__0

00002+__00001O*__00001Q+__00001P*__00001R+__000001*__000002+__00001Q*__00001Q+__00001R*

__00001R+__00001S*__00001S-__000002+__000002*__000002)*lnpe+__00000L*t+__00001B*d1+__00
       001C*d2+__00001D*d3+__00001E*d4+__00001F*d5+__00001G*d6

-       }
       end nlsurwellbehav ---

What's the problem?

Jingjing















Quoting Eva Poen <eva.poen@gmail.com>:

Good.

This means that it looks like you have found all typos and mistakes.
The only thing I can see that you should correct concerns `if'. There
needs to be a space before the `if' in each of the four -replace-
statements at the end of the program. I.e.

... `bl6'*`d6' `if'

instead of

... `bl6'*`d6'`if'

Other than that, with this version of your -nlsurwellbehav- program,
the actual -nlsur- call should run as intended. For -nlsur- , use your
original left hand side variables, not the copies.

Eva


2009/3/16  <jjc.li@utoronto.ca>:

Hi Eva,

I did it again.


program nlsurwellbehav

       version 10.1

....................

matrix a = J(1,62,0)

gen copy_lnc = lnc

gen copy_sl = sl

gen copy_sm = sm

gen copy_sk = sk

set trace on


nlsurwellbehav copy_lnc copy_sl copy_sk copy_sm lnpl lnpk lnpm lnpe
lnq t d1 d2 d3 d4 d5 d6, at (a)

This time it goes to the end:


       begin nlsurwellbehav ---
-       version 10.1
-       syntax varlist(min=16 max=16) [if], at(name)
-       tokenize `varlist'
=       tokenize copy_lnc copy_sl copy_sk copy_sm lnpl lnpk lnpm lnpe lnq
t
d1 d2 d3 d4 d5 d6
-       args lnc sl sk sm lnpl lnpk lnpm lnpe lnq t d1 d2 d3 d4 d5 d6
-       tempname al ak am ae
-       scalar `al'=`at'[1,1]
=       scalar __00001U=a[1,1]
-       scalar `ak'=`at'[1,2]
=       scalar __00001V=a[1,2]
-       scalar `am'=`at'[1,3]
=       scalar __00001W=a[1,3]
-       scalar `ae'=1-`al'-`ak'-`am'
=       scalar __00001X=1-__00001U-__00001V-__00001W
-       tempname aq aqq
-       scalar `aq'=`at'[1,4]
=       scalar __00001Y=a[1,4]
-       scalar `aqq'=`at'[1,5]
=       scalar __00001Z=a[1,5]
-       tempname blq bkq bmq beq
-       scalar `blq'=`at'[1,6]
=       scalar __000020=a[1,6]
-       scalar `bkq'=`at'[1,7]
=       scalar __000021=a[1,7]
-       scalar `bmq'=`at'[1,8]
=       scalar __000022=a[1,8]
-       scalar `beq'=-`blq'-`bkq'-`bmq'
=       scalar __000023=-__000020-__000021-__000022
-       tempname bqt bq1 bq2 bq3 bq4 bq5 bq6
-       scalar `bqt'=`at'[1,9]
=       scalar __000024=a[1,9]
-       scalar `bq1'=`at'[1,10]
=       scalar __000025=a[1,10]
-       scalar `bq2'=`at'[1,11]
=       scalar __000026=a[1,11]
-       scalar `bq3'=`at'[1,12]
=       scalar __000027=a[1,12]
-       scalar `bq4'=`at'[1,13]
=       scalar __000028=a[1,13]
-       scalar `bq5'=`at'[1,14]
=       scalar __000029=a[1,14]
-       scalar `bq6'=`at'[1,15]
=       scalar __00002A=a[1,15]
-       tempname ai aii blt bkt bmt bet
-       scalar `ai'=`at'[1,16]
=       scalar __00002B=a[1,16]
-       scalar `aii'=`at'[1,17]
=       scalar __00002C=a[1,17]
-       scalar `blt'=`at'[1,18]
=       scalar __00002D=a[1,18]
-       scalar `bkt'=`at'[1,19]
=       scalar __00002E=a[1,19]
-       scalar `bmt'=`at'[1,20]
=       scalar __00002F=a[1,20]
-       scalar `bet'=-`blt'-`bkt'-`bmt'
=       scalar __00002G=-__00002D-__00002E-__00002F
-       tempname ad1 ad2 ad3 ad4 ad5 ad6 bt1 bt2 bt3 bt4 bt5 bt6
-       scalar `ad1'=`at'[1,21]
=       scalar __00002H=a[1,21]
-       scalar `ad2'=`at'[1,22]
=       scalar __00002I=a[1,22]
-       scalar `ad3'=`at'[1,23]
=       scalar __00002J=a[1,23]
-       scalar `ad4'=`at'[1,24]
=       scalar __00002K=a[1,24]
-       scalar `ad5'=`at'[1,25]
=       scalar __00002L=a[1,25]
-       scalar `ad6'=`at'[1,26]
=       scalar __00002M=a[1,26]
-       scalar `bt1'=`at'[1,27]
=       scalar __00002N=a[1,27]
-       scalar `bt2'=`at'[1,28]
=       scalar __00002O=a[1,28]
-       scalar `bt3'=`at'[1,29]
=       scalar __00002P=a[1,29]
-       scalar `bt4'=`at'[1,30]
=       scalar __00002Q=a[1,30]
-       scalar `bt5'=`at'[1,31]
=       scalar __00002R=a[1,31]
-       scalar `bt6'=`at'[1,32]
=       scalar __00002S=a[1,32]
-       tempname bl1 bl2 bl3 bl4 bl5 bl6 bk1 bk2 bk3 bk4 bk5 bk6 bm1 bm2
bm3
bm4 bm5 bm6 be1 be

       2 be3 be4 be5 be6

-       scalar `bl1'=`at'[1,33]
=       scalar __00002T=a[1,33]
-       scalar `bl2'=`at'[1,34]
=       scalar __00002U=a[1,34]
-       scalar `bl3'=`at'[1,35]
=       scalar __00002V=a[1,35]
-       scalar `bl4'=`at'[1,36]
=       scalar __00002W=a[1,36]
-       scalar `bl5'=`at'[1,37]
=       scalar __00002X=a[1,37]
-       scalar `bl6'=`at'[1,38]
=       scalar __00002Y=a[1,38]
-       scalar `bk1'=`at'[1,39]
=       scalar __00002Z=a[1,39]
-       scalar `bk2'=`at'[1,40]
=       scalar __000030=a[1,40]
-       scalar `bk3'=`at'[1,41]
=       scalar __000031=a[1,41]
-       scalar `bk4'=`at'[1,42]
=       scalar __000032=a[1,42]
-       scalar `bk5'=`at'[1,43]
=       scalar __000033=a[1,43]
-       scalar `bk6'=`at'[1,44]
=       scalar __000034=a[1,44]
-       scalar `bm1'=`at'[1,45]
=       scalar __000035=a[1,45]
-       scalar `bm2'=`at'[1,46]
=       scalar __000036=a[1,46]
-       scalar `bm3'=`at'[1,47]
=       scalar __000037=a[1,47]
-       scalar `bm4'=`at'[1,48]
=       scalar __000038=a[1,48]
-       scalar `bm5'=`at'[1,49]
=       scalar __000039=a[1,49]
-       scalar `bm6'=`at'[1,50]
=       scalar __00003A=a[1,50]
-       scalar `be1'=`at'[1,51]
=       scalar __00003B=a[1,51]
-       scalar `be2'=`at'[1,52]
=       scalar __00003C=a[1,52]
-       scalar `be3'=`at'[1,53]
=       scalar __00003D=a[1,53]
-       scalar `be4'=`at'[1,54]
=       scalar __00003E=a[1,54]
-       scalar `be5'=`at'[1,55]
=       scalar __00003F=a[1,55]
-       scalar

`be6'=-`bl1'-`bl2'-`bl3'-`bl4'-`bl5'-`bl6'-`bm1'-`bm2'-`bm3'-`bm4'-`bm5'-`bm6'-`

       bk1'-`bk2'-`bk3'-`bk4'-`bk5'-`bk6'-`be1'-`be2'-`be3'-`be4'-`be5'

=       scalar

__00003G=-__00002T-__00002U-__00002V-__00002W-__00002X-__00002Y-__000035-__00003



6-__000037-__000038-__000039-__00003A-__00002Z-__000030-__000031-__000032-__000033-__00
       0034-__00003B-__00003C-__00003D-__00003E-__00003F

-       tempname dll dkl dkk dml dmk dmm a0
-       scalar `dll'=`at'[1,56]
=       scalar __00003H=a[1,56]
-       scalar `dkl'=`at'[1,57]
=       scalar __00003I=a[1,57]
-       scalar `dkk'=`at'[1,58]
=       scalar __00003J=a[1,58]
-       scalar `dml'=`at'[1,59]
=       scalar __00003K=a[1,59]
-       scalar `dmk'=`at'[1,60]
=       scalar __00003L=a[1,60]
-       scalar `dmm'=`at'[1,61]
=       scalar __00003M=a[1,61]
-       scalar `a0'=`at'[1,62]
=       scalar __00003N=a[1,62]
-       quietly {
-       replace `lnc'
=`a0'+`aq'*`lnq'+`al'*`lnpl'+`ak'*`lnpk'+`am'*`lnpm'+(1-`al'-`ak'-`am')*`



lnpe'+`blq'*`lnpl'*`lnq'+`bkq'*`lnpk'*`lnq'+`bmq'*`lnpm'*`lnq'+(-`blq'-`bkq'-`bmq')*`ln


pe'*`lnq'+`aqq'*0.5*`lnq'*`lnq'+`ai'*`t'+`aii'*0.5*`t'*`t'+`blt'*`lnpl'*`t'+`bkt'*`lnpk


'*`t'+`bmt'*`lnpm'*`t'+(-`blt'-`bkt'-`bmt')*`lnpe'*`t'+`bqt'*`lnq'*`t'+`ad1'*`d1'+`ad2'


*`d2'+`ad3'*`d3'+`ad4'*`d4'+`ad5'*`d5'+`ad6'*`d6'+`bl1'*`lnpl'*`d1'+`bl2'*`lnpl'*`d2'+`


bl3'*`lnpl'*`d3'+`bl4'*`lnpl'*`d4'+`bl5'*`lnpl'*`d5'+`bl6'*`lnpl'*`d6'+`bk1'*`lnpk'*`d1


'+`bk2'*`lnpk'*`d2'+`bk3'*`lnpk'*`d3'+`bk4'*`lnpk'*`d4'+`bk5'*`lnpk'*`d5'+`bk6'*`lnpk'*


`d6'+`bm1'*`lnpm'*`d1'+`bm2'*`lnpm'*`d2'+`bm3'*`lnpm'*`d3'+`bm4'*`lnpm'*`d4'+`bm5'*`lnp


m'*`d5'+`bm6'*`lnpm'*`d6'+`be1'*`lnpe'*`d1'+`be2'*`lnpe'*`d2'+`be3'*`lnpe'*`d3'+`be4'*`


lnpe'*`d4'+`be5'*`lnpe'*`d5'+(-`bl1'-`bl2'-`bl3'-`bl4'-`bl5'-`bl6'-`bm1'-`bm2'-`bm3'-`b


m4'-`bm5'-`bm6'-`bk1'-`bk2'-`bk3'-`bk4'-`bk5'-`bk6'-`be1'-`be2'-`be3'-`be4'-`be5')*`lnp


e'*`d6'+`bq1'*`lnq'*`d1'+`bq2'*`lnq'*`d2'+`bq3'*`lnq'*`d3'+`bq4'*`lnq'*`d4'+`bq5'*`lnq'


*`d5'+`bq6'*`lnq'*`d6'+`bt1'*`t'*`d1'+`bt2'*`t'*`d2'+`bt3'*`t'*`d3'+`bt4'*`t'*`d4'+`bt5


'*`t'*`d5'+`bt6'*`t'*`d6'+(-`dll'*`dll'+`al'-`al'*`al')*0.5*`lnpl'*`lnpl'+(-`dll'*`dkl'


-`al'*`ak')*`lnpl'*`lnpk'+(-`dll'*`dml'-`al'*`am')*`lnpl'*`lnpm'+(`dll'*`dll'-`al'+`al'


*`al'+`dll'*`dkl'+`al'*`ak'+`dll'*`dml'+`al'*`am')*`lnpl'*`lnpe'+(-`dkl'*`dkl'-`dkk'*`d


kk'+`ak'-`ak'*`ak')*0.5*`lnpk'*`lnpk'+(-`dkl'*`dml'-`dkk'*`dmk'-`ak'*`am')*`lnpk'*`lnpm


'+(`dll'*`dkl'+`al'*`ak'+`dkl'*`dkl'+`dkk'*`dkk'-`ak'+`ak'*`ak'+`dkl'*`dml'+`dkk'*`dmk'


+`ak'*`am')*`lnpk'*`lnpe'+(-`dml'*`dml'-`dmk'*`dmk'-`dmm'*`dmm'+`am'-`am'*`am')*0.5*`ln


pm'*`lnpm'+(`dll'*`dml'+`al'*`am'+`dkl'*`dml'+`dkk'*`dmk'+`ak'*`am'+`dml'*`dml'+`dmk'*`


dmk'+`dmm'*`dmm'-`am'+`am'*`am')*`lnpm'*`lnpe'+(-`dll'*`dll'+`al'-`al'*`al'-`dll'*`dkl'


-`al'*`ak'-`dll'*`dml'-`al'*`am'-`dll'*`dml'-`al'*`am'-`dkl'*`dml'-`dkk'*`dmk'-`ak'*`am


'-`dml'*`dml'-`dmk'*`dmk'-`dmm'*`dmm'+`am'-`am'*`am'-`dll'*`dkl'-`al'*`ak'-`dkl'*`dkl'-


`dkk'*`dkk'+`ak'-`ak'*`ak'-`dkl'*`dml'-`dkk'*`dmk'-`ak'*`am')*0.5*`lnpe'*`lnpe'`if'

=       replace copy_lnc
=__00003N+__00001Y*lnq+__00001U*lnpl+__00001V*lnpk+__00001W*lnpm+(1-__



00001U-__00001V-__00001W)*lnpe+__000020*lnpl*lnq+__000021*lnpk*lnq+__000022*lnpm*lnq+(-


__000020-__000021-__000022)*lnpe*lnq+__00001Z*0.5*lnq*lnq+__00002B*t+__00002C*0.5*t*t+_


_00002D*lnpl*t+__00002E*lnpk*t+__00002F*lnpm*t+(-__00002D-__00002E-__00002F)*lnpe*t+__0


00024*lnq*t+__00002H*d1+__00002I*d2+__00002J*d3+__00002K*d4+__00002L*d5+__00002M*d6+__0


0002T*lnpl*d1+__00002U*lnpl*d2+__00002V*lnpl*d3+__00002W*lnpl*d4+__00002X*lnpl*d5+__000


02Y*lnpl*d6+__00002Z*lnpk*d1+__000030*lnpk*d2+__000031*lnpk*d3+__000032*lnpk*d4+__00003


3*lnpk*d5+__000034*lnpk*d6+__000035*lnpm*d1+__000036*lnpm*d2+__000037*lnpm*d3+__000038*


lnpm*d4+__000039*lnpm*d5+__00003A*lnpm*d6+__00003B*lnpe*d1+__00003C*lnpe*d2+__00003D*ln


pe*d3+__00003E*lnpe*d4+__00003F*lnpe*d5+(-__00002T-__00002U-__00002V-__00002W-__00002X-


__00002Y-__000035-__000036-__000037-__000038-__000039-__00003A-__00002Z-__000030-__0000


31-__000032-__000033-__000034-__00003B-__00003C-__00003D-__00003E-__00003F)*lnpe*d6+__0


00025*lnq*d1+__000026*lnq*d2+__000027*lnq*d3+__000028*lnq*d4+__000029*lnq*d5+__00002A*l


nq*d6+__00002N*t*d1+__00002O*t*d2+__00002P*t*d3+__00002Q*t*d4+__00002R*t*d5+__00002S*t*


d6+(-__00003H*__00003H+__00001U-__00001U*__00001U)*0.5*lnpl*lnpl+(-__00003H*__00003I-__


00001U*__00001V)*lnpl*lnpk+(-__00003H*__00003K-__00001U*__00001W)*lnpl*lnpm+(__00003H*_


_00003H-__00001U+__00001U*__00001U+__00003H*__00003I+__00001U*__00001V+__00003H*__00003


K+__00001U*__00001W)*lnpl*lnpe+(-__00003I*__00003I-__00003J*__00003J+__00001V-__00001V*


__00001V)*0.5*lnpk*lnpk+(-__00003I*__00003K-__00003J*__00003L-__00001V*__00001W)*lnpk*l


npm+(__00003H*__00003I+__00001U*__00001V+__00003I*__00003I+__00003J*__00003J-__00001V+_


_00001V*__00001V+__00003I*__00003K+__00003J*__00003L+__00001V*__00001W)*lnpk*lnpe+(-__0


0003K*__00003K-__00003L*__00003L-__00003M*__00003M+__00001W-__00001W*__00001W)*0.5*lnpm


*lnpm+(__00003H*__00003K+__00001U*__00001W+__00003I*__00003K+__00003J*__00003L+__00001V


*__00001W+__00003K*__00003K+__00003L*__00003L+__00003M*__00003M-__00001W+__00001W*__000


01W)*lnpm*lnpe+(-__00003H*__00003H+__00001U-__00001U*__00001U-__00003H*__00003I-__00001


U*__00001V-__00003H*__00003K-__00001U*__00001W-__00003H*__00003K-__00001U*__00001W-__00


003I*__00003K-__00003J*__00003L-__00001V*__00001W-__00003K*__00003K-__00003L*__00003L-_


_00003M*__00003M+__00001W-__00001W*__00001W-__00003H*__00003I-__00001U*__00001V-__00003


I*__00003I-__00003J*__00003J+__00001V-__00001V*__00001V-__00003I*__00003K-__00003J*__00
       003L-__00001V*__00001W)*0.5*lnpe*lnpe

-       replace `sl' =
`al'+`blq'*`lnq'+(-`dll'+`al'-`al'*`al')*`lnpl'+(-`dll'*`dkl'-`al'*`ak')



*`lnpk'+(-`dll'*`dml'-`al'*`am')*`lnpm'+(`dll'-`al'+`al'*`al'+`dll'*`dkl'+`al'*`ak'+`dl


l'*`dml'+`al'*`am')*`lnpe'+`blt'*`t'+`bl1'*`d1'+`bl2'*`d2'+`bl3'*`d3'+`bl4'*`d4'+`bl5'*
       `d5'+`bl6'*`d6'`if'

=       replace copy_sl =
__00001U+__000020*lnq+(-__00003H+__00001U-__00001U*__00001U)*lnpl+(-_



_00003H*__00003I-__00001U*__00001V)*lnpk+(-__00003H*__00003K-__00001U*__00001W)*lnpm+(_


_00003H-__00001U+__00001U*__00001U+__00003H*__00003I+__00001U*__00001V+__00003H*__00003


K+__00001U*__00001W)*lnpe+__00002D*t+__00002T*d1+__00002U*d2+__00002V*d3+__00002W*d4+__
       00002X*d5+__00002Y*d6

-       replace `sk' =
`ak'+`bkq'*`lnq'+(-`dll'*`dkl'-`al'*`ak')*`lnpl'+(-`dkl'*`dkl'-`dkk'*`dk



k'+`ak'-`ak'*`ak')*`lnpk'+(-`dkl'*`dml'-`dkk'*`dmk'-`ak'*`am')*`lnpm'+(`dll'*`dkl'+`al'


*`ak'+`dkl'*`dkl'+`dkk'*`dkk'-`ak'+`ak'*`ak'+`dkl'*`dml'+`dkk'*`dmk'+`ak'*`am')*`lnpe'+


`bkt'*`t'+`bk1'*`d1'+`bk2'*`d2'+`bk3'*`d3'+`bk4'*`d4'+`bk5'*`d5'+`bk6'*`d6'`if'

=       replace copy_sk =
__00001V+__000021*lnq+(-__00003H*__00003I-__00001U*__00001V)*lnpl+(-_



_00003I*__00003I-__00003J*__00003J+__00001V-__00001V*__00001V)*lnpk+(-__00003I*__00003K


-__00003J*__00003L-__00001V*__00001W)*lnpm+(__00003H*__00003I+__00001U*__00001V+__00003


I*__00003I+__00003J*__00003J-__00001V+__00001V*__00001V+__00003I*__00003K+__00003J*__00


003L+__00001V*__00001W)*lnpe+__00002E*t+__00002Z*d1+__000030*d2+__000031*d3+__000032*d4
       +__000033*d5+__000034*d6

-       replace `sm' =
`am'+`bmq'*`lnq'+(-`dll'*`dml'-`al'*`am')*`lnpl'+(-`dkl'*`dml'-`dkk'*`dm



k'-`ak'*`am')*`lnpk'+(-`dml'*`dml'-`dmk'*`dmk'-`dmm'*`dmm'+`am'-`am'*`am')*`lnpm'+(`dll


'*`dml'+`al'*`am'+`dkl'*`dml'+`dkk'*`dmk'+`ak'*`am'+`dml'*`dml'+`dmk'*`dmk'+`dmm'*`dmm'


-`am'+`am'*`am')*`lnpe'+`bmt'*`t'+`bm1'*`d1'+`bm2'*`d2'+`bm3'*`d3'+`bm4'*`d4'+`bm5'*`d5
       '+`bm6'*`d6'`if'

=       replace copy_sm =
__00001W+__000022*lnq+(-__00003H*__00003K-__00001U*__00001W)*lnpl+(-_



_00003I*__00003K-__00003J*__00003L-__00001V*__00001W)*lnpk+(-__00003K*__00003K-__00003L


*__00003L-__00003M*__00003M+__00001W-__00001W*__00001W)*lnpm+(__00003H*__00003K+__00001


U*__00001W+__00003I*__00003K+__00003J*__00003L+__00001V*__00001W+__00003K*__00003K+__00


003L*__00003L+__00003M*__00003M-__00001W+__00001W*__00001W)*lnpe+__00002F*t+__000035*d1
       +__000036*d2+__000037*d3+__000038*d4+__000039*d5+__00003A*d6

-       }
       end nlsurwellbehav ---

.

Jingjing

*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/





*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/



© Copyright 1996–2014 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index