Statalist


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

st: r(198) again


From   [email protected]
To   [email protected]
Subject   st: r(198) again
Date   Mon, 16 Mar 2009 12:24:11 -0400

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/



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