删除了一处明显错误：第一个和第二个孩子的初婚年龄小于等于5的样本。
剩下的几个孩子的样本等用到了再删吧，其实我还是保守了，我感觉初婚年龄小于等于12的都可以删。

然后用第二个孩子跑数据试了试，现在显著了。样本量为8483。

首先是总模型，既包括核心解释变量也包括中介变量：

而且初婚年龄也显著了，教育变量的符号也正确，中专可以提高孩子数量

然后以初婚年龄为因变量，重新跑模型：

又显著了，而且系数符号正确，中专确实可以减小初婚年龄。

然后是只有解释变量和被解释变量的模型结果：

也显著。

因为是部分中介，所以接下来就要检查中介效应占总效应的比例：

先总的看一下三个模型：

参考这位同学：

来自b站upchanges431

因为bd和c是同号，所以通过bd/a来计算，结果大概是20%左右。

接下来是对结果的完善：

首先尝试加入所在县的固定效应：

实现方法

暂时没有找到所在县这个变量，先搁置。

接下来的任务就是：
首先按照周兴的论文把缺的变量补齐（包括剩下的解释变量），看看会不会使结果更稳定。
然后学习2sls回归和probit模型，寻找增加论文丰富度的方法。
然后找方法解决内生性问题。