R語(yǔ) 1RR1980S的一個(gè)分支?？梢哉J(rèn)為R是S語(yǔ)言的一種實(shí)現(xiàn)。而S語(yǔ)言是由AT&T貝爾 據(jù)探索、統(tǒng)計(jì)分析和作圖的解釋型語(yǔ)言。最初S語(yǔ)言的實(shí)現(xiàn)版本主要是S-PLUS。AucklandRobertGentleman和RossIhaka及其他志愿人員開發(fā)了一個(gè)R系統(tǒng)。手冊(cè)。所以有人說(shuō)：R，是S-PLUS的一個(gè)“克隆”。 了：R是免費(fèi)的(R )GNU SplusC1/20。R語(yǔ)言 R1.5R語(yǔ)言 TheComprehensiveRArchiveNetwork簡(jiǎn)稱CRAN，由世界幾十個(gè)鏡像 組成網(wǎng)絡(luò)，提供 安裝程序和相應(yīng)軟件包。各鏡像更新頻率一般為1-2天推薦鏡像：中國(guó)的鏡像:數(shù)學(xué) 即時(shí)更新的CRAN1.5R語(yǔ)言 Windows CRAN:RBinaries>windows>baseR-3.0.1-win.exe1.5R語(yǔ)言 Mac CRAN:RBinaries>macosxR-3.0.1.pkg光標(biāo)等待輸菜單欄光標(biāo)等待輸控制臺(tái)R(Windows版路徑：開始>RRGuiFileRGui的EditRGuiPackagesRGuiHelp程序包是什么？Window下的Rzip 為什么要安裝程序包？apevegan包等。常用R

Bootstrap檢驗(yàn)常用R

常用R程序包(

DNA序列分析在CRAN提供了每個(gè)包的源代碼和編譯好的程序包以vegan包為例，CRANPackagesource:vegan_2.0-7.tar.gzMacOSXbinary:vegan_2.0-7.tgzWindowsbinary:vegan_2.0-7.zipReferencemanual:vegan.pdfWindow下程序包為zipMactgz文件，安裝時(shí)不要解壓縮。用函數(shù)install.packages()， 例如：要安裝veganinstall.packages("安裝本地zipPackages>installpackagesfromlocalfiles選擇本地磁盤 zip包的文件夾library(vegan)library(ape)library(vegan)Thisisvegan2.0-vegan程序包內(nèi)部都有哪些函數(shù)？分別Html幫助查看pdf練R3函數(shù)形式:函數(shù)(輸入數(shù)據(jù)，參數(shù)=)例如:平均值meanx,trim0,na.rmFALSE,...)線性模型lm(y~x, 級(jí)函數(shù)可調(diào)用低級(jí)函數(shù),這里的”高級(jí)”函數(shù)習(xí)慣上稱為泛型函數(shù)。如plot()就是泛型函數(shù)，可以根據(jù)數(shù)據(jù)的類型，調(diào)用底層的函數(shù)，應(yīng)R有哪些函數(shù)?查詢的方法：Help>Htmlhelp>packages

glm(y~x1+x2+x3,binomi函數(shù)的調(diào)用方法,函數(shù)名()如plot(),lm(),3.1.2Rboxplot(day~type,data=bac,col="red",xlab="Virus",ylab="days")day~type，以typedaycol=“red”

數(shù)據(jù)來(lái)源xlab=“ ylab=“days”

橫軸名稱為縱軸名稱為區(qū)分大小寫,注 或結(jié)尾，如2result與result2后者是合法的。1NA,NaN,pi,LETTERS,letters,month.abb, 數(shù)值型Numeric字符型Character邏輯型Logical

如100,0,-4.335如 TRUE,因子型 復(fù)數(shù)型 如：2+對(duì)象的類(class)向量(vector)一系列元素的組合。如c123);caa","bbc")因子(factor)c("a","a","b","b","矩陣(matrix)x<-1:12 dim(x)<-c(3,4)[,1][,2][,3][,147258369對(duì)象的類(class)數(shù)組(array)數(shù)組是k維的數(shù)據(jù)表（kin1:n,n 向量(n=1) 矩陣(n=2) 數(shù)組(n>=3)數(shù)據(jù)框(dataframe)列表( 算運(yùn)算后給出數(shù)值結(jié)+,-,*,/,^比較運(yùn)算運(yùn)算后給出判別結(jié)果(TRUE>,<,<=,>=,==,!,&,&&,|,|4015226002551122635030461120390453812242602049102522033

Entry 行名Row

每個(gè)數(shù)據(jù)表可以看作一個(gè)數(shù)據(jù)框(dataframe)(column)作為一個(gè)向量(vector)(row)作為一個(gè)記錄(entry) 方式是用read.table()函數(shù) 一些R程序包（如foreign）也提供了直接 el,SAS,dbf, ,spss,systat,Minitab文件的函數(shù)。 ","Korea","Japan",K","USA","France", ",numeric<-c(1,3,6,7,3,8,6,logical<-c(T,F,T,F,T,F,F,生成向量的函數(shù)c(),rep(),seq(),”:”seq(from=3,to=21,by=3[1]36912151821 4生成矩陣的函數(shù)dimmatrix()dim()定義矩陣的行列數(shù)，例如：x<-147258369matrix.xmatrix(1:12,nrow=3,byrow=T)t(x轉(zhuǎn)置:data.frame(.frame(),cbind(),##data.frame is.numeric()#返回值為TRUE或FALSEis.logical()#是否為邏輯值is.charactor()is.null()is.na()#是否為naas.numeric()#轉(zhuǎn)換為數(shù)值型as.logical()#轉(zhuǎn)換為邏輯型as.charactor()#轉(zhuǎn)換為字符串a(chǎn)s.matrix()#轉(zhuǎn)換為矩陣as.data.frame()#轉(zhuǎn)換為數(shù)據(jù)框as.factor()#轉(zhuǎn)換為因子intake.pre<-c(5260,5470,5640,6180,6390,15,6805,7515,7515,8230,

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

第1到5去除第357intake<-data.frame(intake.pre, intake.post)數(shù)據(jù)框中的元素 例如intake$intake.pre[,]方括 intake[,1];intake[5,];intake[5,1];i=1:5;intake[i,]

1列5行151到5行3.3.4數(shù)據(jù) 可以用”[[]]” list1[[c(1,2,3)]]上述放于”[]”內(nèi)的數(shù)字，稱為下標(biāo)。通過(guò) 3.3.4數(shù)據(jù)attatch()with()#with函數(shù)內(nèi)部形成一個(gè)空間，在這個(gè)空間中，函數(shù)可以 within()transform()#數(shù)據(jù)的轉(zhuǎn)換,如取logsubset()apply()3.3.5條件篩選是先對(duì)變量否滿足條件進(jìn)行判斷，滿足為TRUE，不滿足為FALSEintake.pre>7000[1]FALSEFALSEFALSEFALSEFALSEFALSEFALSE TRUE TRUE TRUE TRUEintake.pre[intake.pre>7000]intake.pre>7000&intake.pre<=8000[1]FALSEFALSEFALSEFALSEFALSEFALSEFALSE TRUE TRUEFALSEFALSEintake.post[intake.pre>7000&intake.pre<=8000]3.3.5 intake$intake.presort(intake$intake.pre)order() order(intake$post)o<-intake[order(intake$intake.pre),]4.R jpgtiffepsemfpdfpngGhostScript600dpi1200dpi

font=字體lty=線類型lwd=pch=點(diǎn)的類型,xlab=橫坐標(biāo)ylab=縱坐標(biāo)xlim=ylim=參見par()4.4舉例：繪圖xy橫坐標(biāo)”，縱軸命名為“縱坐標(biāo)”plot(xymain散點(diǎn)圖xlab橫坐標(biāo)",ylab縱坐標(biāo)")text(0.6,0.6,"textat(0.6,0.6)")plot(x,y,type="n",xlab="",ylab="",axis(1);axis(at=seq(0.2,1.8,0.2),title(main="Maintitle",sub="subtitle",xlab="x-label",ylab="y-label")11