


全文预览已结束
下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
(a.)Sample-c(3.91,4.85,2.28,4.06,3.70,4.04,5.46,3.53,2.28,1.96,2.53,3.88,2.22,3.47,4.82,2.46,2.99,2.54,0.52,2.50)#The log-likelihood function of the parameter thetalogL-function(theta) L=0 for(x in Sample) L=L+log(1-cos(x-theta)/(2*pi) return(L) #Then plot the curve of the log-likelihood functioncurve(logL,-pi,pi)(b.)The method of moment for theta:xbar-mean(Sample)theta theta1 0.05844061(c.)#The derivative of the log-lokelihood function with respect to thetadlogL-function(theta) dL-0 for(x in Sample) dL=dL-sin(x-theta)/(1-cos(x-theta) return(dL)#The Newton methodNewton - function(f, tol=1E-12,x0=1,N=100) h - 0.001 i - 1 x1 - x0 while (i=N) df_dx - (f(x0+h)-f(x0)/h x1 - (x0 - (f(x0)/df_dx) root - x1 i - i + 1 if (abs(x1-x0) tol) break x0 Newton(dlogL,x0=-2.7)1 -2.6667 Newton(dlogL,x0=2.7)1 2.873095(d.)startpoints roots for(i in 1:200)+ root-Newton(dlogL,x0=startpointsi)+ roots roots unique(roots) 1 -3.0930 -2.7860 -2.6670 -2.5080 -2.3880 -2.2970 -2.2320 -1.6580 -1.4470 -0.9533 -0.0120 0.7906 2.0040 2.2360 2.3610 2.4750 2.5140 2.8730 3.1900 根据不同的迭代初值设置,共收敛到19个不同的根。说明对数似然方程的导数在所在区间内有不止一个零点。 distinstart for(i in 1:length(unique(roots)+ n-which(roots=unique(roots)i)+ distinstart distinstart 1 -3.1415927 -2.7942809 -2.7311333 -2.5732643 -2.3838216 -2.3522478 -2.2259526 -2.1943788 -1.4366077 -1.4050339 -0.8051318 0.5209676 1.9733622 2.225952615 2.2891002 2.4785429 2.5101167 2.5416905 3.015297519个不同的根将迭代初值分为19个小区间,这19个小区间的起始位置如上所示。(e.) min(diff(distinstart)1 0.0315738 which.min(diff(distinstart)1 16 distinstar
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 理财从业考试题目及答案解析
- 考点攻克人教版八年级上册物理机械运动《运动的快慢》专题练习试题(含解析)
- 考点解析-苏科版八年级物理下册《物质的物理属性》重点解析试题
- 难点解析-人教版八年级上册物理光现象《光的反射》专题测评试题(解析卷)
- 预防近视国旗下讲话稿
- 考点攻克人教版八年级上册物理《长度和时间的测量》综合训练练习题(详解)
- 2025年大学《越南语》专业题库- 越南语方言词汇的反义词考辨
- 2025年大学《捷克语》专业题库- 捷克语音乐与舞蹈传统
- 2025年大学《捷克语》专业题库- 捷克语专业学生的社会实践活动
- 计算机工程监理考试题及答案
- 路基路面工程现场检测技术培训课件
- 中考复习《中考经典错题集》课件
- 2型糖尿病防治指南解读共120张课件
- 大学生计算思维调查问卷
- 铁路下穿高速公路安全评价报告
- 《中国民间故事》整本书阅读交流展示课ppt课件(完美版) 小学语文五年级必读书目快乐读书吧
- 运动营养学全套教学课件
- 商旅服务方案(详细版)
- 麻醉机注册技术审查指导原则(2019年 )
- 激光切割机使用说明书
- 《今天是你的生日》(课堂PPT)
评论
0/150
提交评论