軟件學習語言初步

—1報告—R簡二函數(shù)與對象 四R繪圖寫函數(shù)六數(shù)據(jù)保2—R簡3什么是RTheRProjectforStatistical圖1R4簡史R語言是從S統(tǒng)計繪圖語言演變而來，可看作S的“方言 ，由RickBecker,JohnChambers,Wilks開發(fā)1995年由新西蘭Auckland大學統(tǒng)計系的RobertGentleman和RossIhaka，基于S5R的特目前在 6為什么選擇GNU軟7R與其他統(tǒng)計軟件比 運行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文件，安裝時不要解壓縮安裝程序用函數(shù)序?qū)⒆?picante包,路徑：Packages>installpackagesfromlocal 程序包使Thisisvegan1.17-查看程序包幫助文vegan查看pdf查看函數(shù)的幫助文函數(shù)的默認值是什么？怎么使用？使用時需要注意什需要查詢函數(shù)的幫RGui>Help>Html查看R包pdf手11R幫助文件的內(nèi)Fitting 安裝R并導入程序安裝R輸入?plot.phyloR圖形界面：RRcommander是R的圖形界面之一，是JohnFox教授編寫的，適用于不希望安裝RR將自 并安裝Rcmdr所需的所有程序圖形界面之一：R界面操作代結(jié)果輸圖12Rcommander界面 為什么要學習編程 要學習R編程，首先要了解R數(shù)與對象R的函函數(shù)(輸入數(shù)據(jù)，參數(shù)平均值mean(x,trim0na.rmFALSE,線性模型lm(y~xR的函 R有哪些函數(shù)查詢的方

glm(y~x1+x2+x3,binomial)R函數(shù)調(diào)用及其選圖13箱線圖修飾前后(左：默認值，右：修改屬性后函數(shù)的調(diào)用方法函數(shù)名如plotlm(),R函數(shù)調(diào)用及其選boxplot(day~type,data=bac,col="red",xlab="

練習二：查詢函數(shù)幫查看最后的 boxplot(count~spray,data=InsectSprays,col參數(shù)更boxplot(count~spray,data=InsectSprays,col=xlab="spray",賦值與注2+a<-<-也可用甚至->b<-c<-a+b如何為區(qū)分大小寫,注 尾，如2result與result2，后者是合法作為間隔，例如NA,NaN,pi,LETTERS,letters,month.abb,元素的類數(shù)值型 如100,0,-字符型Character如“ 邏輯型Logical 如TRUE,FALSE因子型Factor 復數(shù)型 如：2+對象的類向量(vector如c(1,2,3因子(factorx<-1:12;dim(x)<-147258369對象數(shù)組數(shù)組是k維的數(shù)據(jù)表（kin1:n,n為正整數(shù)）。向量(n=1)矩陣(n=2) 數(shù)組(n>=3)數(shù)據(jù)框 運算 算運算后給出數(shù)值結(jié)*^冪(TRUE>,<,<=,>=,==,!,&,&&,|,數(shù)據(jù)表的行與表1

樣方樣方樣方 行名Rownames每列可看做帶名稱

數(shù)據(jù)框(1)從外部數(shù)外部數(shù)最為常用的數(shù) Excel,SAS,spsssystatMinitabread.csv()的使read.table(file=file.choose(),header=T)可以彈出 例：從數(shù)據(jù)輸入到t檢現(xiàn)有是否等于22.。123456 式：從控制臺輸入數(shù)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)第二種方式從外 數(shù)第1步將Excel中的數(shù)據(jù)另存為.txt格式（制表符間隔）或.csv圖14在Excel中將數(shù)據(jù)存為txt例：t檢驗(續(xù) bmi<-t.test(bmimu=22.5)#t練習三：數(shù) 和t檢 t.test.data<-read.table(“X:/t.test.txt”,例：單因素方差分 a,b,c分別接種于10，9，和11只小白鼠上，a菌株：243247725b菌株：5685107126c菌株：7,11,6,6,7951063準備數(shù)據(jù)圖15day和type例：方差分析(續(xù) ba.an<-aov(lm(day~type,data=bac))圖16 向量、矩陣和數(shù)據(jù)框的創(chuàng)向量的創(chuàng) ","Korea","Japan","UK","France", ",復數(shù)型略向量的創(chuàng)生成向量的函數(shù)seq(from=3,to=21,by=3[1]369121518[1]1234567891011121314向量創(chuàng)建——產(chǎn)生隨機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,矩陣的創(chuàng)dim()和matrix()dim()定義矩陣的行列數(shù)，例如：x<-dim(x)<-147258369 數(shù)據(jù)框cbind()#按列組合成數(shù)據(jù)框rbind()#按行組合成數(shù)據(jù)框data.frame()#生成數(shù)據(jù)框head()#默認 列表的創(chuàng)類的判mode()判 classis.logical()#是否為邏輯值is.charactor()#是否為字符串is.null()#是否為空is.na#是否為類的as.numeric()#轉(zhuǎn)換為數(shù)值型as.logical()#轉(zhuǎn)換為邏輯型as.matrix()#轉(zhuǎn)換為矩陣as.factor()#轉(zhuǎn)換為因子向量內(nèi)的### 例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);

數(shù)據(jù)框內(nèi)元素intake<-data.frame(intake.pre,數(shù)據(jù)框中的 例如 i=1:5;

列表內(nèi)元列表內(nèi)元素 可以用如上述放于”[]”內(nèi)的數(shù)字，稱為下標。通過下標的變化，可以方便的向量、數(shù)據(jù)框內(nèi)的元 transform#數(shù)據(jù)的轉(zhuǎn)換,如取subset#條件intake.pre>[1]FALSEFALSEFALSEFALSETRUETRUEintake.pre[intake.pre>

FALSEFALSEintake.pre>7000&intake.pre<=[1]FALSEFALSEFALSEFALSEFALSEFALSEFALSETRUEFALSEintake.post[intake.pre>7000&intake.pre<=練習五：下標和條件創(chuàng)建一個2到50的向量2,4,6,8,...,48,vector1<-seq(from=2,to=50, 選取vector1中的第10，15，20個元排序ordero<-工作空間lsrm()刪除工作空間中的對象save.image()保存工作鏡像sink將運行結(jié)果保存到指定文件中g(shù)etwd()顯示當前工作文件夾setwd()設定工作文件夾輸入歷輸 history(Inf查看打開R ：了解工作路1查看當前R 3 編語 常見 Windows批處理程序PHPPerlPythonRubyJavaScript集成開發(fā)環(huán)很多計算機語言都有IDE(integrateddevelopmentenvironment)境,如Windows的VisualStudio,VisualBasic等等但是R通常無需集成開發(fā)環(huán)境 ++,Linux下的Vi,VimEmacs 編輯 Editplus TinnR)Ultraedit）Emacs）（）記事本或?qū)懽职宓鹊葓D17TinnR對 圖18NotePad++對 顯例-線性回表3攪拌速度對涂料中雜質(zhì)的影轉(zhuǎn)%舉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軟件的結(jié)合，可以生成600dpi，1200dpi的等各種分辨率繪圖函數(shù) 高級繪

圖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繪制的圖形。五編寫函編程函數(shù)名function(數(shù)據(jù),參數(shù)1默認值{}圖28函數(shù)實例

程序流程控if(條件if(條件表達式1p={

}

print("p<=print("p>循環(huán)forfor(變量in向量for(iin1:10)while(條件i<-{i<-i+}

返回異常處inputdatana.omit(inputdata)cat(“NAsarefoundintheinputdata,andhasstop("NAsarenot函數(shù)rcal<-{z<-x^2+y^2result<-sqrt(z)}練習九編寫函{}數(shù)據(jù)Sweave該函數(shù)能 主要講了什么 R編寫函數(shù)無 推薦CrawleyStatisticsanintroductionusingPeterDalgaardIntroduc