R語言課件(簡略版)

1、2021/3/271 R語言語言 2021/3/27計(jì)算機(jī)與R語言2 目錄目錄 2021/3/273 第二章 R語言基礎(chǔ) 一、一個(gè)簡短的一、一個(gè)簡短的R會(huì)話會(huì)話 二、二、R的基本語法的基本語法 三、三、R的數(shù)據(jù)結(jié)構(gòu)的數(shù)據(jù)結(jié)構(gòu) 2021/3/274 一、一個(gè)簡短的一、一個(gè)簡短的R會(huì)話會(huì)話 數(shù)據(jù)的描述 ?mtcars mtcar:美國Motor Trend雜志收集的32輛汽車 （1973-1974車型）的11項(xiàng)指標(biāo)。 2021/3/275 一、一個(gè)簡短的一、一個(gè)簡短的R會(huì)話會(huì)話 數(shù)據(jù)的瀏覽 所有數(shù)據(jù) mtcars 前6個(gè)觀測(cè)值 head(mtcars) 后6個(gè)觀測(cè)值 tail(mtcars) 20

2、21/3/276 一、一個(gè)簡短的一、一個(gè)簡短的R會(huì)話會(huì)話 數(shù)據(jù)的編輯 方法一: data.entry(mtcars) 方法二: MTcarsfix(mtcars) 區(qū)別:edit()修改結(jié)果不存入mtcars中;fix()結(jié) 果保存至mtcars中。 2021/3/277 一、一個(gè)簡短的一、一個(gè)簡短的R會(huì)話會(huì)話 瀏覽變量信息,例如mpg先激活數(shù)據(jù)集mtcars attach(mtcars) mpg 1 21.0 21.0 22.8 21.4 18.7 18.1 14.3 24.4 22.8 19.2 17.8 16.4 17.3 15.2 15 10.4 10.4 14.7 32.4 30.4

3、 33.9 21.5 15.5 15.2 13.3 19.2 27.3 26.0 30.4 29 15.8 19.7 15.0 21.4 2021/3/278 二、R語言的基本語法 1、符號(hào) 1） 命令或運(yùn)算提示符 2） + 續(xù)行符 3） # 說明語句字符 2、賦值符號(hào) 1） 大于 2） =大于等于 4）= 小于等于 5）= 等于 6） != 不等于 比較運(yùn)算得到的結(jié)果是邏輯常量TRUE（真）和FALSE （假）。 2021/3/2711 5、邏輯運(yùn)算符 1） ! x邏輯非 當(dāng)前邏輯值取反運(yùn)算。 2）x mode(compare)-TRUE; mode(z)-1i 1 character 2

4、logical 3 complex 2021/3/2721 三、R語言的數(shù)據(jù)結(jié)構(gòu) 思考題思考題: : 下列不是數(shù)值型數(shù)據(jù)的是（）。 A123B.12.3C. 123 D.12e2.3 下列是字符型數(shù)據(jù)的是（）。 A. TRUEB.1e2C. 123 D.5+1i 下列字符型定界符用法錯(cuò)誤的是（）。 A ab c B. ab c C. abc D. ab c 下列是邏輯型數(shù)據(jù)的是（）。 A TrueB.FALSEC.InfD.NaN Inf是（）型的數(shù)據(jù)。 A.數(shù)值型B.NA型C.字符型D.邏輯型 2021/3/2722 三、R語言的數(shù)據(jù)結(jié)構(gòu) 統(tǒng)計(jì)實(shí)例 例:一家保險(xiǎn)公司收集到 由36個(gè)投保個(gè)人組

5、成的 隨機(jī)樣本,得到每個(gè)投保 人的年齡（周歲）數(shù)據(jù) 如表所示。試確定投保 人平均年齡90%的置信區(qū) 間。 233539273644 364246433133 425345544724 342839364440 394938344850 343945484532 2021/3/2723 三、R語言的數(shù)據(jù)結(jié)構(gòu) 注:置信區(qū)間是指數(shù)據(jù)可 靠程度的范圍。90%指 置信水平。 求置信水平下的置信區(qū) 間公式為:即x的均值加 減常用置信水平的正態(tài) 分布的臨界值乘以樣本 標(biāo)準(zhǔn)差與樣本個(gè)數(shù)開方 的差。 90%- =1.645 95%- =1.96 99%- =2.58 n s x 2021/3/2724 三、R語

