一般线性方程因变量是twtp,自变量是perdisincome和swc 以及常数项,想用最大似然估计，写出来的程序怎么运行不了？请教下

program my1

args todo b lnfj

tempvar xb

quietly generate double  ‘xb’ = ‘b’[1,1]*perdisincome + ’b’[1,2]*swc + ‘b’[1,3]

quietly replace ‘lnfj’ = ln(normalden((twtp- ‘xb’) / ‘b’[1,4]))- ln(‘b’[1,4])

end

ml model gf0 my1 (twtp=perdisincome swc)/sigma

出现错误提示：unknown method gf0

然后我又改为运行：

ml model lf my1 (twtp=perdisincome swc)/sigma

再运行

ml maximize出现提示错误：' invalid name 完成的程序应该如何写呢 （同时我想估计出sigma） 谢谢！！