从你的指令来看，并没有计算lnw_n_f。

感谢您的回复， 我发现少贴了一部分结果。这次您能帮我找到问题吗？


sum lnw if female==1
return list
local lnw_mean_f=r(mean)
local lnw_se_f=r(sd)
local lnw_n_f=r(N)
sum lnw if female==0
return list
local lnw_mean_m=r(mean)
local lnw_se_m=r(sd)
local lnw_n_m=r(N)

local `diff_lnw=lnw_mean_m-lnw_mean_f
local `se_diff_lnw=lnw_se_m-lnw_se_f
local `t_diff_lnw=diff_lnw/(lnw_se_f^2/lnw_n_f+lnw_se_m^2/lnw_n_m)^0.5
di "`t_diff_lnw'"

sum data if female==1
return list
local lnw_mean_f=r(mean)
local lnw_se_f=r(sd)
local lnw_n_f=r(N)
sum data if female==0
return list
local lnw_mean_m=r(mean)
local lnw_se_m=r(sd)
local lnw_n_m=r(N)

local diff_lnw=`lnw_mean_m'-`lnw_mean_f'
local se_diff_lnw=`lnw_se_m'-`lnw_se_f'
local t_diff_lnw=`diff_lnw'/(`lnw_se_f'^2/`lnw_n_f'+`lnw_se_m'^2/`lnw_n_m')^0.5
di `t_diff_lnw'

这些运行可以用collapse命令，更方便

谢谢！