用R软件作参数估计_第1页
用R软件作参数估计_第2页
用R软件作参数估计_第3页
用R软件作参数估计_第4页
用R软件作参数估计_第5页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、 数学与计算科学学院统计建模与R软件实验报告题 目： 用R软件作参数估计 专 业： 信息与计算科学 学 号： 姓 名： 指导教师： 成 绩： 二零一五年 四 月 十四日统计建模与R软件实 验 报 告 （ 5 ）学号： 姓名： 班级： 成绩：实验名称：用R软件作参数估计实验地点：所使用的工具软件及环境： R软件一、实验目的：理解参数估计原理，掌握R软件对数据进行参数估计。二、实验步骤： 1、参数估计原理2、能够编写参数估计有关的R函数三、实验内容：1、 矩估计2、极大似然法3、区间估计解：例1.5：编写相应的R函数（程序名：moment_fun.R）monent_fun-function(p)f

2、-c(p1*p2-A1,p1*p2-p1*p22-M2)J-matrix(c(p2,p1,p2-p22,p1-2*p1*p2),nrow=2,byrow=T)list(f=f,J=J)程序名：Newtons.RNewtons-function(fun,x,ep=1e-5,it_max=100)index-0;k-1while(k=it_max)x1-x;obj-fun(x);x-x-solve(obj$J,obj$f);norm-sqrt(x-x1)%*%(x-x1)if(normep)index-1;breakk-k+1obj-fun(x);list(root=x,it=k,index=in

3、dex,FunVal=obj$f)考虑Newton法求解非线性方程组，其中样本取值由随机数产生。建立矩估计的R函数。x-rbinom(100,20,0.7);n A1-mean(x);M2 source(moment_fun.R)source(Newtons.R) source(moment_fun.R) p x f out out求根函数uniroot（），其结果是：$root1 1.$f.root1 -0.$iter1 6$estim.prec1 6.e-05 loglike out out$minimum1 1.$objective1 1400.774例3.1：编写程序建立R文件（interval_estimate1.R）interval_estimate1-function(x,sigma=-1,alpha=0.05)n-length(x);xb=0)tmp-sigma/sqrt(n)*qnorm(1-alpha/2); df-nelsetmp-sd(x)/sqrt(n)*qt(1-alpha/2,n-1);df source(interval_estimate1.R) X interval_estimate1(X,sigma=0.2) mean df a b1 14.96 5 1

人人文库> 全部分类> 教育资料 > 辅导培训

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论