統(tǒng)計學(xué)R語言實驗

1、文檔編碼 : CY4E4L10G6M5 HO10G1D5Q9P5 ZD1L3P5K5A3. -試驗 5R 繪圖二一、試驗?zāi)康模?. 嫻熟把握描述性統(tǒng)計分析中常用的統(tǒng)計量；2. 把握 R 語言的高水平作圖命令；3. 把握 R 語言的低水平作圖命令；4. 把握多元數(shù)據(jù)的三個數(shù)據(jù)特點：均值向量、協(xié)方差矩陣、相關(guān)系數(shù)矩陣；二、試驗容：練習(xí)：要求：完成練習(xí)并粘貼運行截圖到文檔相應(yīng)位置截圖方法見下 ,并將全部自己輸入文字的字體顏色設(shè)為 紅色 包括后面的摸索及小結(jié),答復(fù)摸索題,簡要書寫試驗小結(jié)；修改本文檔名為“ 本人完整學(xué)號 1,其中 1 表示第 1 次試驗,以后更換為 2,3,.；如文件名為“130554

2、3109立 1,表示學(xué)號為 1305543109的立同學(xué)的第 1 次試驗,留意文件名中沒有空格及任何其它字符；最終連同數(shù)據(jù)文件、源程序文件等假如有的話,本次試驗沒有,一起壓縮打包發(fā)給課代表,壓縮包的文件名同上；截圖方法：法 1：調(diào)整需要截圖的窗口至適宜的大小,并使該窗口為當(dāng)前激活窗口即該窗口在屏幕最前方,按住鍵盤 Alt 鍵空格鍵兩側(cè)各有一個不放,再按鍵盤右上角的截圖鍵通常印有“ 印屏幕或“Pr Scrn等字符,即完成截圖；再粘貼到 word 文檔的相應(yīng)位置即可；法 2：利用 QQ 輸入法的截屏工具；點擊 QQ 輸入法工具條最右邊的“ 扳手圖標(biāo),選擇其中的“ 截屏工具； 1. 自行完成教材 P

3、130 頁開場的 3.3-3.4 節(jié)中的例題；2.習(xí)題 3.5小白鼠在接種了 3 種不同菌型的傷寒桿菌后的存活天數(shù)如下表所示,試?yán)L出數(shù)據(jù)的箱線圖接受兩種方法,一種是 plot 語句,另一種是 boxplot 語句來判定小白鼠被注射3 種菌型后的平均存活天數(shù)有無顯著性差異？白鼠試驗數(shù)據(jù)菌型存活天數(shù)4 1 2 4 3 2 4 7 7 2 2 5 2 5 6 8 5 10 7 12 12 6 6 3 10 3 7 11 6 6 7 9 5 5 10 6 注：平均存活天數(shù)其實應(yīng)當(dāng)是均值,但當(dāng)不知道總體的分布時,通常是比擬中位數(shù)；因此,這里考慮用箱線圖中的中位數(shù)來進展比擬；解：源代碼：y-c2,4,3,

4、2,4,7,7,2,2,5,4, 5,6,8,5,10,7,12,12,6,6, 7,11,6,6,7,9,5,5,10,6,3,10 f-factorcrep1,11,rep2,10,rep3,12 plotf,y . . word.zl-. -A-c2,4,3,2,4,7,7,2,2,5,4 B-c5,6,8,5,10,7,12,12,6,6 C df df pairsdf 運行截圖：結(jié)論：具有相關(guān)關(guān)系的兩個變量的散點圖要么是從左下角到右上角即正相關(guān)關(guān)系,要么是從左上角到右下角即負(fù)相關(guān)關(guān)系 ；但從上圖可知全部的圖中都沒有這樣的趨勢,故均不相關(guān)；4.習(xí)題 3.7某校測得 19 名同學(xué)的四項指

5、標(biāo),性別、年齡、 身高cm和體重 lb,具體數(shù)據(jù)如下表所示；. （1）試?yán)L出體重對于身高的散點圖；coplotab|c （2）繪出不同性別情形下,體重與身高的散點圖；. word.zl-. -（3）繪出不同年齡段的體重與身高的散點圖；coplotab|c+d （4）繪出不同性別和不同年齡段的體重與身高的散點圖；學(xué)號性別年齡身高體重01 Alice F 13 56.5 84.0 02 Becka F 13 65.3 98.0 03 Gail F 14 64.3 90.0 04 Karen F 12 56.3 77.0 05 Kathy F 12 59.8 84.5 06 Mary F 15 66

