豆粕研究员看猪料需求：从大数定律到生猪出栏
伯努利大数定律是说，事件A发生的概率是p，那么假设独立重复的进行n次实验，f_n是n次实验中事件A发生的次数，那么：
也就是说n次实验中事件A发生的次数收敛于单个事件A发生的概率。
从伯努利大数定律到能繁母猪
【能繁母猪】除后备外能够正常生产繁殖的母猪统称能繁母猪，初次配上种的母猪即算作能繁母猪。母猪的怀孕周期大概为4个月，哺乳期大约21天，空怀期大约7天，配种后观察大约24天左右确定妊娠，也就是说两胎之间间隔大约为2个月。如此计算一头母猪一年12月大约能产两胎。
【从伯努利大数定律到能繁母猪】一头能繁母猪，只能有三个状态，怀孕0~1个月，怀孕1~4个月，哺乳+空怀+配种观察期2个月。假设事件A为一头能繁母猪处于怀孕0~1个月的状态，那么事件A发生的概率是1/6，同理，处于怀孕1~4个月的状态发生的概率是1/2，处于哺乳+空怀+配种观察期的状态发生的概率是1/3。根据伯努利大数定律，对我国能繁母猪群体进行抽样实验，事件A发生的概率是1/6，当这个群体足够大的时候，那么这个群体中处于怀孕0~1个月的状态的母猪占比就是1/6。
接下来我们按照这个概率来推算2011年至2017年的生猪出栏量。XX年XX月的出栏量对应着前置10个月的能繁母猪中处于怀孕0~1个月状态的母猪，再根据MSY计算能出栏的育肥猪。根据这个方法推算的出栏量和官方公布的出栏量进行对比的情况：
2011年至2013年拟合的很好，2014年和2017年有差距、差距相对小，但是2015年和2016年推算的数据和官方的数据劈叉较大，这两年官方出栏的数据很大，但细看能繁母猪存栏是呈现比较尖锐的下降趋势的。根据笔者的计算方法，只可能是1/6的系数不对、MSY不对、官方的能繁母猪存栏量不对或者是官方的生猪出栏量数据不对。如果是MSY不对的话，低的能繁母猪存栏和高的出栏量意味着MSY同比大幅增加了，这不太符合事物发展的规律。1/6这个系数不对的话，那又不能解释其它年份拟合程度好的现象。因此可能是能繁母猪或者生猪出栏量统计方面的问题。我国从2014年开始加强禽畜规模养殖防治，在环保高压下，可能存在一些禁养区域的生猪提前出栏，导致出栏量大，以及能繁母猪从南向北转移导致其未纳入统计。
图1 生猪出栏量：根据伯努利大数定律和生猪繁育周期推算
数据来源：wind 中粮期货研究院
图2 生猪出栏量：官方数据与推算数据对比