STATA 读取PYTHON计算结果

Posted 2023-04-22 myrj

// sfi.Macro提供对Stata宏的访问。
//
// 参数 内容 getGlobal(Name) 获取全局宏的内容 getLocal(Name) 获取局部宏的内容 setGlobal(name, value[, vtype]) 设置全局宏的值 setLocal(name, value) 设置局部宏的值
python
from sfi import Data
from sfi import Datetime
from sfi import Frame
from sfi import Macro
from sfi import ValueLabel
cca="北京"
ac=cca.split(\' \')
for ii in ac:
    Macro.setLocal(\'x\',ii)
end
disp "`x\'"

local jieguo  ""
use cgss\\cgss2015,clear
local k=_N
ds,has(type string)
local bbb="`r(varlist)\'"
local cc ""
foreach var of varlist `bbb\'
    gen ab=0
    replace ab=ab+1 if strrpos(`var\',"`x\'") & (!missing(`var\'))
    sum ab
    if `r(sum)\'>0
        local cc="`cc\'"+" "+"`var\'"
        disp "`cc\'"
        //continue,break
    
    drop ab


disp "`cc\'"

 

Stata输出统计结果到Excel或word

目录

一、安装外部包

二、相关命令

三、实例

1. 描述性统计结果输出

2. 相关性结果输入

3. 回归结果输出

（1）单模型结果

（2）多模型结果

 参考

---

一、安装外部包

在Stata内安装外部包：estout和logout

ssc install estout
ssc install logout

---

二、相关命令

tabstat varlist [if] [in] [weight] [, options]

• 例：tabstat y x1 x2, stat(count mean sd max min) col(stat) format(%10.2f)
• 说明：stat表示要显示的统计量，col(stat)表示行列转置，format表示数据格式及小数点位数

logout, [options : command]

• 例：logout, save(test) excel replace: stat(count mean sd max min)
• 说明：save括号内为文件名，excel可以替换为word，replace:后跟描述性统计的其他命令即可

estout [ what ] [ using filename ] [ , options ]

• 例：esttab using reg.rtf, title(reg) cells(b(star fmt(3)) t(par fmt(3)))
• 说明：输出回归结果，rtf为word可以打开的格式，title为标题，b为截距，t为t值​​​​​​​

esttab [ namelist ] [ using filename ] [ , options ]

• 例：esttab using reg.rtf, r2 ar2 se replace nogap
• 说明：esttab为estout包内的命令，r2表示R方，ar2表示调整的R方，se表示标准误，replace表示每次替代原文件，nogap可去空格行

---

三、实例

下文以lce为被解释变量，lur、lpop、lpg为解释变量

1. 描述性统计结果输出

输入以下命令导出描述性统计结果

logout, save(test1) word replace: tabstat lce lur lpop lpg, stat(count mean sd min max) col(stat) format(%10.2f)

可得到rtf和txt格式保存的表格，如下：

2. 相关性结果输入

输入以下命令导出相关系数结果

logout, save(corr) word replace: pwcorr lce lur lpop lpg

可得到rtf和txt格式保存的表格，如下：

3. 回归结果输出

（1）单模型结果

输出结果之前要先运行回归，接着再使用以下命令

estout, title(reg) cells(b(star fmt(3)) t(par fmt(3)))

 可得到Stata内的结果如图：

若要输出成Word的表格则使用命令esttab

esttab using regn.rtf, r2 ar2 se replace nogap

 可得到Word表格结果如图：

（2）多模型结果

 若有多个模型，则需先保存不同模型的结果，再用esttab

xtreg y x1 x2 x3, fe
est store m1
xtreg y x1 x2 x3 x4, fe
est store m2
esttab m1 m2 using reg2.rtf, r2 ar2 se replace nogap

 可得到Word表格结果如图：

注意：此处显著性水平为5%，1%和0.1%

在命令中加入star(* 0.1 ** 0.05 *** 0.01)即可更改显著性水平

esttab m1 m2 using reg2.rtf, r2 ar2 se star(* 0.1 ** 0.05 *** 0.01) replace nogap

结果如图：

---

 参考

Stata官方help文件