软件学习语言初步

—1报告—R简二函数与对象 四R绘图写函数六数据保2—R简3什么是RTheRProjectforStatistical图1R4简史R语言是从S统计绘图语言演变而来，可看作S的“方言 ，由RickBecker,JohnChambers,Wilks开发1995年由新西兰Auckland大学统计系的RobertGentleman和RossIhaka，基于S5R的特目前在 6为什么选择GNU软7R与其他统计软件比 运行S语言，具有复杂的界面，与R完全8R的缺即时编译，约相当于C语言的1/209图2RTheComprehensiveRArchive简称CRAN，由世界几十个镜 安装程序和相应软件包 和安装CRAN:图3 完成后，双击R-2.11.0-win32.exe光标:等光标:等待输快捷按控制4R登陆界面(Windows版图5RGui的File菜 图6RGui的Edit菜 图7RGui的Packages菜 图8RGui的Help菜 R程序包（R程序每个程 为什么要安装程序

常用R程序包

常用R程序包空间对象 和处

常用R程序包CRANTask10veganR程序Packagesource:vegan_1.17-2.tar.gzMacOSXbinary:vegan_1.17-2.tgzWindowsbinary:vegan_1.17-2.zipReferencemanual:vegan.pdfWindow下程序包为zip文件，安装时不要解压缩安装程序用函数序将自 picante包,路径：Packages>installpackagesfromlocal 程序包使Thisisvegan1.17-查看程序包帮助文vegan查看pdf查看函数的帮助文函数的默认值是什么？怎么使用？使用时需要注意什需要查询函数的帮RGui>Help>Html查看R包pdf手11R帮助文件的内Fitting 安装R并导入程序安装R输入?plot.phyloR图形界面：RRcommander是R的图形界面之一，是JohnFox教授编写的，适用于不希望安装RR将自 并安装Rcmdr所需的所有程序图形界面之一：R界面操作代结果输图12Rcommander界面 为什么要学习编程 要学习R编程，首先要了解R数与对象R的函函数(输入数据，参数平均值mean(x,trim0na.rmFALSE,线性模型lm(y~xR的函 R有哪些函数查询的方

glm(y~x1+x2+x3,binomial)R函数调用及其选图13箱线图修饰前后(左：默认值，右：修改属性后函数的调用方法函数名如plotlm(),R函数调用及其选boxplot(day~type,data=bac,col="red",xlab="

练习二：查询函数帮查看最后的 boxplot(count~spray,data=InsectSprays,col参数更boxplot(count~spray,data=InsectSprays,col=xlab="spray",赋值与注2+a<-<-也可用甚至->b<-c<-a+b如何为区分大小写,注 尾，如2result与result2，后者是合法作为间隔，例如NA,NaN,pi,LETTERS,letters,month.abb,元素的类数值型 如100,0,-字符型Character如“ 逻辑型Logical 如TRUE,FALSE因子型Factor 复数型 如：2+对象的类向量(vector如c(1,2,3因子(factorx<-1:12;dim(x)<-147258369对象数组数组是k维的数据表（kin1:n,n为正整数）。向量(n=1)矩阵(n=2) 数组(n>=3)数据框 运算 算运算后给出数值结*^幂(TRUE>,<,<=,>=,==,!,&,&&,|,数据表的行与表1

样方样方样方 行名Rownames每列可看做带名称

数据框(1)从外部数外部数最为常用的数 Excel,SAS,spsssystatMinitabread.csv()的使read.table(file=file.choose(),header=T)可以弹出 例：从数据输入到t检现有是否等于22.。123456 式：从控制台输入数height<-c(1.75,1.80,1.65,1.90,1.74,weight<-c(60,72,57,90,95,t.test(ratio,mu=22.5)第二种方式从外 数第1步将Excel中的数据另存为.txt格式（制表符间隔）或.csv图14在Excel中将数据存为txt例：t检验(续 bmi<-t.test(bmimu=22.5)#t练习三：数 和t检 t.test.data<-read.table(“X:/t.test.txt”,例：单因素方差分 a,b,c分别接种于10，9，和11只小白鼠上，a菌株：243247725b菌株：5685107126c菌株：7,11,6,6,7951063准备数据图15day和type例：方差分析(续 ba.an<-aov(lm(day~type,data=bac))图16 向量、矩阵和数据框的创向量的创 ","Korea","Japan","UK","France", ",复数型略向量的创生成向量的函数seq(from=3,to=21,by=3[1]369121518[1]1234567891011121314向量创建——产生随机runif(10,min=0,max=rnorm(10,mean=0,sd=dunif(x,min=0,max=1,log=punif(q,min=0,max=1,qunif(p,min=0,max=1,runif(n,min=0,矩阵的创dim()和matrix()dim()定义矩阵的行列数，例如：x<-dim(x)<-147258369 数据框cbind()#按列组合成数据框rbind()#按行组合成数据框data.frame()#生成数据框head()#默认 列表的创类的判mode()判 classis.logical()#是否为逻辑值is.charactor()#是否为字符串is.null()#是否为空is.na#是否为类的as.numeric()#转换为数值型as.logical()#转换为逻辑型as.matrix()#转换为矩阵as.factor()#转换为因子向量内的### 例intake.pre<-c(5260,5470,5640,6180,6390,6515,7515,7515,8230,intake.post<-c(3910,4220,3885,5160,5645,4680,5975,6790,6900,

v<-c(3,5,7);

数据框内元素intake<-data.frame(intake.pre,数据框中的 例如 i=1:5;

列表内元列表内元素 可以用如上述放于”[]”内的数字，称为下标。通过下标的变化，可以方便的向量、数据框内的元 transform#数据的转换,如取subset#条件intake.pre>[1]FALSEFALSEFALSEFALSETRUETRUEintake.pre[intake.pre>

FALSEFALSEintake.pre>7000&intake.pre<=[1]FALSEFALSEFALSEFALSEFALSEFALSEFALSETRUEFALSEintake.post[intake.pre>7000&intake.pre<=练习五：下标和条件创建一个2到50的向量2,4,6,8,...,48,vector1<-seq(from=2,to=50, 选取vector1中的第10，15，20个元排序ordero<-工作空间lsrm()删除工作空间中的对象save.image()保存工作镜像sink将运行结果保存到指定文件中getwd()显示当前工作文件夹setwd()设定工作文件夹输入历输 history(Inf查看打开R ：了解工作路1查看当前R 3 编语 常见 Windows批处理程序PHPPerlPythonRubyJavaScript集成开发环很多计算机语言都有IDE(integrateddevelopmentenvironment)境,如Windows的VisualStudio,VisualBasic等等但是R通常无需集成开发环境 ++,Linux下的Vi,VimEmacs 编辑 Editplus TinnR)Ultraedit）Emacs）（）记事本或写字板等等图17TinnR对 图18NotePad++对 显例-线性回表3搅拌速度对涂料中杂质的影转%举rate<-c(20,22,24,26,28,30,32,34,36,38,40,impurity<-c(8.4,9.5,11.8,10.4,13.3,14.8,13.2,16.4,16.5,18.9,运通过 路径：RGui>File>Open #Ctrl+R运ctrl+c,ctrl+v练习七 运 R绘图19R绘制图20左图maps包map()右图 包绘制图21fields包实 spdep包实空间分图22gstatR绘图功 通过与GhostScript软件的结合，可以生成600dpi，1200dpi的等各种分辨率绘图函数 高级绘

图23

低级绘

绘图font=字体lty=线类型lwd=线宽度pch点的类型,xlabylabxlim

ylim参见

举例：绘x<-y<-plot(x,ymain="散点图xlab="横坐标ylab="纵坐标text(0.6,0.6,"textat图24

例：分步plot(x,y,type="n",xlab="",ylab="",123456图25

一页图26一页多图在原有图形上添加元x<-###h<-hist(x,#ylimrange(0,h$densitydnorm(0#设定纵轴的取值范围hist(x,freq=F,ylim=ylim)#绘制直方图curve(dnorm(x),add=T,col="red")#添加曲线图27练习八练习1将Rplots.r中的代码拷贝到R控制台中，查看R绘制的图形。五编写函编程函数名function(数据,参数1默认值{}图28函数实例

程序流程控if(条件if(条件表达式1p={

}

print("p<=print("p>循环forfor(变量in向量for(iin1:10)while(条件i<-{i<-i+}

返回异常处inputdatana.omit(inputdata)cat(“NAsarefoundintheinputdata,andhasstop("NAsarenot函数rcal<-{z<-x^2+y^2result<-sqrt(z)}练习九编写函{}数据Sweave该函数能 主要讲了什么 R编写函数无 推荐CrawleyStatisticsanintroductionusingPeterDalgaardIntroduc