性二型分析

2020-10-16

2020-10-16
性二型分析

  性二型是指在雌雄异体的有性生物中，反映结构和功能变化的某些变量，例如形态特征、行为模式等，在两种性别之间出现固有和明显的差别，使得人们能够以此为根据判断一个个体性别的现象。

  在古生物学中，对于化石性二型的研究有助于帮助古生物学者了解性演化的过程，继而更好地理解贯穿地质历史时期的性选择在生物演化方面所起到的至关重要的作用。然而除了一些保存及其精美的特异埋藏化石外，大多数生物化石既缺乏反映性征的明确证据，如缺少软组织性器官的保存，又不能直接观察到像交配或者繁育等带有鲜明性别差异的性行为，因此在化石研究中，尤其是对于已经灭绝的分类单元或者保存程度较差、数量较少的大化石，个体性别的识别是非常困难的。

  性二型普遍被认为是通过性选择作用演化而来的。通常，介形虫在成虫壳体的形态上表现出性二型的现象，然而以往的研究并没有涉及这种形态性二型（sexual shape dimorphism） 在演化时间尺度上的变化，也没有研究形态性二型与选择压力和幼年个体形态之间的关系。本研究识别出采集自古新世沉积物中的介形类化石Krithe dolichodeira性二型形态特征随时间的变化规律，以讨论性选择作用与性二型之间的相互作用关系。

  注意：

  以下的所有数据中，命名的字段变量中不能存在特殊的符号，例如空格、中文括号、$、#、中文横杆等都需要对数据的表头先进行处理，后续才能对数据进行分析。

1 形态测量与和密度分析

  高斯混合模型就是用高斯概率密度函数（正态分布曲线）精确地量化事物，它是一个将事物分解为若干的基于高斯概率密度函数（正态分布曲线）形成的模型。运用高斯混合模型的前提假设是同一发育阶段的个体呈现正态分布。如下所示：

## [1] 12.31486

  图A表示壳体面积随时间变化的散点图，虚线代表ln（LxH）=12.32 图B为 Krithe dolichodeira成年个体壳体面积的频率分布直方图，曲线是通过高斯混合模型（GMM）估计的密度曲线，利用贝叶斯信息准则（BIC）最低值可得出能将成年个体分为两个亚组分，被视为Krithe dolichodeira的形态变异（morphological variation），分别是morphotype B和morphotype S。为进一步对形态变异进行分析，进行贝叶斯的准则分析。

2 形态变异与贝叶斯准则分析

  贝叶斯信息准则是主观贝叶斯派归纳理论的重要组成部分，就是在不完全情报下，对部分未知的状态用主观概率估计，然后用贝叶斯公式对发生概率进行修正，最后再利用期望值和修正概率做出最优决策。贝叶斯信息准则分为E步骤和M步骤，也就是采用期望值最大化（Expectation Maximization：EM），EM算法就是这样，假设我们估计知道A和B两个参数，在开始状态下二者都是未知的，并且知道了A的信息就可以得到B的信息，反过来知道了B也就得到了A。可以考虑首先赋予A某种初值，以此得到B的估计值，然后从B的当前值出发，重新估计A的取值，这个过程一直持续到收敛为止。如下所示：