# Stata：用负二项分布预测蚊子存活率

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

## 1. 数据生成

``````clear
set obs 12000
``````

``````gen mort_rate=mod(_n,6)/200+0.01
tab mort_rate
``````

``````gen living=1
``````

``````gen survived=0
``````

## 2. 数值模拟

``````forvalue i=1/500{
*随机模拟伯努利试验，生成临时变量`died`表示特定某一天蚊子是否死亡
qui gen died=rbinomial(1,mort_rate) if living==1
*如果蚊子死亡，则替换其存活状态变量
qui replace living=0 if living==1 & died == 1
*如果蚊子还活着，替换survived变量来表示当前是第几天
qui replace survived=`i' if living==1
*每次循环报告出当前循环是第几天，以及当天蚊子尚存活的比例
qui sum living
di "Round `i' :" r(mean)
*删除临时变量`died`
drop died
}
``````

### 2.1 绘制图像

``````hist survived, by(mort_rate)
``````

### 2.2 不同死亡率下的平均存活率

``````bysort mort_rate: sum survived
``````
``````-----------------------------------------------------------------------
-> mort_rate = .01

Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
survived |      2,000     100.424     96.9187          0        500

-----------------------------------------------------------------------
-> mort_rate = .015

Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
survived |      2,000     66.1455    67.45336          0        470

-----------------------------------------------------------------------
-> mort_rate = .02

Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
survived |      2,000      49.517    50.21277          0        408

-----------------------------------------------------------------------
-> mort_rate = .025

Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
survived |      2,000        40.9    40.77116          0        302

-----------------------------------------------------------------------
-> mort_rate = .03

Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
survived |      2,000     32.2755    33.19746          0        276

-----------------------------------------------------------------------
-> mort_rate = .035

Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
survived |      2,000      28.033    28.43426          0        215

``````

## 相关课程

http://lianxh.duanshu.com

### 课程一览

Stata数据清洗 游万海 直播, 2 小时，已上线

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

#### 关于我们

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

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