Stata新命令：ereplace-egenmore-egenmisc

Stata 连享会   主页 || 视频 || 推文

温馨提示： 定期 清理浏览器缓存，可以获得最佳浏览体验。

✨ 课程主页： https://gitee.com/arlionn/Course

作者：连玉君 (中山大学)
E-Mail: arlionn@163.com

目录

1. 简介

• 强大的 egen 和 egenmore 命令都只能生成新变量，然而，在很多情况下，我们需要替换已有变量。此时就可以采用 ereplace 命令。
• egenmisc 提供了几个新函数：nacorr, nacov， navar 可以计算全样本或分组子样本中变量的相关系数、协方差和方差。pick 可以挑选出符合条件的观察值。上述命令均支持 by(varlist) 选项，这很重要。

2. 安装方法

ssc install ereplace, replace
ssc install egenmisc, replace

3. 说明文档

• egen 命令详解 egen (PDF)

• egenmore 命令简介和补充范例 (功能异常强大)

• egenmisc 命令的 Github 项目主页

4. Stata 范例

4.1 ereplace 命令

help ereplace

  sysuse "nlsw88.dta", clear
  
  egen sd_wage = sd(wage)
  list wage sd_wage in 1/10, clean
  ereplace sd_wage = sd(wage), by(industry)
  list wage sd_wage industry in 1/10, clean

结果

.   sysuse "nlsw88.dta", clear
(NLSW, 1988 extract)

.   
.   egen sd_wage = sd(wage)

.   list wage sd_wage in 1/10, clean

           wage     sd_wage  
  1.   11.73913   5.7555229  
  2.   6.400963   5.7555229  
  3.   5.016723   5.7555229  
  4.   9.033813   5.7555229  
  5.   8.083731   5.7555229  
  6.    4.62963   5.7555229  
  7.   10.49114   5.7555229  
  8.   17.20612   5.7555229  
  9.   13.08374   5.7555229  
 10.   7.745568   5.7555229  

.   ereplace sd_wage = sd(wage), by(industry)
(2,246 real changes made)

.   list wage sd_wage industry in 1/10, clean

           wage     sd_wage                 industry  
  1.   11.73913   6.1277447   Transport/Comm/Utility  
  2.   6.400963   5.3684145            Manufacturing  
  3.   5.016723   5.3684145            Manufacturing  
  4.   9.033813   5.1140115    Professional Services  
  5.   8.083731   5.3684145            Manufacturing  
  6.    4.62963   5.1140115    Professional Services  
  7.   10.49114   6.1277447   Transport/Comm/Utility  
  8.   17.20612   5.1140115    Professional Services  
  9.   13.08374   5.1140115    Professional Services  
 10.   7.745568   5.1140115    Professional Services

4.2 egenmisc 命令

sysuse  "nlsw88.dta", clear
egen temp1 = fastxtile(hours) if married == 0, by(race) nq(10)
egen temp2 = xtile(hours) if married == 0, by(race) nq(10)
assert temp1 == temp2
drop temp*

egen temp1 = sd(wage) if married == 0, by(race) 
egen temp2 = fastsd(wage) if married == 0, by(race) 
assert float(temp1) == float(temp2)
drop temp*

sort race
egen temp1 = corr(hours wage) if married == 0, by(race) 
egen temp2 = fastcorr(hours wage) if married == 0, by(race) 
assert temp1 == temp2 if !missing(hours) & !missing(wage)
drop temp*

相关课程

连享会-直播课 上线了！
http://lianxh.duanshu.com

免费公开课：

• 直击面板数据模型 - 连玉君，时长：1小时40分钟
• Stata 33 讲 - 连玉君, 每讲 15 分钟.
• 部分直播课 课程资料下载 (PPT，dofiles等)

课程一览

支持回看，所有课程可以随时购买观看。

Note: 部分课程的资料，PPT 等可以前往 连享会-直播课 主页查看，下载。

关于我们

• Stata连享会 由中山大学连玉君老师团队创办，定期分享实证分析经验。直播间 有很多视频课程，可以随时观看。
• 连享会-主页 和 知乎专栏，300+ 推文，实证分析不再抓狂。
• 公众号推文分类： 计量专题 | 分类推文 | 资源工具。推文分成 内生性 | 空间计量 | 时序面板 | 结果输出 | 交乘调节 五类，主流方法介绍一目了然：DID, RDD, IV, GMM, FE, Probit 等。
• 公众号关键词搜索/回复 功能已经上线。大家可以在公众号左下角点击键盘图标，输入简要关键词，以便快速呈现历史推文，获取工具软件和数据下载。常见关键词：课程, 直播, 视频, 客服, 模型设定, 研究设计, stata, plus, 绘图, 编程, 面板, 论文重现, 可视化, RDD, DID, PSM, 合成控制法 等

连享会小程序：扫一扫，看推文，看视频……

扫码加入连享会微信群，提问交流更方便

✏ 连享会学习群-常见问题解答汇总：
✨ https://gitee.com/arlionn/WD