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

1、R语言编程入门与绘图简介李浩西安交通大学流行病学教研室1R语言历史 R语言是一个自由、免费、开源的软件，是用于统计计算和统计作图的优秀工具。 1977年，贝尔实验室开发了S语言，推出了商业软件S-PLIS，R就是免费版的S-PLIS。 R为各种数学计算提供计算环境，用户可以根据自己的理解定义所需的数学模型，编程之后存储到“包”中，可以分享给全世界的使用者。 所以，R非常自由、灵活；如果说SAS、SPSS是计算机专业的人员做出来的软件，优先考虑算法的实现；那么很多R软件的包都是统计专业的大牛做出的结果，对于本专业的人很容易理解其设计思路。 但是，R没有很好的数据库管理器，自带的数据管理器不好用，

2、这是最主要的缺陷。2R编辑器 主要有两种，一种是Rstudio，一种是Rgui。 Rstudio界面清晰明了，很像matlab，提供了类似的编程语言快速开发环境，能够对代码随时进行调试，是目前最好用的编译器。安装前好像是要安装Rgui? RGui相对界面比较简单，优点是软件启动快，但是界面不是很直观，但是两者可以实现的功能基本没有差别。对于初学者还是建议使用简单明了的Rstudio。 菜单类编辑器：R Commander。这个软件可以提供类似SPSS一样的界面和操作，可以实现一些简单的统计模型的实现。 Rstudio下载地址： 3R数据对象 向量 a - c(1,

3、 2, 3, 4, 5); b - 6:10 ; c - 100:90; d - seq(1, 10, by=2); e =34R中自定义函数 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) else ss = sum(x-m

4、ean(x)2); return(sqrt(ss/length(x) 回顾前面所述：UG1 = UG 1:2, 4:5 其实就是一个选择的过程。6练习：求分段函数积分123x07plot()函数绘图简介 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.00

5、05*(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=merged ASPre)lines(Age,Pos/(Pos+Neg)8图像导出与编辑 作图完成后一般都需要放到word中，这时候需要注意导出为eps(矢量图)格式，这样放到word中才不会发虚。 图像导出之后如果需要做一些适当的编辑可以在AI中操作，操作完成依旧可以导入到word中，但是对中文支持不是很好。 ex:坐

6、标离断的散点图制作：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)9R语言进行蒙特卡罗(Monte Carlo)模拟10R语言进行蒙特卡罗(Monte Carlo)模拟 思路：首先利用多次抽样求出图形的平均高度，然后用平均高度乘以宽度就是面积了。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) return(b-a)*height)l(1000, -1, 1)11R语言的灵活性