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

R语言编程入门与绘图简介李浩西安交通大学流行病学教研163.comR语言历史R语言是一个自由、免费、开源的软件，是用于统计计算和统计作图的优秀工具。1977年，贝尔实验室开发了S语言，推出了商业软件S-PLIS，R就是免费版的S-PLIS。R为各种数学计算提供计算环境，用户可以根据自己的理解定义所需的数学模型，编程之后存储到“包”中，可以分享给全世界的使用者。所以，R非常自由、灵活；如果说SAS、SPSS是计算机专业的人员做出来的软件，优先考虑算法的实现；那么很多R软件的包都是统计专业的大牛做出的结果，对于本专业的人很容易理解其设计思路。但是，R没有很好的数据库管理器，自带的数据管理器不好用，这是最主要的缺陷。R编辑器主要有两种，一种是Rstudio，一种是Rgui。Rstudio界面清晰明了，很像matlab，提供了类似的编程语言快速开发环境，能够对代码随时进行调试，是目前最好用的编译器。安装前好像是要安装Rgui??RGui相对界面比较简单，优点是软件启动快，但是界面不是很直观，但是两者可以实现的功能基本没有差别。对于初学者还是建议使用简单明了的Rstudio。菜单类编辑器：R Commander。这个软件可以提供类似SPSS一样的界面和操作，可以实现一些简单的统计模型的实现。Rstudio下载地址： R数据对象向量 a - c(1, 2, 3, 4, 5); b - 6:10 ; c - 100:90; d - seq(1, 10, by=2); e - seq(1, 10, length=20);数据框 数据框可以看作是多个向量所组成的一种数据集合。 UG = data.frame(a, b[1:5], c[1:5], d[1:5], e[1:5]); UG1 = UG[ 1:2, 4:5]; UG2 = UG$a+UG$b; UG3 = UG[ , UG$a=3]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);逻辑语句选择顺序循环语句在R语言中都可以实现，这里简单说一下if选择语句。SD - function(x){ if(length(x)=1) { return(‘error’)} else{ ss = sum((x-mean(x))^2); return(sqrt(ss/length(x))) }}回顾前面所述：UG1 = UG[ 1:2, 4:5] 其实就是一个选择的过程。练习：求分段函数积分?编辑一个求已知分段函数定积分的函数，分段函数如下： ，求；也就是说求不同x时下图的红线左边面积。calculus - function(x){ if(x3 | x=0){return(error)} if(x0 x=1){return(x*1)} if(x1 x=2){return(1*1+2*(x-1))} if(x2 x=3){return(1*1+2*1+3*(x-2))} }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=merged ASPre)lines(Age,Pos/(Pos+Neg))图像导出与编辑作图完成后一般都需要放到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, t