R语言编程入门与与绘图简介ppt课件

R语言编程入门与绘图简介,李浩西安交通大学流行病学教研1,R语言历史,R语言是一个自由、免费、开源的软件，是用于统计计算和统计作图的优秀工具。1977年，贝尔实验室开发了S语言，推出了商业软件S-PLIS，R就是免费版的S-PLIS。R为各种数学计算提供计算环境，用户可以根据自己的理解定义所需的数学模型，编程之后存储到“包”中，可以分享给全世界的使用者。所以，R非常自由、灵活；如果说SAS、SPSS是计算机专业的人员做出来的软件，优先考虑算法的实现；那么很多R软件的包都是统计专业的大牛做出的结果，对于本专业的人很容易理解其设计思路。但是，R没有很好的数据库管理器，自带的数据管理器不好用，这是最主要的缺陷。,2,R编辑器,主要有两种，一种是Rstudio，一种是Rgui。Rstudio界面清晰明了，很像matlab，提供了类似的编程语言快速开发环境，能够对代码随时进行调试，是目前最好用的编译器。安装前好像是要安装Rgui?RGui相对界面比较简单，优点是软件启动快，但是界面不是很直观，但是两者可以实现的功能基本没有差别。对于初学者还是建议使用简单明了的Rstudio。菜单类编辑器：RCommander。这个软件可以提供类似SPSS一样的界面和操作，可以实现一些简单的统计模型的实现。Rstudio下载地址：,3,R数据对象,向量a=3,4,R中自定义函数,ex:自定义一个求标准差的函数标准差的计算公式：如下：rm(list=ls()SD-function(x)ss=sum(x-mean(x)2);return(sqrt(ss/length(x)a-1:5;b-SD(a);,5,逻辑语句,选择顺序循环语句在R语言中都可以实现，这里简单说一下if选择语句。SD-function(x)if(length(x)=1)return(error)elsess=sum(x-mean(x)2);return(sqrt(ss/length(x)回顾前面所述：UG1=UG1:2,4:5其实就是一个选择的过程。,6,练习：求分段函数积分,编辑一个求已知分段函数定积分的函数，分段函数如下：=1,(0,12,(1,23,(2,3，求calculus=0；也就是说求不同x时下图的红线左边面积。calculus3|x0Pos-hbv.data$V1;Neg-hbv.data$V2;par(mar=c(4,4,4,4)plot(Age,Pos/(Pos+Neg),cex=0.0005*(Pos+Neg),ylim=c(0,0.7),axes=F,ann=F)axis(side=1,at=seq(0,60,by=10);axis(side=2)title(xlab=Age.merg,ylab=ASPre,main=mergedASPre)lines(Age,Pos/(Pos+Neg),8,图像导出与编辑,作图完成后一般都需要放到word中，这时候需要注意导出为eps(矢量图)格式，这样放到word中才不会发虚。图像导出之后如果需要做一些适当的编辑可以在AI中操作，操作完成依旧可以导入到word中，但是对中文支持不是很好。ex:坐标离断的散点图制作：vac-seq(0,0.8,by=0.1)R0-c(2.0617,1.9265,1.7912,1.6559,1.5205,1.385,1.2494,1.1135,0.9775)plot(vac,R0,type=o,pch=15)text(0.5,1.58,(0.4,1.5205)abline(h=1.5205,v=0.4,lty=3),9,R语言进行蒙特卡罗(MonteCarlo)模拟,蒙特卡罗（MonteCarlo）方法，又称随机抽样或统计试验方法，它是在本世纪四十年代中期为了适应当时原子能事业的发展而发展起来的。是一种通过设定随机过程，反复生成时间序列，计算参数估计量和统计量，进而研究其分布特征的方法。ex:利用MonteCarlo方法求定积分：111222；也就是求下图-1,1范围内曲线下面积。,10,R语言进行蒙特卡罗(MonteCarlo)模拟,思路：首先利用多次抽样求出图形的平均高度，然后用平均高度乘以宽度就是面积了。l-function(n,a,b)g-function(x)1/sqrt(2*pi)*exp(-x2/2)x=runif(n)height=mean(g(a+