Statalist


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

Re: st: r(198) again


From   "Martin Weiss" <Martin.Weiss1@gmx.de>
To   statalist@hsphsun2.harvard.edu
Subject   Re: st: r(198) again
Date   Mon, 16 Mar 2009 18:04:13 +0100

<>

-set more off-, -permanently- if you never want it to bother you again...

HTH
Martin
-------- Original-Nachricht --------
> Datum: Mon, 16 Mar 2009 12:55:48 -0400
> Von: jjc.li@utoronto.ca
> An: statalist@hsphsun2.harvard.edu
> Betreff: Re: st: r(198) again

> 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/

-- 
Martin Weiss
Hackersteigle 3
72076 Tübingen
00497071/793535
00491784597218

Psssst! Schon vom neuen GMX MultiMessenger gehört? Der kann`s mit allen: http://www.gmx.net/de/go/multimessenger01
*
*   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