【正文】 ? 在 R中充分使用向量化，因?yàn)?R會(huì)立即調(diào)用 C進(jìn)行運(yùn)算，因而大大提高計(jì)算的效率！ ? 例子 –見 KoKang Wang’s ―R Programming Workshop‖, pp911 (Gamma函數(shù)作圖 ) ? 沒有使用向量和循環(huán)的源程序 ? 使用 for循環(huán)改進(jìn)后的源程序 ? 使用向量化后源程序 R語言與統(tǒng)計(jì)分析 – 上海師范大學(xué) (2022年 2月 ) 58 ? 函數(shù) ? 函數(shù)是一系列語句的組合，在 R中可以寫出自己的函數(shù) ? 形式 : 變量名 = function( 變量列表 ) 函數(shù)體 ? 函數(shù)引用 : 變量名 (變量的值 ) ? 函數(shù)可以遞歸引用，但不提倡！ ? 例子 – 使用 gamma函數(shù)求 n! factorial = function(n) { + if (n=0) gamma(n+1) + else print(―Please input a positive integer!) + } factorial(6) factorial(6) R語言與統(tǒng)計(jì)分析 – 上海師范大學(xué) (2022年 2月 ) 59 ? 用于處理錯(cuò)誤的函數(shù) – 用于處理用戶輸入不正確的類型而可能出現(xiàn)的錯(cuò)誤 ? warning() 若錯(cuò)誤不嚴(yán)重以至影響整個(gè)計(jì)算 ? stop() 若錯(cuò)誤可能導(dǎo)致計(jì)算中止 ? print() 顯示必要的信息 ? formatC() 數(shù)值作為字符串輸出 ? cat() 字符串聯(lián)，可以插入 \n(換行 )及 \t(tab鍵 ) ? paste() 字符粘貼 (非字符型自動(dòng)轉(zhuǎn)換 ) 例子： cat(R, is,a good, software.\n) formatC(1/3, format = f, digits = 4) formatC(1/3, format = e, digits = 4) paste(1:12) 與 (1:12)等價(jià) paste(A, 1:6, sep = ) paste(today is, date()) R語言與統(tǒng)計(jì)分析 – 上海師范大學(xué) (2022年 2月 ) 60 R中的編程 ? 編程的重要性 ? 一個(gè)統(tǒng)計(jì)軟件 (包 )應(yīng)該包羅萬象 —所有統(tǒng)計(jì)方法或技術(shù)，可以完成所有的任務(wù)？？ ? 這是一種錯(cuò)誤的觀點(diǎn)！！ ? 太貴??！ ? 費(fèi)時(shí)?。? ? SAS是一個(gè)世界上最為優(yōu)秀的統(tǒng)計(jì)分析軟件之一，但 …… ！！當(dāng)然 SAS具體編程能力，而 SPSS更糟！ ? 通過編程可以實(shí)現(xiàn) (開發(fā) )滿足自己需要的函數(shù)或宏包 R語言與統(tǒng)計(jì)分析 – 上海師范大學(xué) (2022年 2月 ) 61 ? 好的編程習(xí)慣 ? 為了他人，更為你本人 ! 你的程序應(yīng)該具有 ? 可讀性 (readability) ? 可理解性 (understandability) ? 習(xí)慣之一：行前縮進(jìn) (Indentation),在此推薦使用軟件 ? WinEdt = 已開發(fā)了基于 WinEdt的 R使用平臺(tái) ? (X)Emacs =已開發(fā)了基于 (X)Emacs的 R使用平臺(tái) ? UltraEdit ? 習(xí)慣之二：增加注釋 (Commenting),它是你的幫助 R中使用 作為注釋語句的開始 . ? 習(xí)慣之三：變量的命名 ,使用意義明確的名字，切忌使用人或?qū)櫸锏拿? R語言與統(tǒng)計(jì)分析 – 上海師范大學(xué) (2022年 2月 ) 62 例子 – 縮進(jìn) for(i in 1:1522) { + aov(Intensity ~ A + T + A * T, + sub = G == i) + [i, ] residuals() + [i, ] coef() + } 或 for(i in 1:1522) { + aov(Intensity ~ A + T + A * T, + sub = G == i) + [i, ] residuals() + [i, ] coef() + } R語言與統(tǒng)計(jì)分析 – 上海師范大學(xué) (2022年 2月 ) 63 ? 實(shí)例演示 ? 1． Julian Faraway提供的一個(gè)簡短的入門 Introduction to R ? 使用 R GUI ? 使用 R Commander ? 使用 R_XEmacs Further reference: Moore39。 ? 建立數(shù)據(jù)框架的方法 (函數(shù) ) ? 若你的數(shù)據(jù)本身保存在一個(gè)文件中，則可以使用 ? () 僅接受帶有分界符的 ASCII數(shù)據(jù) 如果數(shù)據(jù)是電子報(bào)表的形式，則采用下面的兩種變型 ? () 先將數(shù)據(jù)另存為帶逗號(hào)的數(shù)據(jù) (Comma Seperated values) ? () 先將數(shù)據(jù)另存為用 tab作為分界符的數(shù)據(jù) 注：若數(shù)據(jù)集很大 (如 1,000,000觀測 x200變量 )，則可由 ODBC聯(lián)接由數(shù)據(jù)庫讀入 . R語言與統(tǒng)計(jì)分析 – 上海師范大學(xué) (2022年 2月 ) 49 ? 若你在 R中建立了一些向量并試圖想由它們生成框架，則可以使用 ()，但需要同時(shí) cbind(). 例子： x=c(42,7,64,9) y=1:4 =(cbind(INDEX = y, VALUE = x)) INDEX VALUE 1 1 42 2 2 7 3 3 64 4 4 9 注 ? .df 只是為了提醒自己 ? INDEX和 VALUE是重新命名的向量名字 R語言與統(tǒng)計(jì)分析 – 上海師范大學(xué) (2022年 2月 ) 50 ? 數(shù)據(jù)框架子集的提出取 ? 格式 foo[row, column] 其中 foo 數(shù)據(jù)框架的名稱 row 需要提出取的行號(hào) column 需要提出取的行號(hào) ? 例子 (續(xù)上一例 ) [1,] 提取第一個(gè)觀測 (第一行 ) INDEX VALUE 1 1 42 [,1] 提取第一列 (第一個(gè)變量的值 ) [1] 1 2 3 4 [2,1] 提取第二行第一列的元素 [1] 2 R語言與統(tǒng)計(jì)分析 – 上海師范大學(xué) (2022年 2月 ) 51 ? 列表 (list) ? 復(fù)雜的數(shù)據(jù)分析時(shí)，僅有向量與數(shù)據(jù)框架還不夠 ? 有時(shí)需要生成包含不同類型的對象 ? R的列表 (list)就是包含任何類型的對象 例子 : foo = list(x = 1:6, y = matrix(1:4, nrow = 2)) foo $x [1] 1 2 3 4 5 6 $y [,1] [,2] [1,] 1 3 [2,] 2 4 R語言與統(tǒng)計(jì)分析 – 上海師范大學(xué) (2022年 2月 ) 52 ? 列表子集的提出取 ? 提取一個(gè)子對象如 foo的 x,下面三種方式等價(jià) foo$x foo[1] foo[[1]] 例子 foo$y foo[2] foo[[2]] foo[[1]][2] foo$y[2] foo$y[4] R語言與統(tǒng)計(jì)分析 – 上海師范大學(xué) (2022年 2月 ) 53 ? 條件語句 ? 作用 : 避免除零或負(fù)數(shù)的對數(shù)等數(shù)學(xué)問題 ? 形式 1： if (條件 ) 表達(dá)式 1 else 表達(dá)式 2 ? 形式 2 – 常優(yōu)于形式 1！ ifelse(條件 , yes, no) 試比較下面的三個(gè)結(jié)果： x = c(6:4) sqrt(ifelse(x = 0, x, NA)) ifelse(x = 0, sqrt(x), NA) if (x = 0) sqrt(x) else NA R語言與統(tǒng)計(jì)分析 – 上海師范大學(xué) (2022年 2月 ) 54 ? 循環(huán) (loops) ? for() 若知道終止條件 for (變量 in 向量 ) 表達(dá)式 ? while() 若無法知道運(yùn)行次數(shù) while(條件 ) 表達(dá)式 ? 兩者通?？梢赞D(zhuǎn)換 例 1—試比較兩種方法 for (i in 1:5) print (1:i) i=1 while(i = 5) { + print(1:i) + i = i+1 + } R語言與統(tǒng)計(jì)分析 – 上海師范大學(xué) (2022年 2月 ) 55 例 2 – 見 KoKang Wang’s ―R Programming Workshop‖, pp68 ? Suppose we ge