用Gumbel分布推算南方水稻基因隔离距离极值.pdfVIP

专 论 与研 兜 中 霸举 米 2010，l6(3)：16—18 用Gumbe1分布推算南方水稻基因隔离距离极值 李建业 t 李仁忠z (浙江省龙游县气象局 ，浙江龙游 324400；浙江省气候中心，浙江杭州 310017) 摘 要：利用Gumbe1分布对南方水稻最大基因隔离距离进行了极值推算 ，得到南方稻区不同熟制水稻的50年 和 100年一遇最大基因隔离距离。从基因隔离距离等值线图可以看出，高值区集中在沿海和东部一带，低值区在内 陆及山区盆地，具有从东到西，从南到北递减或先递减再递增的特征。另外，不同熟制水稻具有晚稻单季稻早稻的 特点。 关键词：Gumbe1分布；基因隔离距离；水稻 过去 10年间，全球转基因作物的种植面积增加了 计算得到基因漂流隔离距离。 50多倍，2005年达到9000万h --。转基因作物在给 1．2 理论方法 人类带来巨大效益的同时 ，也可能给生态环境带来潜 本文采用Gumbe1分布理论来推算极值iT--s]，并用最 在风险，必须通过以科学为基础的安全性评价，提供扎 小二乘法来估算分布中的参数，因为对于样本资料较 实的科学数据和安全管理措施 。转基因水稻基因漂 短的情况，直接用样本平均数和标准差来估算参数，比 流可能产生的环境风险是 目前普遍关注的重点之一 _2_ 炬法准确8[1，比其它方法简单。 4】 。 防止水稻基因逃逸最好的方法是确保水稻制种的纯 将极值样本视为随机变量，可采用耿贝尔分布，其 度，在源头上保证水稻种子的纯度。水稻种子田的隔离 分布形式为 ： 方法之一是空间隔离。空间隔离的距离一般是根据当 )=尸(x x)：exp{一exp一【a(x—b)】J(a0，一o。b+∞) 地的风向和风速来确定。如果同类作物在上风方，隔离 (1) 区就要大一些，反之可以小一些。一般地，空间隔离的 其超过保证率函数为： 距离应不少于200In。但是不同纬度 、不同海拔及不同 x)=1一exp{一exp[一a(x—b)】) (2) 的杂交组合，在制种隔离技术上不能搞一刀切，不能按 可用最小二乘法估算耿贝尔分布的参数。这种方 照一种模式去操作 ，各地应结合当地的气候条件和不 法 ^fx)用类型来进行渐近分布推算。若 n相当大时，有 同的杂交组合对一些技术措施和指标作一些适当的调 近似式 ： 整 。 尸=[a(考 )(x】 ／＼(x) (3) 本文对水稻基因漂流模型计算的基因漂流隔离距 即： 离进行了极值统计 ，用最／1,,-乘法估算 Gumbe1分布参 x)=p(∈．x)=^[a(x—b)]=exp{一exp[一a(x—b)]} (4) 数 ，推算今后50年和 100年一遇的最大基因漂流隔离 因而求考 的渐近分布就归结为由考 的观测值来估计 距离，为制种和食品安全等提供参考。 参数a、b的问题。记为： y=-ln[-lnFn(x)]=a(x_b) (5) 1 资料与方法 而 ：j／(m+1)，In为序列的个数 ，故yi=一In[一ln(Y 1．1 资料 (m+1)】与对应起来 (这里为极值样本从小到大排列的数 本文选用水稻基因漂流模型 计算的结果，取每