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

精选,1,R语言编程入门与绘图简介,李浩西安交通大学流行病学教研精选,2,R语言历史,R语言是一个自由、免费、开源的软件，是用于统计计算和统计作图的优秀工具。1977年，贝尔实验室开发了S语言，推出了商业软件S-PLIS，R就是免费版的S-PLIS。R为各种数学计算提供计算环境，用户可以根据自己的理解定义所需的数学模型，编程之后存储到“包”中，可以分享给全世界的使用者。所以，R非常自由、灵活；如果说SAS、SPSS是计算机专业的人员做出来的软件，优先考虑算法的实现；那么很多R软件的包都是统计专业的大牛做出的结果，对于本专业的人很容易理解其设计思路。但是，R没有很好的数据库管理器，自带的数据管理器不好用，这是最主要的缺陷。,精选,3,R编辑器,主要有两种，一种是Rstudio，一种是Rgui。Rstudio界面清晰明了，很像matlab，提供了类似的编程语言快速开发环境，能够对代码随时进行调试，是目前最好用的编译器。安装前好像是要安装Rgui?RGui相对界面比较简单，优点是软件启动快，但是界面不是很直观，但是两者可以实现的功能基本没有差别。对于初学者还是建议使用简单明了的Rstudio。菜单类编辑器：RCommander。这个软件可以提供类似SPSS一样的界面和操作，可以实现一些简单的统计模型的实现。Rstudio下载地址：,精选,4,R数据对象,向量a=3,精选,5,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);,精选,6,逻辑语句,选择顺序循环语句在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其实就是一个选择的过程。,精选,7,练习：求分段函数积分,精选,8,plot()函数绘图简介,R语言提供丰富的绘图包，plot()是一个比较基础的绘图函数。示例如下：setwd(D:/桌面)hbv.data-read.table(D:/桌面/mergdata.dat,header=F)Age-hbv.data$V4;Pos-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),精选,9,图像导出与编辑,作图完成后一般都需要放到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),精选,10,R语言进行蒙特卡罗(MonteCarlo)模拟,精选,11,R语言进行蒙特卡罗(MonteCarlo)模拟,思路：首先利用多次抽样求出图形的平均高度，然后用平均高度乘以宽度就是面积了。l-function(n,a,b)g-function(x)1/sqrt(2*pi)*exp(-x2/2)x=runif(n)height=mean(g(a+(b-a)*x