啊，我明白了，关键的问题在于xtfontier默认的是u--truncated normal 分布，所有还会涉及到mu的检验，汗死！我还以为默认的是half normal呢。多谢四楼高手的指点啊！
好人做到底，高手再多指点我一下吧：）我用xtfrontier命令得到结果后(option为cost)，然后用predict 命令得到TE，但是奇怪的是得到的TE全都是>1，很不可思议，只好每次自己用公式套算te和u。后来我用另一组数据使用froniter命令（option也是cost），也出现了同样的问题，不知道是那一步的命令出错了......
我查看了一下文献，对于product function,te的定义是：te=[f(xi,bi)*e^(vi+ui)]/[f(xi,bi)*e^(vi)]=e^ui，这样，因为对于product function而言，ui总是=0，则te>=1。所以ce=1/te，我又找到xtfront_p.ado文件查看了他的程序，比较了文献上的公式TEi = E[exp(-sui)Iei] ={[1-F(s*sigmai-u/sigmai)]/[1-F(-u/sigmai)]}*exp(-s*u+1/2*sigmai^2)。其中s=1(product function)/=-1(cost function)，u=(-s*e*sigma_u^2+mu*sigma_v^2)/sigmai^2，sigmai^2=sigma_u^2+sigma_v^2，两者公式是一样的。这是不是意味着s的作用只是在估算对数似然统计量和ui的条件期望时起作用，而没有直接得到CE，我得到了TE以后还需要求倒数才是所需要的ce？恳请高手指点一下啊！：）
[此贴子已经被作者于2008-4-13 14:56:01编辑过]