I don't know if this is the source of your problem, but in both the successful and unsuccessful nl
statement, the estimated standard deviation, exp([lnsigma]_b[_cons])
, divides only the last [Tier*]_b[race]
term, and then only in the second argument of binormal
(this perhaps because you constrained the SD to be 1
in the first argument). Here's a guess at what you meant in the second nl
statement, but it's just a guess.
#delim ;
_nl_1:
binormal
((1.733373263151318-(.0438395852394023*[Tier1]_b[race])+(0*[Tier1]_b[race]))
/exp([lnsigma]_b[_cons]),
(7.575546508411257-(.0438395852394023*[Tier2]_b[race])+ (0*[Tier2]_b[race]))
/exp([lnsigma]_b[_cons]),
tanh([tau]_b[_cons])) ;
#delim cr