统计模拟与R相关资料习题答案省公共课一等奖全国赛课获奖课件.pptx
《统计模拟与R相关资料习题答案省公共课一等奖全国赛课获奖课件.pptx》由会员分享,可在线阅读,更多相关《统计模拟与R相关资料习题答案省公共课一等奖全国赛课获奖课件.pptx(162页珍藏版)》请在咨信网上搜索。
1、统计模拟与统计模拟与主讲:朱 强 理学院QQ群群:6052759E-Mail:qiang_Simulation and R第1页课程要求课程要求1.本课程上机实习3-4次2.每次上机需要完成一次试验并提交3.本课程怎么考试?4.总评=考勤+试验作业+考试成绩第2页主要内容主要内容1.R软件介绍2.R统计分析3.统计模拟介绍第3页教材和主要参考书籍教材和主要参考书籍1、肖枝洪、朱强著肖枝洪、朱强著,统计模拟及其统计模拟及其R实现,实现,武汉大学出版社,武汉大学出版社,2、Sheldon.M.Ross,Simulation,(3nd),Elsevier,3、王兆军,陈广雷,邹长亮译,统计模拟王兆军
2、,陈广雷,邹长亮译,统计模拟(第第4版版),人民邮电出版社,人民邮电出版社,.4、薛毅,陈立萍,统计建模与薛毅,陈立萍,统计建模与R软件,软件,清华大学出版社清华大学出版社,.第4页第一篇:第一篇:R介绍介绍R是是一个开放一个开放(GPL)统计编程环境统计编程环境一个语言,是语言一个语言,是语言(由由AT&T Bell试验室试验室Rick Becker,John Chambers,Allan Wilks开发开发)一个一个方言方言(dialect)之一,另一则为之一,另一则为S-plus.一个软件,是集统计分析与图形直观显示于一体统计分一个软件,是集统计分析与图形直观显示于一体统计分析析第5页R
3、作为一个计划作为一个计划(project),最早,最早(1995年年)是由是由Auckland大学统计系大学统计系Robert Gentleman和和Ross Ihaka开始编制,当前由关键开发小组开始编制,当前由关键开发小组(R Development Core Team 以后用以后用R DCT表示表示)维护,他们完全自愿、工作努力负责,维护,他们完全自愿、工作努力负责,并将全球优异统计应用软件打包提供给我们。我们并将全球优异统计应用软件打包提供给我们。我们能够经过计划网站能够经过计划网站(http:/www.r-project.org)了解相关最新信息和使用说明,了解相关最新信息和使用说明
4、,得到最新版本软件和基于应用统计软件包得到最新版本软件和基于应用统计软件包.第6页R无私贡献者Bill Venables Robert Gentleman Ross Ihaka 第7页R无偿无偿 R 资源公开资源公开(不是黑盒子不是黑盒子,也不是吝啬鬼也不是吝啬鬼)R能够在能够在UNIX,Windows和和Macintosh运运行行.R 有优异有优异内在帮助系统内在帮助系统.R有优异有优异画图功效画图功效学生能够轻松地学生能够轻松地转转到商业支持到商业支持 S-Plus程序程序(假假如需要使用商业软件如需要使用商业软件)R语言有一个语言有一个强大强大,轻易学习轻易学习语法语法,有许多内在统有许
5、多内在统计函数计函数.第8页经过用户经过用户自编程序自编程序,R语言很轻易语言很轻易延伸和扩大延伸和扩大.它它就是这么成长就是这么成长.经过语言许多内嵌统计函数,很轻易学习和掌握经过语言许多内嵌统计函数,很轻易学习和掌握R语语言语法言语法.我们能够编制自己函数来扩展现有语言我们能够编制自己函数来扩展现有语言(这就是为何这就是为何它在不停等级完善它在不停等级完善!)R 是计算机编程语言是计算机编程语言.类似于类似于UNIX语言语言,C语言语言,Pascal,Gauss语言等语言等.对于熟练编程者对于熟练编程者,它将以为该语言它将以为该语言比其它语言更比其它语言更熟悉熟悉.而对计算机初学者而对计算
6、机初学者,学习学习R语言语言使得学习下一步使得学习下一步其它编程不那么困难其它编程不那么困难.那些傻瓜软件那些傻瓜软件(SAS,SPSS等等)语言语法则完全语言语法则完全不一样不一样.第9页R缺点缺点不如不如S-Plus在编辑输出画图那样好在编辑输出画图那样好没有商业支持没有商业支持(但有网上支持但有网上支持)需要编程需要编程,不傻瓜不傻瓜.第10页R历史历史S语言在语言在1980年代后期在年代后期在AT&T试验室开发试验室开发.R 项目由项目由Auckland 大学统计系大学统计系Robert Gentleman和和Ross Ihaka于于1995年开始年开始.它很快得到广泛用户欢迎它很快得
7、到广泛用户欢迎.当前它是由当前它是由R关键发关键发展团体维持展团体维持;它是一个由志愿者组成工作努力国它是一个由志愿者组成工作努力国际团体际团体第11页参考资料参考资料 随软件所附随软件所附pdf文档文档(help-manuals),随版随版本更新本更新:W.N.Venables,D.M.Smith and the R DCT:Introduction to R-Notes on R:A Programming Environment for Data Analysis and Graphics,.http:/bayes.math.montana.edu/Rweb/Rnotes/R.htmlR
8、 DCT,The R Environment for Statistical Computing and Graphics-Reference Index,.R DCT,R Data Import/Export,.R DCT,R Language Definition,R DCT,Writing R Extensions,第12页Kurt Hornik,R FAQ,Version 1.8-1,-10-07 B.D.Ripley,R for Windows FAQ,Version for rw1080R Html Help,Statistical Data Analysis 其它其它PDF/HT
9、ML文件:文件:Kickstarting R,http:/cran.r-project.org/doc/contrib/Lemon-kickstart/R examples,Alison Gibbs,http:/www.utstat.toronto.edu/alisong/Teaching/Winter04/Sta248/Rex.html 第13页*Ko-Kang Wang,Introduction to R for Statisticians,J.H.Maindonald,Using R for Data Analysis and Graphics An Introduction,J.H.M
10、aindonald,Using S-PLUS for Data Analysis and Graphics,Emmanuel Paradis,R for beginners Jonathan Baron,R reference card Bret Larget,R for Introductory Statistics,第14页近期畅销书近期畅销书:Peter Dalgaard,Introductory Statistics with R,Springer,John Maindonald,John Braun,Data Analysis and Graphics Using R-An Exam
11、ple-based Approach,Cambridge University Press,John Fox,An R and S-Plus Companion to Applied Regression,Sage Publications,Inc.,W.N.Venables and B.D.Ripley,Modern Applied Statistics with S,4th Ed.,第15页近期畅销书近期畅销书:统计建模与统计建模与R,薛毅等,清华大学出版社,薛毅等,清华大学出版社,R语言统计分析软件教程,王斌会,语言统计分析软件教程,王斌会,R语言与统计分析语言与统计分析,汤银才,高教出
12、版社,汤银才,高教出版社,R语言及语言及Bioconductor在基因组分析中应用,在基因组分析中应用,孙啸等,孙啸等,科学出版社科学出版社第16页O、R安装1.进入网站进入网站http:/www.r-project.org2.点击左边点击左边Download,Packages下下CRAN选择镜像:选择镜像:China有六个镜像点,随意点击一个有六个镜像点,随意点击一个第17页O、R安装3.点击下列图中点击下列图中Windows,再点击再点击base4.点击第一行点击第一行Download R 2.13.1 for Windows(32 megabytes)下载下载第18页一、R运行平台R_G
13、UI开启开启R,我们看到我们看到R GUI(graphic users interface)主窗口主窗口,它由三部分组成它由三部分组成主菜单主菜单工具条工具条R console(R运行窗口运行窗口)R console你主要工作是在这里经过公布命令来完成你主要工作是在这里经过公布命令来完成,包含数据集建包含数据集建立立,数据分析数据分析,作图等作图等.在这里你能够得到在线帮助在这里你能够得到在线帮助help.start()HTML格式关于格式关于R帮助文件帮助文件help()得到对应函数帮助得到对应函数帮助,比如比如help(plot)demo()得到得到R提供几个示例提供几个示例q()退出退出
14、R 同同Matlab类似,用方向键能够重现以前命令类似,用方向键能够重现以前命令第19页二、R语法与数据结构语法语法符号符号 命令或运算提醒符命令或运算提醒符+续行符续行符基本算术运算基本算术运算+加号加号 -减号减号*乘号乘号/除号除号 乘方乘方%*%矩阵乘法矩阵乘法%求余数求余数%/%整数除法整数除法赋值符赋值符=或或 3+53-53/535x=5?plothelp(plot)第21页三、向量三、向量向量是向量是R中最为基本类型中最为基本类型一个向量中元素类型必须相同,包含一个向量中元素类型必须相同,包含数值型数值型整型整型单精度实型单精度实型双精度实型双精度实型逻辑型逻辑型复值型复值型字
15、符型字符型第22页1、数值向量、数值向量定义向量惯用方法定义向量惯用方法(函数函数)seq()或或:若向量若向量(序列序列)含有较为简单规律含有较为简单规律rep()若向量若向量(序列序列)含有较为复杂规律含有较为复杂规律c()若向量若向量(序列序列)没有什么规律没有什么规律例例1:10seq(1,10,by=0.5)seq(1,10,length=21)rep(2:5,2)rep(2:5,rep(2,4)x=c(42,7,64,9)length(x)第23页向量运算向量运算:+,-,*,/,含义是对每个元素进行运算含义是对每个元素进行运算.例例x=1:32*x-111 3 51:表示该行第一
16、个数下标表示该行第一个数下标.%/%:表示整数除法表示整数除法,%表示求余数表示求余数7%/%37%3两个等长向量之间加两个等长向量之间加,减减,乘乘,除除,乘方运算是对应乘方运算是对应元素间四则运算元素间四则运算.第24页1:3+c(5,7,2)/c(5,7,1)12 3 5两个不等长向量之间运算则是两个不等长向量之间运算则是:长度短向量将循长度短向量将循环使用环使用.c(100,200)+c(1,2,3,4)1101 202 103 204第25页向量运算中循环法则(recycling rule)1:2+1:41:4+1:7第26页sqrt(x),log(x),exp(x),sin(x),
17、cos(x),tan(x),abs(x)分别表示平方根分别表示平方根,对数对数,指数指数,三角函数及绝对值三角函数及绝对值.sort(x,decreasing=FALSE)返回按返回按x元素从小到大排元素从小到大排序结果向量序结果向量.order(x)是使得是使得x从小到大排列元素下标向量从小到大排列元素下标向量,sort(x)与与xorder(x)是等效是等效.numeric(n):表示长度为表示长度为n零向量零向量注意:1:n-1与1:(n-1)区分.1:3-11:(3-1)第27页2、逻辑向量逻辑向量向量能够取逻辑值向量能够取逻辑值,如如l=c(TRUE,TRUE,FALSE)x=c(1
18、,4,6.25)l3l1FALSE TRUE TRUE两个向量也可比较两个向量也可比较log(10*x)1 2.302585 3.688879 4.135167log(10*x)x1 TRUE FALSE FALSE比较运算符比较运算符:,=,=(相等相等),!=(不等不等)逻辑向量能够进行与逻辑向量能够进行与(&)表示同时满足表示同时满足,或或(|)二者之一二者之一运算运算.第28页如如:(x=2)&(xall(log(10*x)x)判断是否有真值函数为判断是否有真值函数为anyany(log(10*x)x)函数函数is.na(x)用来判断用来判断x每一个元素是否有缺失值每一个元素是否有缺失
19、值.如如is.na(c(1,2,NA)1 FALSE FALSE TRUE还能够将逻辑值强制转化为整数值,如:将还能够将逻辑值强制转化为整数值,如:将TRUE变成变成1,FALSE变成变成0,c(0,1)(x3)+11 0 1 1(x3)+1#+1表示把逻辑值强制转化为整数值表示把逻辑值强制转化为整数值1,01 1 2 2第29页3、字符型向量字符型向量c1=c(“a”,”b”)c11“a”“b”c2=c(“weight”,”height”,”年纪年纪”)c21”weight”“height”“年纪年纪”paste函数用来把它自变量连成字符串,中间用对应分隔符函数用来把它自变量连成字符串,中间
20、用对应分隔符paste(”My”,”Job”)paste(c(“X”,”Y”),”=“,1:4)paste(“result.”,1:5,sep=“”)#sep表示所用分隔符表示所用分隔符 paste(my,job,sep=)1 myjob paste(my,job,sep=.)1 my.job第30页4、复数向量复数向量复数输入只要像格式复数输入只要像格式5+2i输入即可,以下定义复数向量输入即可,以下定义复数向量x=(0:100)/100*2*piy=sin(x)z=complex(re=x,im=y)complex函数定义复数,函数定义复数,x为实部,为实部,y为虚部为虚部plot(z)画
21、出画出z点图点图也可给定模和辐角定义复数向量:也可给定模和辐角定义复数向量:zz=complex(mod=1,arg=(0:11)/12*2*pi)Re()计算实部,计算实部,Im()计算虚部,计算虚部,Mod()计算复数模,计算复数模,Arg()计算复数辐角,计算复数辐角,Conj()计算共轭。基本数学函数也支持复计算共轭。基本数学函数也支持复数运算。数运算。sqrt(-2+0i)第31页5、向量下标运算向量下标运算R中向量下标从中向量下标从1开始,这与通常统计或数学软件不开始,这与通常统计或数学软件不 一致而象一致而象C语言等计算机高级语言向量下标则从语言等计算机高级语言向量下标则从0开始
22、!开始!访问向量某个元素访问向量某个元素xix=c(42,7,64,9)x2也可单独改变某个元素值也可单独改变某个元素值x3=-1 访问向量一部分几个方法:访问向量一部分几个方法:(1)正整数下标向量正整数下标向量xc(1,4)142 9x2:4(2)负整数下标向量负整数下标向量x-1表示删除第一个元素表示删除第一个元素(3)逻辑下标向量逻辑下标向量xx10表示选取表示选取xxxages=c(Li=23,Zhang=33,Wang=45)ages Li Zhang Wang 23 33 45 访问时能够用通常方法,还能够用元素名访问访问时能够用通常方法,还能够用元素名访问ages“Zhang”
23、Zhang 33还能够定义向量后,再后加上名字还能够定义向量后,再后加上名字age1=c(21,34,56)names(age1)=c(“Zhang”,”Ding”,”Liu”)下面定义一个向量分段函数下面定义一个向量分段函数y=numeric(length(x)yx0=1-xxyx=0=1+xx=0第33页四、矩阵与多维数组四、矩阵与多维数组1、矩阵函数格式、矩阵函数格式matrix(data=NA,nrow=1,ncol=1,byrow=FALSE,dimnames=NULL)data为数组数据向量,为数组数据向量,nrow为函数,为函数,ncol为列数,为列数,byrow=TRUE表示按
24、行写入数据,不然按列。表示按行写入数据,不然按列。R中默认是按列填数据。中默认是按列填数据。A=matrix(1:12,ncol=4,byrow=TRUE)A ,1,2,3,41,1 2 3 42,5 6 7 83,9 10 11 12c(A)显示显示A全部向量,是按列拉直向量。全部向量,是按列拉直向量。cbind()把向量横向拼成一个大矩阵,把向量横向拼成一个大矩阵,rbind()把向量纵向拼把向量纵向拼成一个大矩阵。成一个大矩阵。对角和单位矩阵对角和单位矩阵 diag(1,2,3),diag(1,3)第34页B=rbind(c(1,2),c(3,4)C=cbind(c(11,12),c(1
25、3,14)D=rbind(B,C)E=cbind(B,C)2、访问矩阵元素和子矩阵、访问矩阵元素和子矩阵A2,3访问是矩阵(访问是矩阵(2,3)元素)元素7Ai,访问第访问第i行,行,A,j访问第访问第j列列 表示访问第表示访问第 第第列交叉而成子矩阵列交叉而成子矩阵注意:注意:A=0表示不改变表示不改变A大小只改变元素值为大小只改变元素值为0也能够对矩阵行标和列标重新命名。也能够对矩阵行标和列标重新命名。rownames(A)=c(“a”,”b”,”c”)colnames(A)=paste(“X”,1:4,sep=“”)3、矩阵运算、矩阵运算矩阵能够进行四则运算矩阵能够进行四则运算(+,-,
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 统计 模拟 相关 资料 习题 答案 公共课 一等奖 全国 获奖 课件
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【快乐****生活】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【快乐****生活】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。