《An introduction to R (III)計算機(jī)教學(xué)課件PPT》由會員分享，可在線閱讀，更多相關(guān)《An introduction to R (III)計算機(jī)教學(xué)課件PPT（60頁珍藏版）》請在裝配圖網(wǎng)上搜索。
1、2. An introduction to R (III)2.13 Data input & outputn 鏈接嵌入的數(shù)據(jù)庫鏈接嵌入的數(shù)據(jù)庫 n data() #Data sets in package datasetsn data(USArrests) #創(chuàng)建同名的數(shù)據(jù)框n 若數(shù)據(jù)本身已保存在一個文件中，則可以使用若數(shù)據(jù)本身已保存在一個文件中，則可以使用nread.table() 僅接受帶有分界符的僅接受帶有分界符的ASCII數(shù)據(jù)數(shù)據(jù) 如果數(shù)據(jù)是電子報表的形式，則采用下面的兩種變型如果數(shù)據(jù)是電子報表的形式，則采用下面的兩種變型nread.csv() 先將數(shù)據(jù)另存為帶逗號的文件先將數(shù)據(jù)另存為
2、帶逗號的文件(Comma Seperated Value file) nread.delim() 先將數(shù)據(jù)另存為用先將數(shù)據(jù)另存為用tab作為分界符的文件作為分界符的文件(delimited file) 注：若數(shù)據(jù)集很大注：若數(shù)據(jù)集很大(如如1,000,000觀測觀測x200變量變量)，則可使用，則可使用MySQL讀讀入入. http:/ read.table()n 讀取表格形式純文本文件存為數(shù)據(jù)框n 現(xiàn)有文件houses.datan Price Floor Area Rooms Age Cent.heatn 01 52.00 111.0 830 5 6.2 non 02 54.75 128.0
3、 710 5 7.5 non 03 57.50 101.0 1000 5 4.2 non 04 57.50 131.0 690 6 8.8 non 05 59.75 93.0 900 5 1.9 yesn rt is.data.frame(rt)n 1 TRUEn rt getwd( ) n 使用全路徑使用全路徑n rt mydata cat(TITLE extra line, 2 3 5 7, 11 13 17, file=ex.data, sep=n) n是回車是回車n scan(ex.data, skip = 1)是為了跳過空格是為了跳過空格n #有表頭需要跳過有表頭需要跳過n Read
4、 7 itemsn 1 2 3 5 7 11 13 17n pp scan(ex.data, skip = 1, nlines=1)n Read 4 itemsn 1 2 3 5 7n class(pp)n 1 numeric在當(dāng)前目錄下生成在當(dāng)前目錄下生成ex.data文件文件TITLE extra line 2 3 5 7 11 13 17whatn scan(ex.data, what = list(,)-p;pn Read 4 recordsn 1n 1 TITLE 2 7 17 n 2n 1 extra 3 11 n 3n 1 line 5 13 n 警告信息：警告信息：n In s
5、can(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, :n 讀取的項(xiàng)目數(shù)必需是列數(shù)的倍數(shù)讀取的項(xiàng)目數(shù)必需是列數(shù)的倍數(shù)ex.data內(nèi)容內(nèi)容TITLE extra line 2 3 5 7 11 13 17 scan(ex.data, what = list(,), flush = TRUE) -p; pRead 3 records11 TITLE 2 11 21 extra 3 13 31 line 5 17 class(p)1 list read.table(ex.data,fill=T) TITLE extra l
6、ine2 3 5 711 13 17 NAP45 exam0203.txtn scan(exam0203.txt,what=list(,0,0,0),skip=1)nRead 19 recordsn1n 1 Alice Becka Gail Karen Kathy Mary Sandy n 8 Sharon Tammy Alfred Duke Guido James Jeffreyn15 John Philip Robert Thomas Williamn2n 1 F F F F F F F F F M M M M M M M M Mn19 Mn3n 1 13 13 14 12 12 15 1
7、1 15 14 14 14 15 12 13 12 16 12 11 15n4n 1 56.5 65.3 64.3 56.3 59.8 66.5 51.3 62.5 62.8 69.0 63.5 67.0 57.3 62.5n15 59.0 72.0 64.8 57.5 66.5n5n 1 84.0 98.0 90.0 77.0 84.5 112.0 50.5 112.5 102.5 112.5 102.5 133.0n13 83.0 84.0 99.5 150.0 128.0 85.0 112.0Name Sex Age Height WeightAlice F 13 56.5 84.0Be
8、cka F 13 65.3 98.0Gail F 14 64.3 90.0Karen F 12 56.3 77.0Kathy F 12 59.8 84.5Mary F 15 66.5 112.0Sandy F 11 51.3 50.5Sharon F 15 62.5 112.5n scan還可以用以從鍵盤輸入數(shù)據(jù)還可以用以從鍵盤輸入數(shù)據(jù),格式為格式為:n A An 1 1 2 3 6恢復(fù)輸入形式 n x=scan(f:book1.txt)n 這是按照文本一行一行讀入的數(shù)據(jù)n 如果原先是45矩陣形式,則加用n x=matrix(x,4,5,byrow=T) #或直接用n x= matrix(sc
9、an(f:book1.txt),4,5,b=T)n 如果原先是45有名字的data.frame形式,則用n x=read.table(f:bookww.txt,header=T)n x GM VW HUNDA1993 1 2 31994 6 7 8讀取固定寬度格式數(shù)據(jù)read.fwf()n read.fwf(file, widths, sep=t, as.is = FALSE, skip = 0, row.names,col.names, n = -1, .)n # Read a table of fixed width formatted data into a data.frame. n
10、除了除了widths用來說明讀取字段的寬度外，用來說明讀取字段的寬度外， 選項(xiàng)與選項(xiàng)與read.table()基本相同。基本相同。n 現(xiàn)有一個名為現(xiàn)有一個名為data.txt的文件，內(nèi)容為的文件，內(nèi)容為n A 1.50 1.2n A 1.55 1.3n B 1.60 1.4n B 1.65 1.5n C 1.70 1.6n C 1.75 1.7n mydata mydata V1 V2 V31 A 1.50 1.22 A 1.55 1.33 B 1.60 1.44 B 1.65 1.55 C 1.70 1.66 C 1.75 1.7n ff cat(file=ff, 123456, 98765
11、4, sep=n) n read.fwf(ff, widths=c(1,2,3) n # 1 23 456 9 87 654 n read.fwf(ff, widths=c(1,-2,3) n # 1 456 9 654 n unlink(ff) n cat(file=ff, 123, 987654, sep=n) n read.fwf(ff, widths=c(1,0, 2,3) n # 1 NA 23 NA 9 NA 87 654 n unlink(ff) n cat(file=ff, 123456, 987654, sep=n) n read.fwf(ff, widths=list(c(
12、1,0, 2,3), c(2,2,2) n # 1 NA 23 456 98 76 54 n unlink(ff) V1 V2 V31 1 23 4562 9 87 654123456987654123987654123456987654讀SPSS,SAS等數(shù)據(jù)文件n 基礎(chǔ)包中并不包含所需的這些函數(shù)。首先載入程序包基礎(chǔ)包中并不包含所需的這些函數(shù)。首先載入程序包n library(foreign)n 現(xiàn)有現(xiàn)有SPSS數(shù)據(jù)文件數(shù)據(jù)文件(“educ_scores.sav”)，SAS 數(shù)據(jù)文件數(shù)據(jù)文件(“educ_scores .xpt”)，S-Plus數(shù)據(jù)文件數(shù)據(jù)文件(“educ_scores”)及
13、及Stata數(shù)據(jù)文件數(shù)據(jù)文件(educ_scores .dta).n rs rs rx rs rd read.table(ex.data,fill=T)-p; pn TITLE extra linen2 3 5 7n11 13 17 NAn write.table(p,foo.txt)n row.names(p)n 1 2 11n col.names(p)n 錯誤錯誤: 沒有沒有col.names這個函數(shù)這個函數(shù)n write.table(p,foo.txt,col.names=T)n write.table(p,foo.txt,col.names=F)n write.table(p,foo.
14、txt,col.names=NA) TITLE extra line2” 3 5 711 13 17 NA2 3 5 711 13 17 NA“” TITLE extra line2” 3 5 711 13 17 NAn write.csv(p,fo.txt)n , TITLE, extra, linen 2, 3, 5, 7n 11, 13, 17, NAn #分隔符全是逗號分隔符全是逗號n write.csv(p,fo.csv)n write.table(p,fo1.csv) #只有一列只有一列n write.csv uses . for the decimal point and a c
15、omma for the separator. n write.csv2 uses a comma for the decimal point and a semicolon for the separator, the Excel convention for CSV files in some Western European locales. 附：常用的數(shù)學(xué)函數(shù)n 簡單數(shù)學(xué)運(yùn)算有簡單數(shù)學(xué)運(yùn)算有:+,-,*,/,%*%,%(mod),%/%(整除整除)等等等等n abs , sign , log , log2, log10 , logb, expm1, log1p(x), sqrt , e
16、xp , sin , cos , tan , acos , asin, atan , cosh , sinh(雙曲正弦1/2ex-e(-x) )， tanhn round, floor, ceilingn gamma , lgamma, digamma and trigamma (the first and second derivatives of the logarithm of the gamma function ).n sum, prod, cumsum, cumprod, factorialn max, min, cummax, cummin, pmax, pmin, rangen
17、 mean, length, var, duplicated, uniquen union, intersect, setdiff, setequaln , =, , =, &, |, !還有還有 n 輸入輸出數(shù)據(jù): scan, read.table, save, load, write, write.tablen letters, LETTERS, month.namen list, matrix, array, cbind, rbind, mergen sort, order, sort.list, stack, unstack , reshape 2.14 Creating new fu
18、nctions n 控制流控制流n R是一個表達(dá)式語言，其任何一個語句都可以看成是一個表達(dá)式。是一個表達(dá)式語言，其任何一個語句都可以看成是一個表達(dá)式。n 若干個表達(dá)式放在一起組成一個復(fù)合表達(dá)式，用若干個表達(dá)式放在一起組成一個復(fù)合表達(dá)式，用括起來括起來n 1 分支語句分支語句n 包括包括if/else, switchif/else語句n if(cond) statement_1n if(cond) statement_1 else statement_2if( any(x = 0) ) y - log(1+x) else y - log(x)y - if( any(x x = c(3:-4) s
19、qrt(ifelse(x = 0, x, NA)(否則是否則是NA) ifelse(x = 0, sqrt(x), NA) if (x = 0) sqrt(x) else NAx = c(3:-4)n sqrt(ifelse(x = 0, x, NA) #清洗數(shù)據(jù)（不想要的清清洗數(shù)據(jù)（不想要的清洗為洗為NA）n 1 1.732051 1.414214 1.000000 0.000000 NA NA NA NAn ifelse(x = 0, sqrt(x), NA)n 1 1.732051 1.414214 1.000000 0.000000 NA NA NA NAn 警告信息：警告信息：n I
20、n sqrt(x) : 產(chǎn)生了產(chǎn)生了NaNsn if (x = 0) sqrt(x) else NAn 11.732051 1.414214 1.000000 0.00000 NaN NaN NaN NaNn 警告信息：警告信息：n 1: In if (x = 0) sqrt(x) else NA : 條件的長度大于一，因此只能用其第一元條件的長度大于一，因此只能用其第一元素素n 2: In sqrt(x) : 產(chǎn)生了產(chǎn)生了NaNsswitch語句n switch (statement, list)n x switch(x, 2+2, mean(1:10), rnorm(4)n 1 0.892
21、7328 -0.7827752 1.0772888 1.0632371n switch(2, 2+2, mean(1:10), rnorm(4)n 1 5.5n switch(6, 2+2, mean(1:10), rnorm(4)n NULLn y switch(y,fruit=banana,vegetable=broccoli)n 1 banana2 中止語句break與空語句next 3 循環(huán)語句n for循環(huán)語句循環(huán)語句n for (name in expr_1) expr_2n #構(gòu)造構(gòu)造4階階Hilbert矩陣n n-4; x for (i in 1:n)n for (j in 1
22、:n)n xi,j xn ,1 ,2 ,3 ,4n 1, 1.0000000 0.5000000 0.3333333 0.2500000n 2, 0.5000000 0.3333333 0.2500000 0.2000000n 3, 0.3333333 0.2500000 0.2000000 0.1666667n 4, 0.2500000 0.2000000 0.1666667 0.1428571n f - function(x) n + for(i in 1:10) n + res x=5; f(x)n cat(x,3) #output to the screenn x 3n cat(x,3
23、)-an x 3n anNULL5 * 1 = 5 5 * 2 = 10 5 * 3 = 15 5 * 4 = 20 5 * 5 = 25 5 * 6 = 30 5 * 7 = 35 5 * 8 = 40 5 * 9 = 45 5 * 10 = 50 while循環(huán)語句n while (condition) exprn編寫1000以內(nèi)的Fibonacci數(shù)n f-1f1=1; f2-1; i while (fi+fi+11000) nfi+2-fi+fi+1ni fn1 1 1 2 3 5 8 13 21 34 55 89 144n13 233 377 610 987n 循環(huán)循環(huán)(loops)
24、nfor() #若知道終止條件若知道終止條件for (變量變量 in 向量向量) 表達(dá)式表達(dá)式nwhile() #若無法知道運(yùn)行次數(shù)若無法知道運(yùn)行次數(shù)while(條件條件) 表達(dá)式表達(dá)式n兩者通?？梢赞D(zhuǎn)換兩者通?？梢赞D(zhuǎn)換試比較兩種方法試比較兩種方法for (i in 1:5) print (1:i) i=1while(i repeat exprnrepeat依賴依賴break跳出循環(huán)跳出循環(huán)n編寫1000以內(nèi)的Fibonacci數(shù)n f-1; f2-1; i repeat nfi+2-fi+fi+1ni=1000) breakn# if (fi+fi+1 sen Error: Object s
25、e not found2sn編輯自己的函數(shù)n se - function(x) n v - var(x)n n - length(x)n return(sqrt(v/n)n n General formn function() n 調(diào)用調(diào)用n se(c(45,2,3,5,76,2,4)n 1 11.10310n 注注:如果函數(shù)已保存在文件如果函數(shù)已保存在文件code.R里，里，”運(yùn)行運(yùn)行R腳本文件腳本文件”,或者將該文或者將該文件打開在主窗口運(yùn)行一下件打開在主窗口運(yùn)行一下,或或source(code.R)2sn編輯自己的函數(shù)n n basic.stats(c(45,2,4,46,43,65,N
26、A,6,-213,-3,-45)n $nn 1 11n $nNAsn 1 1n $meann 1 -5n $stdn 1 79.87768n $medn 1 5 basic.stats(c(45,2,4,46,43,65,NA,6,-213,-3,-45),more=T)$n1 11$nNAs1 1$mean1 -5$std1 79.87768$med1 5$skew1 -1.638217$kurt1 1.708149二分法n 編寫一個用二分法求非線性方程根的函數(shù)編寫一個用二分法求非線性方程根的函數(shù),并求方程并求方程n x3-x-1=0n 在區(qū)間在區(qū)間1,2內(nèi)的根內(nèi)的根.n 解解:先編寫一個用
27、二分法求非線性方程根的函數(shù)先編寫一個用二分法求非線性方程根的函數(shù),函數(shù)名函數(shù)名:fzero,但文件名但文件名可以任意，如可以任意，如bisect.R函數(shù)fzeron fzero0)n list(fail=finding root failes!) n elsen repeat n if (abs(b-a)eps) break;n x-(a+b)/2;n if (f(a)*f(x)0) b-x else a-xn n list(root=(a+b)/2, fun=f(x)n n 編寫要求根的函數(shù),用fzero求根n 參照參照fzero f fzero(f, 1, 2, 1e-6)n $rootn
28、 1 1.324718n $funn 1 -1.857576e-06n 如果在如果在0,1內(nèi)求根內(nèi)求根,則則n fzero(f,0,1)n $failn 1 finding root fails!直接求根直接求根 uniroot(f, c(1,2)$root1 1.324718$f.root1 -5.634261e-07$iter1 7$estim.prec1 6.103516e-05兩總體的T檢驗(yàn)n 兩樣本的方差相同，且未知，兩樣本的方差相同，且未知，T統(tǒng)計量為統(tǒng)計量為n twosam - function(y1, y2) n n1 - length(y1); n2 - length(y2)
29、n yb1 - mean(y1); yb2 - mean(y2)n s1 - var(y1); s2 - var(y2)n s - (n1-1)*s1 + (n2-1)*s2)/(n1+n2-2)n tst - (yb1 - yb2)/sqrt(s*(1/n1 + 1/n2)n tstn 2221122121211211()()(),nSnSXYTSnnSnn 其其中中n 參照參照 twosam A B twosam(A,B) n 結(jié)果為結(jié)果為n 1 3.472245定義新的二元函數(shù) %anything%n 定義定義 = exp(|x y|2/2),用運(yùn)算符號用運(yùn)算符號%!% 表示表示n %!
30、% fun1 ans ans ans - fun1(data=d, limit=20, graph=TRUE, data.frame=df)n 三者等價三者等價Newton法求根n牛頓法求方程的根，是將方程牛頓法求方程的根，是將方程f(x)=0在在xk泰勒展泰勒展開開,去掉去掉2階及階及2階以上的項(xiàng)得階以上的項(xiàng)得f(x)=f(xk)+f (xk)(x- xk),從而迭代公式為從而迭代公式為xk+1= xk- f(xk)/ f (xk).方程組的情形為方程組的情形為xk+1= xk J(xk)1f(xk), k= 0, 1, ,Newtons法程序n 參照參照 xk+1= xk J(xk)1f(
31、xk)n Newtons-function (fun, x, ep=1e-5, it_max=100)n index-0; k-1n while (k=it_max)n x1 - x; obj - fun(x);n x - x - solve(obj$J, obj$f);n norm - sqrt(x-x1) %*% (x-x1)n if (normep)n index-1; breakn n k-k+1n n obj funs-function(x)n f-c(x12+x22-5, (x1+1)*x2-(3*x1+1) J Newtons(funs,c(0,1)n $rootn 1 1 2n
32、 $itn 1 6n $indexn 1 1n $FunValn 1 1.598721e-14 6.217249e-15 4 遞歸函數(shù)n 遞歸即在函數(shù)自身內(nèi)定義遞歸即在函數(shù)自身內(nèi)定義/調(diào)用函數(shù)自身調(diào)用函數(shù)自身.n 計算數(shù)值積分計算數(shù)值積分 精度要求為精度要求為e=10-5n 運(yùn)算過程是先將區(qū)間分成二等分運(yùn)算過程是先將區(qū)間分成二等分,對二部分分別以函對二部分分別以函數(shù)值為上下底的梯形求面積數(shù)值為上下底的梯形求面積,然后相加然后相加,直致前后二次直致前后二次的結(jié)果差小于給定誤差限或遞歸次數(shù)達(dá)到要求時停的結(jié)果差小于給定誤差限或遞歸次數(shù)達(dá)到要求時停止止.51dxx 函數(shù)arean area - fun
33、ction(f, a, b, eps = 1.0e-06, lim = 10) n fun1 - function(f, a, b, fa, fb, a0, eps, lim, fun) n d - (a + b)/2; h - (b - a)/4; fd - f(d)n a1 - h * (fa + fd); a2 - h * (fd + fb)n if(abs(a0 - a1 - a2) eps | lim = 0)n return(a1+a2)n else n return(fun(f, a, d, fa, fd, a1, eps, lim - 1, fun)n + fun(f, d,
34、b, fd, fb, a2, eps, lim - 1, fun)n n n fa - f(a); fb - f(b); a0 f quad-area(f,1,5); quadn 1 1.609452適用于種群動態(tài)變化的非線性模型n Ricker 模型模型n Nt+1 = Nt exp r(1- Nt /k )n 這個模型廣泛地用于種群動態(tài)變化（這個模型廣泛地用于種群動態(tài)變化（population dynamics）的研究）的研究,特別是魚類的種群變化特別是魚類的種群變化. n 我們想用一個函數(shù)去模擬這個模型關(guān)于增長率我們想用一個函數(shù)去模擬這個模型關(guān)于增長率r和和初始群體大小初始群體大小N0的
35、變化情況的變化情況(承載能力承載能力K 常常設(shè)定常常設(shè)定為為1且以這個值作為默認(rèn)值且以這個值作為默認(rèn)值)的影響的影響; 結(jié)果將以種群結(jié)果將以種群大小相對時間的圖表示大小相對時間的圖表示. 我們還將設(shè)定一個可選項(xiàng)我們還將設(shè)定一個可選項(xiàng)允許用戶只顯示最后若干步中種群大小允許用戶只顯示最后若干步中種群大小(默認(rèn)所有默認(rèn)所有結(jié)果都會被繪制出來結(jié)果都會被繪制出來). Ricker函數(shù)的數(shù)值模擬n # Nt+1 = Nt exp r(1- Nt /k)n ricker - function(nzero, r, K=1, time=100, from=0, to=time) n N - numeric(ti
36、me+1) #產(chǎn)生序列，不是產(chǎn)生序列，不是asn N1 - nzeron for (i in 1:time) Ni+1 - Ni*exp(r*(1 - Ni/K)n Time - from:ton Population layout(matrix(1:3, 3, 1)n ricker(0.1, 1)n title(r = 1)n ricker(0.1, 2)n title(r = 2)n ricker(0.1, 3)n title(r = 3)函數(shù)小結(jié)n函數(shù)是一系列語句的組合，在函數(shù)是一系列語句的組合，在R中可以寫出自己的函數(shù)中可以寫出自己的函數(shù)n形式形式: 變量名變量名 = function
37、( 變量列表變量列表 ) 函數(shù)體函數(shù)體n函數(shù)引用函數(shù)引用: 變量名變量名(變量的值變量的值)n函數(shù)可以遞歸引用，但不提倡！函數(shù)可以遞歸引用，但不提倡！n例子例子 使用使用gamma函數(shù)求函數(shù)求n!factorial = function(n) + if (n=0) gamma(n+1) + else print(“Please input a positive integer!) + factorial(6) factorial(-6)附n用于處理錯誤的函數(shù)用于處理錯誤的函數(shù) 用于處理用戶輸入不正確的類用于處理用戶輸入不正確的類型而可能出現(xiàn)的錯誤型而可能出現(xiàn)的錯誤nwarning() 若錯誤不
38、嚴(yán)重以至影響整個計算若錯誤不嚴(yán)重以至影響整個計算nstop() 若錯誤可能導(dǎo)致計算中止若錯誤可能導(dǎo)致計算中止nprint() 顯示必要的信息顯示必要的信息nformatC() 數(shù)值作為字符串輸出數(shù)值作為字符串輸出ncat() 字符串聯(lián)，可以插入字符串聯(lián)，可以插入n(換行換行)及及t(tab鍵鍵)npaste() 字符粘貼字符粘貼(非字符型自動轉(zhuǎn)換非字符型自動轉(zhuǎn)換)n例子：例子：cat(R, is,a good, software.n) formatC(1/3, format = f, digits = 4) formatC(1/3, format = e, digits = 4)paste(1
39、:12) # 與與as.character(1:12)等價等價paste(A, 1:6, sep = )paste(today is, date()R中的編程n編程的重要性編程的重要性n一個統(tǒng)計軟件一個統(tǒng)計軟件(包包)應(yīng)該包羅萬象應(yīng)該包羅萬象所有統(tǒng)計方所有統(tǒng)計方法或技術(shù)，可以完成所有的任務(wù)？法或技術(shù)，可以完成所有的任務(wù)？n這是一種錯誤的觀點(diǎn)！這是一種錯誤的觀點(diǎn)！n太貴！太貴！n費(fèi)時！費(fèi)時！nSAS是一個世界上最為優(yōu)秀的統(tǒng)計分析軟件之一，是一個世界上最為優(yōu)秀的統(tǒng)計分析軟件之一，但但！當(dāng)然！當(dāng)然SAS具體編程能力，而具體編程能力，而SPSS更糟！更糟！n通過編程可以實(shí)現(xiàn)通過編程可以實(shí)現(xiàn)(開發(fā)開發(fā))
40、滿足自己需要的函數(shù)或宏滿足自己需要的函數(shù)或宏包包 n好的編程習(xí)慣好的編程習(xí)慣n 為了他人，更為你本人為了他人，更為你本人! 你的程序應(yīng)該具有你的程序應(yīng)該具有n可讀性可讀性(readability)n可理解性可理解性(understandability)n 習(xí)慣之一：行前縮進(jìn)習(xí)慣之一：行前縮進(jìn)(Indentation),在此推薦使用軟件在此推薦使用軟件nWinEdt = 已開發(fā)了基于已開發(fā)了基于WinEdt的的R使用平臺使用平臺n(X)Emacs =已開發(fā)了基于已開發(fā)了基于(X)Emacs的的R使用平臺使用平臺nUltraEditn 習(xí)慣之二：增加注釋習(xí)慣之二：增加注釋(Commenting),
41、它是你的幫助它是你的幫助 R中使用中使用#作為注釋語句的開始作為注釋語句的開始.n 習(xí)慣之三：變量的命名習(xí)慣之三：變量的命名,使用意義明確的名字，切忌使用人或?qū)櫸锏拿质褂靡饬x明確的名字，切忌使用人或?qū)櫸锏拿掷永?縮進(jìn)縮進(jìn) for(i in 1:1522) + gene.aov - aov(Intensity A + T + A * T,+ sub = G = i)+ res.mati, - residuals(gene.aov)+ coef.mati, for(i in 1:1522) + gene.aov - aov(Intensity A + T + A * T,+ sub = G
42、 = i)+ res.mati, - residuals(gene.aov)+ coef.mati, - coef(gene.aov)+ n實(shí)例演示實(shí)例演示n1Julian Faraway提供的一個簡短的入門提供的一個簡短的入門Introduction to R n使用使用 R GUIn使用使用 R Commandern使用使用 R_XEmacs Further reference: Moores The Basic Practice of Statistics, Second Edition (Freeman, 2000), n2. John Verzani 提供的統(tǒng)計學(xué)入門提供的統(tǒng)計學(xué)入門SimpleR - Using R for Introductory Statistics n使用使用R_WinEdt作業(yè)nP105 2.1, 2.2, P106 2.6, 2.7n試編出矩陣的試編出矩陣的Kronecker乘積的函數(shù)乘積的函數(shù).n 提示：矩陣提示：矩陣ApqBrs,定義定義A與與B的的Kronecker product為為Cprqs，其中，其中p,q,r,s均為正整數(shù)均為正整數(shù)1,11,21,2,12,22,1,2,.ssrrr sAbAbAbAbAbAbCAbAbAb
An introduction to R (III)計算機(jī)教學(xué)課件PPT
