乌龟性别与温度的LOGIT模型-数学模型.doc

《乌龟性别与温度的LOGIT模型-数学模型.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、数学建模第四次作业第一小组乌龟性别的决定因素问题（习题10——10）钟杰20107122赵新20107134卢昌20107127摘要本文针对幼龟性别温度关系进行建立模型，在建模时建立温度与雄龟比例的模型，由于雄龟比例只可能在【0,1】之间取值，所以不能建立一般的统计回归模型，所以可以建立Logit模型。利用MATLAB统计工具箱中的命令glmfit求解。再求出各个系数估计值后用EXCEL表格进行处理算出每个温度下雄龟比例的估计值，对模型进行验证，为了提高模型的拟合效果可以在Logit模型中添加t的高次方，并用利用MATLAB统计工具箱中的命令regress求解,求出各个系数估计值，置信区

2、间和模型相关系数。为了能直观的比较模型拟合效果可以利用MATLAB中的plot工具会出（t,logit）的散点图，并用利用Tools下BasicFitting工具找出一条拟合效果最好的曲线，并得出其对应的系数。最后可以通过logit基本模型得出当温度为27.7329度幼龟比例为1:1，而且温度每增加一度雄龟与雌龟比例扩大9.12484倍。关键词：logit模型，置信区间，相关系数一、问题提出人类的性别是由基因决定的，乌龟的性别主要有什么因素决定的呢？科学研究表明，决定幼龟性别的最关键的因素是乌龟孵化时的温度。为了研究温度是如何影响幼龟的雌雄比例，美国科学家对某一类乌龟的孵化过程做了实验，

3、试验在五个不同恒定温度下进行，每个温度下分别观察3批乌龟蛋的孵化过程，得到的数据如下：温度乌龟蛋个数雄龟个数雌龟个数雄龟比例27.2101910%8080%91811.1%27.7107370%64266.7%86275%28.313130100%96366.7%87187.5%28.4107370%85362.5%97277.8%29.91110190.9%880100%990100%二、基本假设假设1；幼龟性别只与温度有关三、符号说明符号意义单位备注t乌龟孵化时的温度℃P（x）雄龟比例S乌龟蛋总数个X雄龟个数个C雌龟个数个四、问题分析在本题由于是求温度与性别比例的模型，在数据表中每个

4、温度都记录了三批乌龟，所以首先要对数据进行处理，利用excel算出每个温度对应下乌龟的总数，雄龟总个数，雌龟总个数和雄龟比例。温度乌龟蛋个数雄龟个数雌龟个数雄龟比例27.2272257.41%27.72417770.83%28.33026486.67%28.42719870.37%29.92827196.43%为了使得运算更简单可以把温度进行预处理，把27.2看成0，以此类推可得温度乌龟蛋个数雄龟个数雌龟个数雄龟比例0272257.41%0.52417770.83%1.13026486.67%1.22719870.37%2.72827196.43%为了更直观观察其回归关系，利用MATLA

5、B绘制出散点图。从图中可以看出回归曲线是一条近于3次样条的多项式回归曲线，其回归模型为（1）然而在这个问题中（1）是回归方程中P(x)的取值不一定在[0,1]中，即使P（x）取值在[0,1]中，有意在给定t是，误差项也只能取0,1两个值，显然不具有正态性，而且的方差依赖于与t，具有异方差性，这些都违反了普通回归分析的前提条件，因此，该题不能用用普通回归分析。由于P（x）在[0,1]之间取值，可以使用Logit模型。五、模型的建立与求解5.1模型的建立5.2模型的求解Logit模型是一种广义线性模型，可利用MATLAB统计工具箱总的命令glmfit求解。参数参数估计值标准差B0-10178

6、30.3739B12.21100.4309所以估计值为即所以当幼龟比例为1:1时计算出温度为27.7329度令odd(t)为雄龟与雌龟比例，故有当温度增加1度时odd比为于是由于b1=2.2110,所以温度每增加一度时雄龟与雌龟比例增加到原来的9.12484倍5.3结果的分析及验证对模型各个温度进行验证温度乌龟蛋个数雄龟个数雌龟个数雄龟比例雄龟估计值0272257.41%0.2353570.52417770.83%0.4818121.13026486.67%0.7779551.22719870.37%0.8138012.72827196.43%0.9917685.3模型改进从以上结果可知

7、拟合偏差太大，不适合于做为最终结果。由于模型的右端是温度t的线性函数，可以考虑加入t的二次项后，看是否能提高模型的拟合程度。即考虑模型为执行以下程序>>[b2,dev2]=glmfit([tt.^2],[xs],'binomial','logit');b2,pval=1-chi2cdf(dev-dev2,1)b2=-1.65823.7840-0.7745pval=0.0304计算出b0,b1,b2的估计值为-1.6582、3.784