6、.5 112.0 07 Sandy F 11 51.3 50.5 08 Sharon F 15 62.5 112.5 09 Tammy F 14 62.8 102.5 10 Alfred M 14 69.0 112.5 11 Duke M 14 63.5 102.5 12 Guido M 15 67.0 133.0 13 James M 12 57.3 83.0 14 Jeffrey M 13 62.5 84.0 15 John M 12 59.0 99.5 16 Philip M 16 72.0 150.0 17 Robert M 12 64.8 128.0 18 Thomas M 11 5

7、7.5 85.0 19 William M 15 66.5 112.0 解：源代碼： ef attachef plot 體重 身高 coplot 體重 身高 | 性別 coplot 體重 身高 | 年齡 coplot 體重 身高 | 年齡 + 性別 運行截圖：. . word.zl-. -. . word.zl-. -. . word.zl-. -5.習(xí)題 3.8畫出函數(shù) z = x 4- 2x 2y + x 2- 2xy + 2y 2 + 9 x/2 - 4y + 4 在區(qū)域 - 2 x 3,- 1 y 7 上的三維網(wǎng)格曲面和二維等值線,其中 x 與 y各點之間的間隔為 0.05,等值線的值

8、分別為 0, 1, 2, 3, 4, 5, 10, 15, 20, 30, 40, 50, 60, 80, 100,共 15 條；注：在三維圖形中選擇適宜的角度解：源代碼：x-seq-2,3,0.05 y-seq-1,7,0.05 f-functionx,y x4-2*x2*y+x2-2*x*y+2*y2+9*x/2-4*y+4 z perspx,y,z,theta=30,phi=30,expand=0.7,col=lightblue 運行截圖：. . word.zl-. -6.習(xí)題 3.9用 Pearson相關(guān)檢驗法檢驗習(xí)題3.7 中的身高與體重是否相關(guān)；解：源代碼：cor.test身高 ,

9、體重 運行截圖：結(jié)論：其 P 值為 7.887e-070.05,拒絕原假設(shè),認(rèn)為變量身高與體重相關(guān)；7.習(xí)題 3.10繪出例 3.17 中 48 號求職者數(shù)據(jù)的星圖； 數(shù)據(jù)見文檔最終（1）以 15 項自變量 FL,APP, ,SUIT 為星圖的軸；（2）以 G1,G2, , G5 為星圖的軸； 提示：先通過這些星圖, 能否說明應(yīng)選哪 在作圖中可適當(dāng)調(diào)整各種參數(shù)；6 名應(yīng)聘者；為使星圖能夠充分反映應(yīng)聘者的情形,提示：第 2小題,先利用教材P153 頁的公式,運算出48 個求職者這5 個變量的值,公式如下：G1-SC+LC+SMS+DRV+AMB+GSP+POT/7 . . word.zl-. -

10、G2-FL+EXP+SUIT/3 G3-LA+HON+KJ/3 G4-AA G5 X starsX attachX 48 行 5 列數(shù)據(jù)賦值一個新的數(shù)據(jù)框,再按 X$G1X$G2 X$G3 X$G4 X$G5rtstarsrt 運行截圖：. . word.zl-. -摸索：1. R 語言的高水平作圖函數(shù)與低水平作圖函數(shù)的主要區(qū)分是什么？所謂高水平函數(shù)是與低水平作圖函數(shù)相對應(yīng)的,即全部的繪圖函數(shù)均可產(chǎn)生圖形,可以有坐標(biāo)軸,以及圖和坐標(biāo)軸的說明文字；2. plot函數(shù)除了能畫出散點圖外,仍能畫什么圖？能畫出曲線圖嗎？曲線圖,箱線圖；能3. 多元數(shù)據(jù)主要有哪三個數(shù)據(jù)特點指標(biāo)？均值向量、協(xié)方差矩陣、相

