【正文】 Ripley, B. D. Modern Applied Statistics with S 123 R網(wǎng)絡(luò)資源 R主頁(yè) : R資源列表 NCEAS R Graphical Manual 統(tǒng)計(jì)之都 : QuikR 丁國(guó)徽的 R文檔 : R語(yǔ)言中文論壇 124 謝謝 敬請(qǐng)指正 ! 。 4. R編寫(xiě)函數(shù)無(wú)需聲明變量的類(lèi)型，能利用循環(huán)、條件語(yǔ)句，控制程序的流程。 2. R中的向量、列表、數(shù)組、函數(shù)等都是對(duì)象，可以方便的查詢(xún)和引用，并進(jìn)行條件篩選。 } 119 六 數(shù)據(jù)保存 120 () () () sink() unlink() 若有 LaTeX基礎(chǔ)，可以用 Sweave() 函數(shù) 該函數(shù)能將腳本、程序說(shuō)明和運(yùn)算結(jié)果直接保存成 .tex文件，用LaTeX編譯成 pdf文件。 z3=z1+z2。 sqtestfunction(x, y){ z1=x^2。 ? 警告的寫(xiě)法 if(any((inputdata))) inputdata (inputdata) cat(“NAs are found in the input data, and has been removed.\n) ? 終止的寫(xiě)法 if(any((xx))) stop(NAs are not allowed!\n) 117 函數(shù)舉例 問(wèn)題 :輸入直角三角形的兩個(gè)邊長(zhǎng)，求其斜邊長(zhǎng)。 ? 但是一個(gè)函數(shù)的返回的對(duì)象只有一個(gè)。 ? 如果函數(shù)的結(jié)果需要有多個(gè)返回值，可以創(chuàng)建一個(gè) list()，并返回該對(duì)象。 ? 返回值必須是一個(gè)對(duì)象。 } 函數(shù)內(nèi)部也可用 添加注釋 112 函數(shù)實(shí)例 ? 圖 28 函數(shù)實(shí)例 data2mat() 113 程序流程控制 if if(條件 ) 表達(dá)式 if(條件 ) 表達(dá)式 1 else 表達(dá)式 2 舉例 p = { if(p=) print(p = !) else print(p !) } 114 循環(huán) for, while for(變量 in 向量 ) 表達(dá)式 用法： for(i in 1:10) print(i) while(條件 ) 表達(dá)式 用法 : i 1 while(i10){ print(i)。 表達(dá)式 (循環(huán) /判別 )。編程時(shí)無(wú)需聲明變量的類(lèi)型，這與 C,C++等語(yǔ)言不同。 練習(xí) 2 對(duì)例進(jìn)行回歸分析，并繪制散點(diǎn)圖，并為散點(diǎn)圖添加回歸直線。 font = 字體 ,lty = 線類(lèi)型 ,lwd = 線寬度 ,pch = 點(diǎn)的類(lèi)型 , xlab = 橫坐標(biāo) ,ylab = 縱坐標(biāo) , xlim = 橫坐標(biāo)范圍 , ylim = 縱坐標(biāo)范圍 , 也可以對(duì)整個(gè)要繪制圖形的各種參數(shù)進(jìn)行設(shè)定 參見(jiàn) par() 102 舉例：繪圖 生成 0到 2之間的 50個(gè)隨機(jī)數(shù)，分別命名為 x,y x runif(50,0,2) y runif(50,0,2) 繪圖：將主標(biāo)題命名為“散點(diǎn)圖”， 橫軸命名為”橫坐標(biāo)”， 縱軸命名為“縱坐標(biāo)” plot(x, y, main=散點(diǎn)圖 , xlab=橫坐標(biāo) , ylab=縱坐標(biāo) ) text(,text at (,)) abline(h=.6,v=.6) 103 圖 24 繪圖舉例 104 例：分步繪圖 ，不繪制任何對(duì)象 plot(x, y, type=n, xlab=, ylab=, axes=F) points(x,y) axis(1)。 繪圖函數(shù)包括 高級(jí)繪圖函數(shù) 和 低級(jí)繪圖函數(shù) 。 通過(guò)與 GhostScript軟件的結(jié)合，可以生成 600dpi， 1200dpi的等各種分辨率和尺寸的圖形。 92 四 R繪圖 93 圖 19 R繪制的圖形 94 繪制地圖 ? 圖 20 左圖 maps包 map()右圖 PBSmapping包 addBubbles() 95 繪制地圖 ? 圖 21 fields 包實(shí)例 spdep 包實(shí)例 96 空間分析繪圖 圖 22 gstat程序包實(shí)例 97 R繪圖功能 R具備卓越的繪圖功能，通過(guò)參數(shù)設(shè)置對(duì)圖形進(jìn)行精確控制。 ratec(20, 22, 24, 26, 28, 30, 32, 34, 36, 38, 40, 42) impurity c(, , , , , , , , , , , ) plot(impurity~rate) reglm(impurity~rate) abline(reg,col=red) summary(reg) 90 運(yùn)行腳本 三種運(yùn)行方式 1 通過(guò) source()函數(shù)運(yùn)行 source(d:/) 2 通過(guò) R腳本編輯器運(yùn)行 路徑： RGuiFileOpen Script Ctrl+R運(yùn)行 3 直接粘貼到 R控制臺(tái) ctrl+c, ctrl+v 第三種最為簡(jiǎn)單直接 91 練習(xí)十 :R腳本運(yùn)行 將 R命令粘貼到記事本中，另存為 。 在代碼較多的情況下，常需要對(duì)行數(shù)、函數(shù)、括號(hào)、函數(shù)選項(xiàng)等進(jìn)行高亮顯示，設(shè)置成不同的顏色，以減少錯(cuò)誤的發(fā)生。 如 Windows自帶的記事本 ,Notepad++, UltraEdit， TinnR。 84 集成開(kāi)發(fā)環(huán)境 很多計(jì)算機(jī)語(yǔ)言都有 IDE (integrated development environment)，即集成開(kāi)發(fā)環(huán)境 ,如 Windows的 VisualStudio, Visual Basic等等。 特點(diǎn) 腳本語(yǔ)言語(yǔ)法和結(jié)構(gòu)通常比較簡(jiǎn)單，不需要編譯，通過(guò)解釋器對(duì)腳本進(jìn)行解釋?zhuān)瑥亩o出結(jié)果，能用簡(jiǎn)單的代碼完成復(fù)雜的功能，但是速度較慢。 80 輸入歷史 history ? 輸入的命令，在 R中作為歷史 history保存，可輸入函數(shù) ? history() 查看輸入歷史 ? history(Inf) 查看打開(kāi) R之后所有的輸入 ? 可用向上或向下的箭頭切換輸入的行。 order(intake$post) o intake[order(intake$),] 79 工作空間 image R的所有對(duì)象都在計(jì)算機(jī)內(nèi)存的工作空間中。 = 8000] 77 練習(xí)五：下標(biāo)和條件篩選 創(chuàng)建一個(gè) 2到 50的向量 vector1 2, 4, 6, 8, ..., 48, 50 vector1seq(from=2, to=50, by=2) 選取 vector1中的第 20個(gè)元素 vector1[20] 選取 vector1中的第 10， 15， 20個(gè)元素 vector1[c(10,15,20)] 選取 vector1中的第 10到 20個(gè)元素 vector1[10:20] 選取 vector1中值大于 40的元素 vector1[vector40] 78 排 序 將向量中的元素按照一定順序排列。 7000 [1] FALSE FALSE FALSE FALSE FALSE FALSE FALSE TRUE TRUE TRUE TRUE [ 7000] 7000 amp。 within() transform() 數(shù)據(jù)的轉(zhuǎn)換 ,如取 log subset() 取數(shù)據(jù)的子集 apply() 對(duì)數(shù)據(jù)表或矩陣應(yīng)用某個(gè)函數(shù)，可減少循環(huán) 76 條件篩選 條件篩選是先對(duì)變量否滿足條件進(jìn)行判斷，滿足為 TRUE，不滿足為FALSE。熟悉下標(biāo)的用法對(duì)掌握循環(huán)結(jié)構(gòu)是非常重要的。 intake[i,] 引用 1到 5行 74 列表內(nèi)元素的引用 列表內(nèi)元素的引用可以用 ” [[]]” 如 list1[[c(1,2,3)]] 上述放于 ” []”內(nèi)的數(shù)字，稱(chēng)為下標(biāo)。 引用第 5行 intake[5,1]。 引用第 1到 5個(gè)元素 [c(3,5,7