统计学R语言实验

1、文档编码 : CY4E4L10G6M5 HO10G1D5Q9P5 ZD1L3P5K5A3. -试验 5R 绘图二一、试验目的：1. 娴熟把握描述性统计分析中常用的统计量；2. 把握 R 语言的高水平作图命令；3. 把握 R 语言的低水平作图命令；4. 把握多元数据的三个数据特点：均值向量、协方差矩阵、相关系数矩阵；二、试验容：练习：要求：完成练习并粘贴运行截图到文档相应位置截图方法见下 ,并将全部自己输入文字的字体颜色设为 红色 包括后面的摸索及小结,答复摸索题,简要书写试验小结；修改本文档名为“ 本人完整学号 1,其中 1 表示第 1 次试验,以后更换为 2,3,.；如文件名为“130554

2、3109立 1,表示学号为 1305543109的立同学的第 1 次试验,留意文件名中没有空格及任何其它字符；最终连同数据文件、源程序文件等假如有的话,本次试验没有,一起压缩打包发给课代表,压缩包的文件名同上；截图方法：法 1：调整需要截图的窗口至适宜的大小,并使该窗口为当前激活窗口即该窗口在屏幕最前方,按住键盘 Alt 键空格键两侧各有一个不放,再按键盘右上角的截图键通常印有“ 印屏幕或“Pr Scrn等字符,即完成截图；再粘贴到 word 文档的相应位置即可；法 2：利用 QQ 输入法的截屏工具；点击 QQ 输入法工具条最右边的“ 扳手图标,选择其中的“ 截屏工具； 1. 自行完成教材 P

3、130 页开场的 3.3-3.4 节中的例题；2.习题 3.5小白鼠在接种了 3 种不同菌型的伤寒杆菌后的存活天数如下表所示,试绘出数据的箱线图接受两种方法,一种是 plot 语句,另一种是 boxplot 语句来判定小白鼠被注射3 种菌型后的平均存活天数有无显著性差异？白鼠试验数据菌型存活天数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 注：平均存活天数其实应当是均值,但当不知道总体的分布时,通常是比拟中位数；因此,这里考虑用箱线图中的中位数来进展比拟；解：源代码：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 运行截图：结论：具有相关关系的两个变量的散点图要么是从左下角到右上角即正相关关系,要么是从左上角到右下角即负相关关系 ；但从上图可知全部的图中都没有这样的趋势,故均不相关；4.习题 3.7某校测得 19 名同学的四项指

5、标,性别、年龄、 身高cm和体重 lb,具体数据如下表所示；. （1）试绘出体重对于身高的散点图；coplotab|c （2）绘出不同性别情形下,体重与身高的散点图；. word.zl-. -（3）绘出不同年龄段的体重与身高的散点图；coplotab|c+d （4）绘出不同性别和不同年龄段的体重与身高的散点图；学号性别年龄身高体重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.习题 3.8画出函数 z = x 4- 2x 2y + x 2- 2xy + 2y 2 + 9 x/2 - 4y + 4 在区域 - 2 x 3,- 1 y 7 上的三维网格曲面和二维等值线,其中 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.习题 3.9用 Pearson相关检验法检验习题3.7 中的身高与体重是否相关；解：源代码：cor.test身高 ,

9、体重 运行截图：结论：其 P 值为 7.887e-070.05,拒绝原假设,认为变量身高与体重相关；7.习题 3.10绘出例 3.17 中 48 号求职者数据的星图； 数据见文档最终（1）以 15 项自变量 FL,APP, ,SUIT 为星图的轴；（2）以 G1,G2, , G5 为星图的轴； 提示：先通过这些星图, 能否说明应选哪 在作图中可适当调整各种参数；6 名应聘者；为使星图能够充分反映应聘者的情形,提示：第 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 列数据赋值一个新的数据框,再按 X$G1X$G2 X$G3 X$G4 X$G5rtstarsrt 运行截图：. . word.zl-. -摸索：1. R 语言的高水平作图函数与低水平作图函数的主要区分是什么？所谓高水平函数是与低水平作图函数相对应的,即全部的绘图函数均可产生图形,可以有坐标轴,以及图和坐标轴的说明文字；2. plot函数除了能画出散点图外,仍能画什么图？能画出曲线图吗？曲线图,箱线图；能3. 多元数据主要有哪三个数据特点指标？均值向量、协方差矩阵、相

11、关系数矩阵4. 相关检验的原假设是什么？总体 X,Y 不相关5. 多元数据的图形表示方法主要有哪几种图？轮廓图,星图,调和曲线图三、试验小结必写,但字数不限本次试验, 主要是在星图方面有些阻碍,用不同的函数读取的结果会不一样,有时仍读不出来,最终仍是把数据存入 excel 中,才能成功；要学会描述性统计分析中常用的统计量,高水平作图函数与低水平图函数的区分；. . word.zl-. -附：48 名应聘者的得分情形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