![Statalist](/includes/images/head_statalist.jpg) |
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
st: r(198) again
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 <[email protected]>:
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 <[email protected]>:
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
Quoting Eva Poen <[email protected]>:
Almost there. Lose the quotes, and call your program directly.
*
* 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/