6、言的數(shù)據(jù)結(jié)構(gòu) x-c(23,35,39,27,36,44, 36,42,46,43,31,33, 42,53,45,54,47,24,34,28,39,36,44,40, 39,49,38,34,48,50,34,39,45,48,45,32) m-mean(x) s-sd(x) z1-m+1.645*(s/sqrt(36) z2-m-1.645*(s/sqrt(36) z 1:10 1 1 2 3 4 5 6 7 8 9 10 a-1;b 1:(10-1) 1 1 2 3 4 5 6 7 8 9 1:10-1 1 0 1 2 3 4 5 6 7 8 9 2021/3/2729 三、R語言的數(shù)

7、據(jù)結(jié)構(gòu) 2、seq(n1,n2,by=n3,length=n4) #生成從n1到n2的向量, n3為步長,n4為生成元素的數(shù) 量 seq(1,10) 1 1 2 3 4 5 6 7 8 9 10 seq(1,5,by=0.5) 1 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 seq(1,10,length=11) 1 1.0 1.9 2.8 3.7 4.6 5.5 6.4 7.3 8.2 9.1 10.0 2021/3/2730 三、R語言的數(shù)據(jù)結(jié)構(gòu) 3、rep(n1,n2)#生成n1重復(fù)n2次的向量 rep(2,3) 1 2 2 2 rep(1:5,2) 1 1

8、2 3 4 5 1 2 3 4 5 rep(1:5,1:5) 1 1 2 2 3 3 3 4 4 4 4 5 5 5 5 5 rep(1:5,rep(2,5) 1 1 1 2 2 3 3 4 4 5 5 2021/3/2731 三、R語言的數(shù)據(jù)結(jié)構(gòu) 4、c()#直接輸出向量 x=c(42,7,64,9) length(x) 5、scan()#通過鍵盤建立向量 scan() 1: 1 3 6 8 #此行末打一空格后回車 5: #冒號(hào)后直接打回車 Read 4 items 1 1 3 6 8 2021/3/2732 三、R語言的數(shù)據(jù)結(jié)構(gòu) 6、sequence(n1)#生成從1開始到n1結(jié)束步長 為

9、1的向量 sequence(5) 1 1 2 3 4 5 sequence(c(5,5) 1 1 2 3 4 5 1 2 3 4 5 sequence(c(5,4) 1 1 2 3 4 5 1 2 3 4 2021/3/2733 思考題思考題: : 下列（）不能生成向量（1 5 4 3 2 1）。 A. 5:1 B. 6-1:5 C.seq(1,5,by=-1) D.6-sequence(5) 下列（）不能生成向量（1 1 2 3 4 5 1 2 3 4 5）。 A. rep(1:5,2) B. sequence(rep(5,2) C. 1:10%6 D. (0:9%5)+1 2021/3/2

10、734 三、R語言的數(shù)據(jù)結(jié)構(gòu) 2) 數(shù)值型向量的運(yùn)算 向量的運(yùn)算方法如下: （1）向量與一個(gè)常數(shù)的加、減、乘、除為向量的每一個(gè)元素與 此常數(shù)進(jìn)行加、減、乘、除; （2）向量的乘方（）與開方(sqrt)為每一個(gè)元素的乘方與開方, 這對(duì)log,exp,sin,cos,tan 等普通的運(yùn)算函數(shù)同樣適用; （3）同樣長度向量的加、減、乘、除等運(yùn)算為對(duì)應(yīng)元素進(jìn)行加、 減、乘、除等; （4）不同長度向量的加、減、乘、除遵從循環(huán)法則（recycling rule）,但要注意這種場(chǎng)合通常要求向量的長度為倍數(shù)關(guān)系,否 則會(huì)出現(xiàn)警告:“長向量并非是短向量的整數(shù)倍”。 2021/3/2735 三、R語言的數(shù)據(jù)結(jié)構(gòu)

11、循環(huán)法則:同一個(gè)表達(dá)式中的向量并不需要具有相同 的長度,如果它們的長度不同,表達(dá)式的結(jié)果是一個(gè) 與表達(dá)式中最長向量有相同長度的向量,表達(dá)式中較 短的向量會(huì)根據(jù)它的長度被重復(fù)使用若干次（不一 定是整數(shù)次）,直到與長度最長的向量相匹配,而常 數(shù)將被不斷重復(fù)。 2021/3/2736 三、R語言的數(shù)據(jù)結(jié)構(gòu) sqrt(c(2,4,25) 1 1.414214 2.000000 5.000000 1:2+1:4 1 2 4 4 6 1:4+1:7 1 2 4 6 8 6 8 10 警告信息: In 1:4 + 1:7 : 長的對(duì)象長 度不是短的對(duì)象長度 的整倍數(shù) 例子: 5+c(4,7,17) 1 9

12、12 22 5*c(4,7,17) 1 20 35 85 c(-1,3,-17)+c(4,7,17) 1 3 10 0 c(2,4,5)2 1 4 16 25 2021/3/2737 統(tǒng)計(jì)函數(shù)作用 Max(x)返回向量x中最大的元素 Min(x)返回向量x中最小的元素 Mean(x)計(jì)算樣本（向量）x的均值 Median(x)計(jì)算樣本（向量）x的中位數(shù) Var(x)計(jì)算樣本（向量）x的方差 Sd(x)計(jì)算向量x的標(biāo)準(zhǔn)差 Length(x)返回向量x的長度 Sum(x)給出向量x的總和 Cumsum(x)返回向量x和累積和（其第i個(gè)元素是從x1到xi的和） Cumprod(x)返回向量x和累積積

13、（其第i個(gè)元素是從x1到xi的積） Cummin(x)返回向量x和累積最小值（其第i個(gè)元素是從x1到xi的最小 值） Cummax(x)返回向量x和累積最大值（其第i個(gè)元素是從x1到xi的最大 值） Cov(x,y)計(jì)算樣本（向量）x與y的協(xié)方差 Cor(x,y)計(jì)算樣本（向量）x與y的相關(guān)系數(shù) 2021/3/2738 統(tǒng)計(jì)實(shí)例 兩變量線形相關(guān)系數(shù)的計(jì)算。 在-1和+1之間波動(dòng),其絕對(duì)值越接近1,兩個(gè)變量間的直線相關(guān) 越密切。 例題:trees數(shù)據(jù)集中樹木高度(Height)和木材體積(Volume)的相 關(guān)關(guān)系分析。 （1）繪制散點(diǎn)圖 （2）計(jì)算相關(guān)系數(shù)值。 的方差的方差 的協(xié)方差和 yx

14、y x 2021/3/2739 三、R語言的數(shù)據(jù)結(jié)構(gòu) 3）建立字符型向量 1. c() 直接輸入向量 例如例如: : c(a,b,1,1e2,TRUE)#數(shù)字、邏輯型數(shù)據(jù)被強(qiáng) 制轉(zhuǎn)換為字符串 1 a b 1 100 TRUE 2021/3/2740 2. paste() 接受任意字符,并順次組合連成字符串,在 默認(rèn)情況下,各字符串由一個(gè)空格分隔,可以通過 sep=string把它更改為其他字符,包括空字符串。 paste(a,1) 1 a 1 paste(c(a,b),1) 1 a 1 b 1 paste(c(a,b),c(1,2) 1 a 1 b 2 paste(app, le,sep=)

15、1 apple 2021/3/2741 三、R語言的數(shù)據(jù)結(jié)構(gòu) 例子 z-c(green,blue sky,-99) paste(1:12) paste(A, 1:6, sep = ) labsrep(a,3) 1 a a a rep(c(a,b),3) 1 a b a b a b rep(c(a,b,c),rep(2,3) 1 a a b b c c 2021/3/2743 4）建立邏輯型向量 1. c() 直接輸入向量 例如例如: : c(TRUE,FALSE,TRUE) 1 TRUE FALSE TRUE 2. rep() #生成重復(fù)的向量 例如例如: : rep(c(TRUE,FALSE

16、),2) 1 TRUE FALSE TRUE FALSE 3. 由運(yùn)算生成 例如例如: : 1:53 1 FALSE FALSE FALSE TRUE TRUE 2021/3/2744 三、R語言的數(shù)據(jù)結(jié)構(gòu) 5）建立因子型向量 一個(gè)因子(factor)或因子向量不僅包括分類變量本身,還包括變 量不同的可能水平。 factor(x,levels = sort(unique(x), na.last = TRUE),labels = levels, exclude =NA, ordered = is.ordered(x) 說明:這里x可以是數(shù)值型或字符型向量。 levels用來指定因子的水平（缺省值

17、是向量x中不同的值）; labels用來指定水平的名字; exclude表示從向量x中剔除的水平值; ordered是一個(gè)邏輯型選項(xiàng),用來指定因子的水平是否有次序。 2021/3/2745 三、R語言的數(shù)據(jù)結(jié)構(gòu) （1）將字符型向量轉(zhuǎn)換為因子 a a a 1 green blue green yellow Levels: blue green yellow 2021/3/2746 三、R語言的數(shù)據(jù)結(jié)構(gòu) （2）將數(shù)值型向量轉(zhuǎn)換成因子 b b b 1 1,2,3,1 Levels: 1 2 3 factor(rep(1:5,2) 1 1 2 3 4 5 1 2 3 4 5 Levels: 1 2 3

18、 4 5 2021/3/2747 三、R語言的數(shù)據(jù)結(jié)構(gòu) （3）將字符型因子轉(zhuǎn)換為數(shù)值型因子 a a levels(a) a 1 2 1 2 3 Levels: 1 2 3 4 ff ff 1 1 2 3 Levels: 1 2 3 2021/3/2748 三、R語言的數(shù)據(jù)結(jié)構(gòu) （4）將數(shù)值型因子轉(zhuǎn)換為字符型因子 b b levels(b) b 1 low middle high low Levels: low middle high ff ff 1 A B C Levels: A B C 注:函數(shù)levels()用來提取一個(gè)因子中可能的水平值, 2021/3/2749 三、R語言的數(shù)據(jù)結(jié)構(gòu) 例

19、如 ff ff 1 2 4 Levels: 2 3 4 5 levels(ff) 1 2 3 4 5 2021/3/2750 三、R語言的數(shù)據(jù)結(jié)構(gòu) （5）函數(shù)gl( )能產(chǎn)生規(guī)則的因子序列。 用法是 gl(k, n, length=產(chǎn)生數(shù)據(jù)的個(gè)數(shù), label=每個(gè)水平因子的名字) 其中k是水平數(shù),n是每個(gè)水平重復(fù)的次數(shù)。 例子 gl(3,5) 1 1 1 1 1 1 2 2 2 2 2 3 3 3 3 3 Levels: 1 2 3 2021/3/2751 三、R語言的數(shù)據(jù)結(jié)構(gòu) gl(3,5,length=30) gl(2,6,label=c(Male,Female) gl(2,10) gl

20、(2,1,length=20) gl(2,2,length=20) 2021/3/2752 統(tǒng)計(jì)函數(shù)統(tǒng)計(jì)函數(shù)作作 用用說明說明 1 1max(x)max(x)返回向量返回向量x x 中最大的元素中最大的元素 2 2min(x)min(x)返回向量返回向量x x 中最小的元素中最小的元素 3 3which.max(x)which.max(x)返回向量返回向量x x中最大元素的下標(biāo)中最大元素的下標(biāo) 4 4which.min(x)which.min(x)返回向量返回向量x x中最小元素的下標(biāo)中最小元素的下標(biāo) 5 5mean(x)mean(x)計(jì)算樣本（向量）計(jì)算樣本（向量）x x的均值的均值 6 6

21、median(x)median(x)計(jì)算樣本（向量）計(jì)算樣本（向量）x x的中位數(shù)的中位數(shù) 對(duì)長度為對(duì)長度為n的系列數(shù)據(jù)的系列數(shù)據(jù),根據(jù)數(shù)據(jù)大根據(jù)數(shù)據(jù)大 小排列得到的位于小排列得到的位于(n+1)/2位置上位置上 的數(shù)據(jù)。的數(shù)據(jù)。 7 7mad(x)mad(x)計(jì)算中位絕對(duì)離差計(jì)算中位絕對(duì)離差 8 8var(x)var(x)計(jì)算樣本（向量）計(jì)算樣本（向量）x x的方差的方差 各個(gè)數(shù)據(jù)與平均數(shù)之差的平方的平各個(gè)數(shù)據(jù)與平均數(shù)之差的平方的平 均數(shù)。均數(shù)。 9 9sd(x)sd(x)計(jì)算向量計(jì)算向量x x的標(biāo)準(zhǔn)差的標(biāo)準(zhǔn)差 各數(shù)據(jù)偏離平均數(shù)的距離的平均數(shù)各數(shù)據(jù)偏離平均數(shù)的距離的平均數(shù), 1010rang

22、e(x)range(x)返回長度為返回長度為2 2的向量的向量: : c(min(x),max(x)c(min(x),max(x) 1111IQR(x)IQR(x)計(jì)算樣本的四分位數(shù)極差計(jì)算樣本的四分位數(shù)極差 1212quantile(x)quantile(x)計(jì)算樣本常用的分位數(shù)計(jì)算樣本常用的分位數(shù) 1313summary(x)summary(x)計(jì)算常用的描述性統(tǒng)計(jì)量（最小、最大、計(jì)算常用的描述性統(tǒng)計(jì)量（最小、最大、 平均值、中位數(shù)和四分位數(shù)）平均值、中位數(shù)和四分位數(shù)） 2021/3/2753 1414length(x)length(x)返回向量返回向量x x的長度的長度 1515sum(

23、x)sum(x)給出向量給出向量x x的總和的總和 1616prod(x)prod(x)給出向量給出向量x x的乘積的乘積 1717rev(x)rev(x)取向量取向量x x的逆序的逆序 1818sort(x)sort(x)將向量將向量x x按升序排序按升序排序, ,選項(xiàng)選項(xiàng) decreasing=TRUEdecreasing=TRUE表示降序表示降序 1919order(x)order(x)返回返回x x的秩（升序）的秩（升序）, ,選項(xiàng)選項(xiàng) decreasing=TRUEdecreasing=TRUE, ,得到降序的秩得到降序的秩 2020rank(x)rank(x)返回返回x x的秩的秩

24、 2121cumsum(x)cumsum(x)返回向量返回向量x x 和累積和和累積和 2222cumprob(x)cumprob(x)返回向量返回向量x x和累積積和累積積 2323cummin(x)cummin(x)返回向量返回向量x x和累積最小值和累積最小值 2424cummax(x)cummax(x)返回向量返回向量x x和累積最大值和累積最大值 2525var(x,y)var(x,y)計(jì)算樣本（向量）計(jì)算樣本（向量）x x與與y y的方差的方差 2626cov(x,y)cov(x,y)計(jì)算樣本（向量）計(jì)算樣本（向量）x x與與y y的協(xié)方差的協(xié)方差 2727cor(x,y)cor(

25、x,y)計(jì)算樣本（向量）計(jì)算樣本（向量）x x與與y y的相關(guān)系數(shù)的相關(guān)系數(shù) 2828outer(x,y)outer(x,y)計(jì)算樣本（向量）計(jì)算樣本（向量）x x 與與y y的外積的外積 29runif(n,min=n1,max=n2) 產(chǎn)生產(chǎn)生n個(gè)個(gè)n1和和n2之間的隨機(jī)數(shù)之間的隨機(jī)數(shù) 2021/3/2754 三、R語言的數(shù)據(jù)結(jié)構(gòu) 6）向量的下標(biāo)與子集的提取 選擇一個(gè)向量的子集可以通過在其名稱后面追加一 個(gè)方括號(hào)中的索引向量來完成。 （1）正整數(shù)向量 x1:10 選取了x的前10個(gè)元素（假設(shè)x的長度不小于10）。 xc(1,4) 取出了向量x的第1和第4個(gè)元素。 xx2 y fruit n

26、ames(fruit) fruit lunch lunch 2021/3/2756 三、R語言的數(shù)據(jù)結(jié)構(gòu) （4）邏輯向量 x x10 xx10 xx10 numeric(0) xx10 x 2021/3/2757 思考題思考題: : 1.已知a-1:100,若要取出向量a中的偶數(shù)位置上的元素,表達(dá) 式應(yīng)為（）。 A. a2:100B. a/50 C.aseq(2,100,by=2) D.a/2 2. 已知a5D.a5 3. 表達(dá)式max(c(1:5,7:10)6:8)的結(jié)果是（）。 A. 1B.6C.7D.9 4. 已知a-c(1:5,7:10),則max(a6:8-3:1)的結(jié)果是（）。 A

27、. 1B.3C.9D.10 5. 已知a array(data,dim,dimnames) 其中data為一向量,其元素用于構(gòu)建數(shù)組;dim為數(shù)組的維數(shù)向 量;dimnames為由各維的名稱構(gòu)成的向量（字符型）,缺省為 空。 2021/3/2759 三、R語言的數(shù)據(jù)結(jié)構(gòu) 例子 A A , , 1 ,1 ,2 ,3 ,4 1, 1 4 7 10 2, 2 5 8 11 3, 3 6 9 12 , , 2 ,1 ,2 ,3 ,4 1, 1 4 7 10 2, 2 5 8 11 3, 3 6 9 12 其中,dim=c(3,4,2)的 第一個(gè)3代表有3行, 第二個(gè)4代表有4列, 第三個(gè)2代表有2組。

28、 A B B dim(A) dimnames(A) A colnames(A) rownames(A) dimnames(A) 2021/3/2761 三、R語言的數(shù)據(jù)結(jié)構(gòu) 例子 A A A A A A 2021/3/2762 三、R語言的數(shù)據(jù)結(jié)構(gòu) 2）矩陣的建立 由于矩陣在數(shù)學(xué)及統(tǒng)計(jì)中的特殊性,在R中最為常用 的是使用命令matrix( )建立矩陣。格 式:matrix(data=NA,nrow=1,ncol=1,byrow=FALSE,dimna mes=NULL) 其中data是一個(gè)向量數(shù)據(jù),nrow是矩陣的行數(shù),ncol是 矩陣的列數(shù),當(dāng)byrow=TRUE時(shí),生成矩陣的數(shù)據(jù)按行 放置

29、,缺省時(shí)相當(dāng)于byrow=FALSE,數(shù)據(jù)按列放置。 dimnames是數(shù)組維的名字,缺省時(shí)為空。 2021/3/2763 三、R語言的數(shù)據(jù)結(jié)構(gòu) 例子 X X X X A A 注意:循環(huán)準(zhǔn)則仍然適用于matrix(),但是要求數(shù)據(jù)項(xiàng) 的個(gè)數(shù)等于矩陣的列數(shù)的倍數(shù),否則會(huì)出現(xiàn)警告。 2021/3/2764 三、R語言的數(shù)據(jù)結(jié)構(gòu) 3）數(shù)組與矩陣的下標(biāo)（index）與子集（元素）的提取 同向量的下標(biāo)一樣,矩陣與數(shù)組的下標(biāo)可以使用正整 數(shù)、負(fù)整數(shù)和邏輯表達(dá)式,從而實(shí)現(xiàn)子集的提取或修 改。 2021/3/2765 三、R語言的數(shù)據(jù)結(jié)構(gòu) X X X2,2 #取出一個(gè)元素 X2, #取出一行元素 X,2 #取

30、出一列元素 X-1, #去掉第1個(gè)行 X,-2 #去掉第2個(gè)列 X,3Xis.na(x) X 2021/3/2766 思考題思考題: : 已知矩陣a X X t(X) 提取對(duì)角元diag( ) X X diag(X) diag(3) #生成3行的單位矩陣 2021/3/2768 三、R語言的數(shù)據(jù)結(jié)構(gòu) 按行合并rbind( )與按列合 并cbind( ) m1 m2 rbind(m1,m2) cbind(m1,m2) 矩陣的逐元乘積“*” m2*m2 ,1 ,2 ,3 ,4 1, 1 1 2 2 2, 1 1 2 2 ,1 ,2 1, 1 1 2, 1 1 3, 2 2 4, 2 2 ,1 ,2

31、 1, 4 4 2, 4 4 ,1 ,2 1, 1 1 2, 1 1 ,1 ,2 1, 2 2 2, 2 2 m1m2 rbind(m1,m2) cbind(m1,m2) m2*m2 2021/3/2769 三、R語言的數(shù)據(jù)結(jié)構(gòu) 矩陣的代數(shù)乘積“%*%” 【注】矩陣的代數(shù)乘法 設(shè) 分別是m n, n p矩陣, 則矩陣A與B的乘積是一m p矩 陣: 2021/3/2770 例1: m3 m4 m3%*%m4 ,1 ,2 1, 11 19 2, 16 28 例2: rbind(m1,m2)%*% cbind(m1,m2) ,1 ,2,3 ,4 1, 2 2 4 4 2, 2 2 4 4 3, 4

32、4 8 8 4, 4 4 8 8 cbind(m1,m2) %*% rbind(m1,m2) ,1 ,2 1, 10 10 2, 10 10 2021/3/2771 三、R語言的數(shù)據(jù)結(jié)構(gòu) 方陣的行列式det( ) X X ,1 ,2 1, 1 3 2, 2 4 det(X) 1 -2 說明:行列式相當(dāng)于主對(duì)角線（左上至右下）元素的乘積 減去副對(duì)角線（右上至左下）元素的乘積。 2021/3/2772 三、R語言的數(shù)據(jù)結(jié)構(gòu) 其它函數(shù) crossprod( ):交叉乘積函數(shù); eigen( ):特征根與特征向量函數(shù); qr( ):QR分解函數(shù)。 solve ():逆矩陣函數(shù)。 例:x-matrix(

33、1:4,2,2) y-solve(x) y 2021/3/2773 思考題思考題: : 已知a-matrix(1:4,2,2);bf=matrix(c(11.3262,36.241,5.1928,3.550,406.82,11.4638, 40.382,5.3782,4.120,415.92,11.5993,45.178,5.7170,5.70 0,429.03,513.2178,2495.299,456.2197,1667.402,786.45, 613.3035,3006.700,542.1962,1778.898,790.48),5,5,byro w=T) f cor(f) #計(jì)算相關(guān)系

34、數(shù) pairs(f) #矩陣散點(diǎn)圖 2021/3/2776 三、R語言的數(shù)據(jù)結(jié)構(gòu) （2）矩陣的統(tǒng)計(jì)運(yùn)算 函數(shù)1:apply(X,MARGIN,FUN) 其中X為參與運(yùn)算的矩陣,FUN為上面的一個(gè)函 數(shù)或“+”、“-”、“*”、“/”（必須放在引號(hào) 中）,MARGIN=1表示按行計(jì)算,MARGIN=2表示按 列計(jì)算。 函數(shù)2:sweep(X,MARGIN,STATS,FUN) 表示從矩陣X中按MATGIN計(jì)算STATS,并從X中 除去。 FUN默認(rèn)是“-”。 2021/3/2777 三、R語言的數(shù)據(jù)結(jié)構(gòu) 例1:求均值,中位數(shù)等 m apply(m,MARGIN=1,FUN=mean) #求各行的

35、均值 apply(m,MARGIN=2,FUN=mean) #求各列的均值 2021/3/2778 三、R語言的數(shù)據(jù)結(jié)構(gòu) 減去中位數(shù) row.med-apply(m,MARGIN=1,FUN=median) sweep(m,MARGIN=1,STATS=row.med,FUN=-) 2021/3/2779 三、R語言的數(shù)據(jù)結(jié)構(gòu) 4、列表與數(shù)據(jù)框的建立 1）列表的建立 列表是一種特別的對(duì)象集合,它的元素也由序 號(hào)（下標(biāo)）區(qū)分,但是各元素的類型可以是任意對(duì) 象,不同元素不必是同一類型。元素本身允許時(shí)其 它復(fù)雜數(shù)據(jù)類型,比如,列表的一個(gè)元素也允許是 列表。 2021/3/2780 三、R語言的數(shù)據(jù)結(jié)

36、構(gòu) 格式: Lst Lst Lst $name 1 Fred $wife 1 Mary $no.children 1 3 $child.ages 1 4 7 9 2021/3/2782 列表元素總可以用“列表名下標(biāo)”的格式引用, Lst2 1 Mary Lst42 1 7 Lstname 1 Fred Lstchild.age NULL Lstchild.ages 1 4 7 9 Lst$name 1 Fred Lst$wife 1 Mary 2021/3/2783 三、R語言的數(shù)據(jù)結(jié)構(gòu) 列表的元素可以修改,只要把元素引用賦值即可。 Lst$name Lst$income list.ABC-c(list.A,list.B,list.C) 2021/3/2784 三、R語言的數(shù)據(jù)結(jié)構(gòu) 2）數(shù)據(jù)框的建立 統(tǒng)計(jì)分析中一個(gè)完整的數(shù)據(jù)集通常是由若干個(gè) 變量的若干個(gè)觀測(cè)值組成的,在R中稱為數(shù)據(jù)框。 維數(shù):二維。 2021/3/2785 三、R語言的數(shù)據(jù)結(jié)構(gòu) 矩陣: 數(shù)據(jù)框: 二維數(shù)據(jù)對(duì)象相同點(diǎn) 不同點(diǎn) 數(shù)據(jù)框 二維 行與列的意義不同:列表示變量;行 表示觀測(cè)值。 矩陣、二維數(shù)組行列意義相同 2021/3/2786 三、R語言的數(shù)據(jù)結(jié)構(gòu) （1）數(shù)據(jù)框的生