悬赏50金币求一个stata程序中一处小错误解释

以下一段stata程序执行到最后提示出错，请告知如何改正，谢谢。
capture program drop hw2question2
program hw2question2
args  lnf a1 a2 a3 a4 a5 a6 a7 g1 g2 g3 g4 g5 g6 g7
tempvar phi1 phi2 phi3 phi4 phi5 phi6 phi7
quietly {
         generate double `phi1' = `a1' * $ML_y2 - `g1'
   generate double `phi2' = (`a1' - `a2') * $ML_y2 - `g1' - `g2'
   generate double `phi3' = (`a1' - `a2' - `a3') * $ML_y2 - `g1' - `g2' - `g3'
   generate double `phi4' = (`a1' - `a2' - `a3' - `a4') * $ML_y2 - `g1' - `g2' - `g3' - `g4'
   generate double `phi5' = (`a1' - `a2' - `a3' - `a4' - `a5') * $ML_y2 - `g1' - `g2' - `g3' - `g4' - `g5'
   generate double `phi6' = (`a1' - `a2' - `a3' - `a4' - `a5'- `a6') * $ML_y2 - `g1' - `g2' - `g3' - `g4' - `g5' - `g6'
   generate double `phi7' = (`a1' - `a2' - `a3' - `a4' - `a5'- `a6' - `a7') * $ML_y2 - `g1' - `g2' - `g3' - `g4' - `g5' - `g6' - `g7'
   replace `lnf' = ln(1 - normal(`phi1')) if $ML_y1 == 0
         replace `lnf' = ln(normal(`phi1') - normal(`phi2')) if $ML_y1 == 1
         replace `lnf' = ln(normal(`phi2') - normal(`phi3')) if $ML_y1 == 2
         replace `lnf' = ln(normal(`phi3') - normal(`phi4')) if $ML_y1 == 3
         replace `lnf' = ln(normal(`phi4') - normal(`phi5')) if $ML_y1 == 4
         replace `lnf' = ln(normal(`phi5') - normal(`phi6')) if $ML_y1 == 5
         replace `lnf' = ln(normal(`phi6') - normal(`phi7')) if $ML_y1 == 6
         replace `lnf' = ln(normal(`phi7')) if $ML_y1 == 7
         }
end
ml model if hw2question2 (nfirm population = ) (a2: ) (a3: ) (a4: ) (a5: ) (a6: ) (a7: ) (g1: ) (g2: ) (g3: ) (g4: ) (g5: ) (g6: ) (g7: )
ml maximize
ml display
program drop hw2question2

执行结果：

. ml model if hw2question2 (nfirm population = ) (a2: ) (a3: ) (a4: ) (a5: ) (a6: ) (a7: ) (g1: ) (g2: ) (g3: ) (g4: ) (g5: ) (g6: ) (g7: )
model required
r(100);
end of do-file

ml model if

不是if而是lf。

另外，最好贴出部分数据，并介绍一下你的任务。

sungmoo 发表于 2012-4-1 06:33
不是if而是lf。

另外，最好贴出部分数据，并介绍一下你的任务。

同意

或许是这里有bug：
ml model lf hw2question2 nfirm population a2 a3 a4 a5 a6 a7 g1 g2 g3 g4 g5 g6 g7

不是if，而应当是lf？
另外，从hw2question2的定义来看，其后面的参数应当是变量。

sungmoo 发表于 2012-3-31 20:12
不是if而是lf。

另外，最好贴出部分数据，并介绍一下你的任务。

谢谢。不过我不知系统何时奉送金币，请稍待。

luhmann 发表于 2012-4-1 08:56
或许是这里有bug：
ml model lf hw2question2 nfirm population a2 a3 a4 a5 a6 a7 g1 g2 g3 g4 g5 g6 g7
...

确实是"if"部分出错，不过由于有人回复在先，悬赏的金币不能相送。但还是很感谢，不知系统是否有可以赠送金币的方式，若有，愿送上20金币。或者以后你有需要下书，可代为下载。谢谢

The problem has been solved. Thank you for all the people responsed.