11、關(guān)系數(shù)矩陣4. 相關(guān)檢驗的原假設(shè)是什么？總體 X,Y 不相關(guān)5. 多元數(shù)據(jù)的圖形表示方法主要有哪幾種圖？輪廓圖,星圖,調(diào)和曲線圖三、試驗小結(jié)必寫,但字?jǐn)?shù)不限本次試驗, 主要是在星圖方面有些阻礙,用不同的函數(shù)讀取的結(jié)果會不一樣,有時仍讀不出來,最終仍是把數(shù)據(jù)存入 excel 中,才能成功；要學(xué)會描述性統(tǒng)計分析中常用的統(tǒng)計量,高水平作圖函數(shù)與低水平圖函數(shù)的區(qū)分；. . word.zl-. -附：48 名應(yīng)聘者的得分情形ID FL APP AA LA SC LC HON SMS EXP DRV AMB GSP POT KJ SUIT 1 6 7 2 5 8 7 8 8 3 8 9 7 5 7 10

12、2 9 10 5 8 10 9 9 10 5 9 9 8 8 8 10 3 7 8 3 6 9 8 9 7 4 9 9 8 6 8 10 4 5 6 8 5 6 5 9 2 8 4 5 8 7 6 5 5 6 8 8 8 4 4 9 5 8 5 5 8 8 7 7 6 7 7 7 6 8 7 10 5 9 6 5 8 6 6 6 7 9 9 8 8 8 8 8 8 10 8 10 8 9 8 10 8 9 9 9 8 9 9 8 8 10 9 10 9 9 9 10 9 9 9 7 8 8 8 8 5 9 8 9 8 8 8 10 10 4 7 10 2 10 10 7 10 3 10 10

13、10 9 3 10 11 4 7 10 0 10 8 3 9 5 9 10 8 10 2 5 12 4 7 10 4 10 10 7 8 2 8 8 10 10 3 7 13 6 9 8 10 5 4 9 4 4 4 5 4 7 6 8 14 8 9 8 9 6 3 8 2 5 2 6 6 7 5 6 15 4 8 8 7 5 4 10 2 7 5 3 6 6 4 6 16 6 9 6 7 8 9 8 9 8 8 7 6 8 6 10 17 8 7 7 7 9 5 8 6 6 7 8 6 6 7 8 18 6 8 8 4 8 8 6 4 3 3 6 7 2 6 4 19 6 7 8 4 7 8

14、 5 4 4 2 6 8 3 5 4 20 4 8 7 8 8 9 10 5 2 6 7 9 8 8 9 21 3 8 6 8 8 8 10 5 3 6 7 8 8 5 8 22 9 8 7 8 9 10 10 10 3 10 8 10 8 10 8 23 7 10 7 9 9 9 10 10 3 9 9 10 9 10 8 24 9 8 7 10 8 10 10 10 2 9 7 9 9 10 8 25 6 9 7 7 4 5 9 3 2 4 4 4 4 5 4 26 7 8 7 8 5 4 8 2 3 4 5 6 5 5 6 27 2 10 7 9 8 9 10 5 3 5 6 7 6 4

15、 5 28 6 3 5 3 5 3 5 0 0 3 3 0 0 5 0 29 4 3 4 3 3 0 0 0 0 4 4 0 0 5 0 30 4 6 5 6 9 4 10 3 1 3 3 2 2 7 3 31 5 5 4 7 8 4 10 3 2 5 5 3 4 8 3 32 3 3 5 7 7 9 10 3 2 5 3 7 5 5 2 33 2 3 5 7 7 9 10 3 2 2 3 6 4 5 2 34 3 4 6 4 3 3 8 1 1 3 3 3 2 5 2 35 6 7 4 3 3 0 9 0 1 0 2 3 1 5 3 36 9 8 5 5 6 6 8 2 2 2 4 5 6 6 3 37 4 9 6 4 10 8 8 9 1 3 9 7 5 3 2 38 4 9 6 6 9 9 7 9 1 2 10 8 5 5 2 . . word.zl-. -39 10 6 9 10 9 10 10 10 10 10 8 10 10 10 10 40 10 6 9 10 9 10 10 10 10 10 10 10 10 10 10 41 10 7 8 0 2 1 2 0 10 2 0 3 0 0 10 42 10 3