




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、案例:头围是反映婴幼儿大脑和颅骨发育程度的重要指标之一,在医学上,对头围的研究具有重要的意义。数据格式如下。根据数据建立头围与年龄间的回归方程令x表示年龄,y表示头围。x和y均为一维变量,同样可以从x和y的散点图上直观地观察它们之间的关系,然后再作进一步的分析从图中可以看出,年龄和头围服从非线性关系,可以考虑做非线性回归分析。根据散点图的走势,可以选取以下函数作为理论回归方程负指数函数双曲线函数幂函数Logistic曲线函数对数函数231xyebbb+=123xyxbbb+=+()312yxbbb=+()3121xyebbb-+=+()123lnyxbbb=+以上函数中都包含多个未知参数,需要
2、由观测数据进行估计,根据需要还可减少或增加未知参数的个数。以上函数都可以呈现出先急速增加,而后区域平缓的趋势,比较适合头围和年龄的观测数据,均可以作为备选的理论回归方程nlinfit函数Matlab统计工具箱中的nlinfit函数用来作一元或多重非线性回归。记Y为因变量观测值向量,矩阵X为自变量观测值矩阵,X的每一列对应一个变量,每一行对应一组观测。为模型的未知参数向量。12nyyYy轾犏犏犏=犏犏犏犏臌( ) ( )( )( ) ()()( ) ()( )111212122212ppnnnpxxxxxxXxxx轾犏犏犏犏=犏犏犏犏臌12nbbbb轾犏犏犏=犏犏犏犏臌beta = nlinfi
3、t(X, y, fun, beta0)返回非线性回归方程式中的未知参数的估计值beta,若回归方程中有k个未知参数,beta为包含k个元素的向量,X为自变量的观测值矩阵,是n*p的矩阵,输入参数y为因变量的观测值向量,是n*1的列向量。输入参数fun是回归方程式所对应的函数的句柄或函数名。应根据回归方程编写如下形式的函数:modelfun为函数名,b为未知参数向量。nlinfit函数的输入参数beta0为用户设定的未知参数的初值,不同的初值可能会有不同的估计结果,故设定初值时最好能够根据实际问题有个提前的预判(),yhatmodelfun b X=beta, r, J, COVB, mse =
4、 nlinfit(X, y, fun, beta0)还返回残差值向量r,雅克比矩阵J,未知参数的协方差矩阵COVB,误差方差2的估计mse(均方误差平方和)。这里的输出可作为其他后续函数的输入,用来计算参数估计值的置信区间,也可用来计算给定x处的预测值及预测值的置信区间。 = nlinfit(X, y, fun, beta0, options)nlinfit函数利用麦夸特(Levenberg-Marquardt)算法求解非线性回归方程中的未知参数,它是一种迭代算法,这里的options用来设定该算法的控制参数。在调用nlinfit函数作一元非线性回归分析之前,应根据观测数据的特点选择合适的理论
5、回归方程,理论回归方程往往不是唯一的,可以有多种选择。有了理论回归方程之后,首先编写理论回归方程对应的M函数。函数应有两个输入参数,一个输出参数。第一个输入为未知参数向量,对于一元回归,第二个输入为自变量观测值向量,而对于多重回归,第二个输入为自变量观测值矩阵。函数的输出为因变量观测值向量,针对所选择的负指数函数,编写M函数如下在种定义方式下,可以将HeadCirl或HeadCirl传递给nlinfitz函数,作为fun参数在种定义方式下,可以将HeadCirl或HeadCirl传递给nlinfitz函数,作为fun参数对于比较简单的理论回归方程,还可以使用符号定义匿名函数,例如可以如下定义负指数函数HeadCir2 = (beta, x)beta(1)*exp(beta(2)./(x+beta(3); 未知参数的选取时一个难点,从散点图上看,随着年龄的增长,人的头围也在增长,但不会一直增长,到了一定的年龄之后,头围就稳定在5055之间。注意到 ,可以选取1的初值为5055之间的一个数,不妨选取为53.2311limxxebbbb+ + =再注意到,初生婴儿的头围在35
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- d-KYFIL-生命科学试剂-MCE
- 固废处理项目可行性分析与发展前景
- 成功项目团队建设的策略试题及答案
- 打造现代农业人才培养新路径与行动方案
- 项目文档管理的试题及答案
- 价值分析2025年特许金融分析师考试试题及答案
- 2025年注册会计师考试学习与应用相结合的策略试题及答案
- 总结经验的证券从业资格试题及答案
- 2025年证券从业资格证考试重要文件解读试题及答案
- 2025年特许金融分析师考试学科交叉分析试题及答案
- 农村急救体系建设
- 仓库搬运工安全操作培训课程
- 广东省地质灾害危险性评估实施细则(2023年修订版)
- 梯子的安全使用课件
- 《非税收入征收管理》课件
- 老年人的口腔知识讲座
- 西格列汀二甲双胍缓释片-药品解读
- 政府采购工作的不足和整改措施
- Unit1+Art+Ancient+Reading+and+Thinking+Chinese+Art+on+show教学设计 高中英语人教选择性必修第三册
- 自驾车出差油费报销单
- 隧道超前地质预报(课件)
评论
0/150
提交评论