应用统计分析实验R软件.ppt

《应用统计分析实验R软件.ppt》由会员分享，可在线阅读，更多相关《应用统计分析实验R软件.ppt（52页珍藏版）》请在三一文库上搜索。

1、应用统计分析实验 R软件 规 太 艰 食 赚 诺 介 泞 哭 枕 脾 萝 孵 换 停 叮 暑 挑 愿 宝 版 奥 痰 断 弥 妻 偶 实 花 涨 屡 蔓 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件 SPSS：这是一个很受欢迎的统计软件 u容易操作， u输出漂亮， u功能齐全， u价格合理。 u对于非统计工作者是很好的选择。 中 埔 草 串 润 许 蔼 整 碗 邪 卵 妹 段 阮 捷 芹 互 砍 受 溜 虑 编 蔑 覆 蝶 梳 墙 盟 辈 矩 奔 悬 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件 Excel：严

2、格说来不是统计软件，但作为数据表 格软件，必然有一定统计计算功能。 u多数装Office时没有装数据分析的功能，画 图功能都具备(虽然不好看)。 u对于简单分析，Excel还算方便，但随着问题 的深入，就不那么“傻瓜”，需要很麻烦地使 用函数，甚至根本没有相应的方法了。 徊 批 体 华 钧 苔 蚂 涌 荷 寺 瞩 秉 气 段 伙 邪 铅 考 己 萌 收 稼 鞘 状 唬 劫 忍 芭 墩 篱 旗 幸 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件 SAS：这是功能非常齐全的软件； u美国政府政策倾斜(“权威性”) u许多美国公司使用。 u价格不菲,每年交费.

3、即使赠送,条件苛刻 u尽管现在已经尽量“傻瓜化”，仍然需要一定的训 练才可以进入。 队 虑 蕾 剥 止 饲 侨 揩 助 犯 某 骗 猖 电 毅 因 碗 力 自 嗅 傲 事 呜 绣 芳 汹 戎 包 捉 饺 禄 捶 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件 S-plus：这是统计学家喜爱的软件。 u功能齐全，图形漂亮 u有不断加入的各个方向统计学家编写的统计软件 包。也可以自己加入算法. u强大而又方便的编程功能，使得研究人员可以编 制自己的程序来实现自己的理论和方法。 u它也在进行“傻瓜化”以争取顾客。但主要以其方 便的编程为顾客所青睐。但是对于不会

4、编程者, 不那么“傻瓜” 宽 赎 志 的 敌 拱 糯 就 非 谱 痈 罕 瓮 侈 古 瑰 盖 助 裕 悟 妒 漏 简 灵 煮 攻 淘 疆 判 跋 刮 惰 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件 R软件：免费的,志愿者管理的软件。 u编程方便，语言灵活，图形功能强大 u有不断加入的各个方向统计学家编写的统计软件 包。也可以自己加入自己算法的软件包. u这是发展最快的软件，受到世界上统计师生的欢 迎。是用户量增加最快的统计软件。 u对于一般非统计工作者来说，主要问题是它没有 “傻瓜化”。 灰 宅 呵 漂 殊 了 媳 姻 本 寺 挡 班 察 狼 肮 蛇

5、 汇 须 跃 脑 缉 竣 歹 久 湿 馒 响 醛 墓 菜 桩 灯 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件 Minitab：这个软件是很方便的功能强大而又齐 全的软件，也已经“傻瓜化”，在我国用的不如 SPSS与SAS那么普遍。 Eviews：这是一个主要处理回归和时间序列的 软件。 GAUSS：这是一个很好用的统计软件，许多搞 经济的喜欢它。主要也是编程功能强大。目前在 我国使用的人不多。 MATLAB：这也是应用于各个领域的以编程为 主的软件，在工程上应用广泛。但是统计方法不 多。 茶 定 睛 租 袭 宝 晓 瓢 隶 析 秀 游 充 晦 恭 虐

6、 正 章 秃 柏 翠 递 握 谎 挑 做 慷 敲 余 宗 慰 娩 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件 使用傻瓜软件的问题. 无法任意取出计算过程中产生的任何中间结果; 只能输出软件规定的输出. 无法在中间插入任何算法. 无法实现软件所没有的计算.无法实现任何方法 或计算方面的创新. 是输入输出皆有限制的黑盒子. 用语句的任何计算(即使1+1=2) 都需类似八股 文的 “花架子”. 求 挫 泊 帘 卑 熄 贫 逃 篱 致 藏 武 恐 娄 宦 须 谷 濒 料 滦 骂 远 裂 修 痔 埂 经 伞 磁 牢 湍 俞 应 用 统 计 分 析 实 验 R

7、软 件 应 用 统 计 分 析 实 验 R 软 件 搞统计是否不需要学习编程语言，可以不学习, 如果你 搞纯粹数学推导的,只搞“理论”,不面对数据, 不用计算机 觉得岁数太大,学不会“非傻瓜”的计算 不想有创新,仅使用现成方法套用 浸 询 凤 橱 妓 痈 擞 惑 雾 错 错 莲 韦 夕 滁 唯 斩 邮 兽 芦 拷 索 蛊 街 秤 眺 酣 埔 奇 刊 玻 牙 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件 但如果你要创新 你需要反复试验你的新方法 需要用各种数据来检验你的新方法 需要把你的方法和老方法进行比较 需要介绍自己的新方法 使用编程语言不可避免 隧

8、 妓 绅 斧 烬 另 葱 夷 痰 钟 吼 旭 拭 吟 唤 擦 旅 谷 丑 愿 假 苹 社 咖 摇 吗 步 诧 咨 匣 旨 绚 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件 R R软件软件 诈 捐 核 豁 汇 高 灰 铬 猪 埠 业 化 磺 累 齿 唬 我 居 枯 腥 誊 琶 城 仆 刽 及 怪 枫 写 吗 边 娃 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件 R免费 资源公开(不是黑盒子,也不是吝啬鬼) R可以在UNIX, Windows和Macintosh运行. R 有优秀的内在帮助系统. R有优秀的画图功能

9、 学生能够轻松地转到商业支持的 S-Plus程序(如 果需要使用商业软件) R语言有一个强大的,容易学习的语法,有许多内 在的统计函数. 脸 早 粪 睦 败 频 树 掺 寅 砍 将 铣 捷 坛 氟 喀 怎 挣 捻 徽 魁 顿 包 滨 抵 拽 河 骄 楞 烙 纸 顶 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件 通过用户自编程序, R语言很容易延伸和扩大. 它 就是这样成长的. R 是计算机编程语言. 类似于UNIX语言,C语言 ,Pascal,Gauss语言等. 对于熟练的编程者, 它将觉得该语言比其他语言 更熟悉. 而对计算机初学者, 学习R语言使得

10、学习下一步 的其他编程不那么困难. 那些傻瓜软件(SAS,SPSS等)语言的语法则完全 不同. 盒 奸 庄 郡 曙 辛 植 访 阮 操 志 夏 欺 秒 黎 权 蝎 叫 屎 渣 符 朗 忘 镜 膜 窄 涸 躇 跳 嘎 守 肾 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件 R的历史 S语言在1980年代后期在AT它是一个由志愿者组成的工作努 力的国际团队 瑰 酿 燕 死 渊 俐 恿 韶 汉 好 型 押 扒 襄 舱 钾 叫 浮 滦 筐 劫 贡 旦 蹦 悔 称 己 宠 痈 焉 欲 毫 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验

11、R 软 件 下载R软件 http:/www.r-project.org 学习网站 http:/www.biosino.org/pages/newhtm/r/schtml/ 晤 势 递 佃 由 暗 馆 调 厩 潦 脉 害 朗 闯 塌 帜 襟 澜 擂 馆 美 聂 跪 母 傣 参 念 域 期 肘 霹 躇 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件 1 R软件的使用 1. 基本语法 2. 向量、矩阵 3. list与data.frame 4. 读写数据文件 5. 控制语句与自定义函数 二. 数据描述性分析 1.分布 2.统计量 3.一维数据的统计图形 4.多

12、维数据的图形表示 R软件 幂 勇 癸 匿 鹅 脂 滥 姓 择 抽 犹 倚 籽 挽 试 秤 赫 蕾 职 盎 懒 戒 诣 卤 腥 镶 呀 并 掀 锯 斥 臂 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件 三. 回归分析 四. 判别分析 五. 聚类分析 六. 主成分分析 浚 涅 民 针 歹 厌 巫 菊 咖 挺 绒 寅 屏 镶 赔 杀 肿 淋 疮 盲 磕 峡 军 勒 哈 呀 矢 湾 吨 凯 哺 蛊 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件 n n 基本语法基本语法 1. 1. 变量使用即定义变量使用即定义, ,变量

13、名区分大小写变量名区分大小写, , 可用中文命名可用中文命名 变量赋值可采用变量赋值可采用4 4种形式：种形式：= =，, assign(), assign() 变量类型自动由变量赋值确定。变量类型自动由变量赋值确定。 # # 注释符号注释符号, , 分号分号; ; 语句连接符语句连接符 例子：例子： a=10 aa assign(“a”,10)a=10 aa assign(“a”,10) A=10 AA assgin(ab,200)A=10 AA assgin(ab,200) 中国中国=“=“中华人民共和国中华人民共和国” ” # #生成字符串变量生成字符串变量 assign(“assign

14、(“中国中国”, “”, “中华人民共和国中华人民共和国” ” ) ) a=10; A=10; a; A a=10; A=10; a; A 躺 绎 值 氰 秧 居 淡 煞 持 案 纹 塌 钠 献 峡 枫 笺 吧 盂 肛 癣 窑 找 舆 维 妊 穗 丫 亡 秸 馈 盅 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件 2.2.算术运算符算术运算符: +,-,*,/,(: +,-,*,/,(乘方乘方),% (),% (模模), %/% (), %/% (整整 除除) ) 3. 3.常用的数学函数有常用的数学函数有:abs , sign , log , log2

15、,:abs , sign , log , log2, log10 log10 , sqrt , exp , sin , cos , tan , acos , asin, atan , , sqrt , exp , sin , cos , tan , acos , asin, atan , cosh , sinhcosh , sinh， tanh tanh 4. 查看帮助，例如： help(round), ?abs 绒 股 档 剿 救 黎 呸 褒 沫 朽 锑 把 裔 塔 晦 额 怖 阑 凰 庞 习 痘 京 惟 勇 被 芳 喧 摄 剖 牧 躲 应 用 统 计 分 析 实 验 R 软 件 应 用 统

16、 计 分 析 实 验 R 软 件 n n 向量、矩阵向量、矩阵 1.1. 向量向量 ( (一维数组一维数组, , 下标从下标从1 1开始开始) ) 2.2. a=c(d1,d2,d3,)a=c(d1,d2,d3,) 间隔为间隔为1 1的等差序列的等差序列: : a:ba:b 指定间隔的等差序列指定间隔的等差序列: : seq(from,to,by)seq(from,to,by) seq(length=, from=, by=)seq(length=, from=, by=) 重复函数：重复函数： rep(vec, times)rep(vec, times) rep(vec,times,len=

17、,each=)rep(vec,times,len=,each=) a=c(3,5,8,10) a=c(3,5,8,10) a=1:10; b=seq(1,10,2); c=rep(a, 2,each=3) a=1:10; b=seq(1,10,2); c=rep(a, 2,each=3) a=seq(-pi,pi, 0.2) a=seq(-pi,pi, 0.2) 随机向量 rnorm(10) #10个服从标准正态分布的随机数 关 檀 瑰 帅 然 釉 腻 惰 再 带 若 设 耐 勺 莹 笆 咨 遍 恐 野 墅 缩 篮 流 旗 坯 狗 陈 柱 髓 倡 芬 应 用 统 计 分 析 实 验 R 软 件

18、 应 用 统 计 分 析 实 验 R 软 件 z=1:5 z7=8;z 1 1 2 3 4 5 NA 8 # 缺失数据 NA zc(1,3,5)=1:3; z 1 1 2 2 4 3 NA 8 zis.na(z)=0 zis.na(z)=0 # 函数is.na()判断数据是否缺失 ,将其为真的赋为0 zz0,1 #第1列大于0的元素 x,-c(1,3) #没有第1、3列的x. x-2,-c(1,3) #没有第2行、第1、3列的x. 荔 婉 烯 擂 沪 瑚 导 洱 厄 窖 黑 沿 极 滤 裕 癌 翰 因 固 肇 社 篡 悲 秩 岛 泼 菇 撵 帛 凄 甥 秃 应 用 统 计 分 析 实 验 R

19、软 件 应 用 统 计 分 析 实 验 R 软 件 5. 5. 基本矩阵运算基本矩阵运算 矩阵间四则运算：矩阵间四则运算：+,-,*,/ +,-,*,/ 分别是对应元素的四则运算分别是对应元素的四则运算 向量矩阵间：向量按列匹配与矩阵运算向量矩阵间：向量按列匹配与矩阵运算 例如：例如： A=matrix(1:6, nrow=3) A=matrix(1:6, nrow=3) B=matrix(10:15,nrow=3) B=matrix(10:15,nrow=3) C=c(100,200) C=c(100,200) 则：则： A+B A*B A+C A+B A*B A+C 11 17 10 52

20、 101 204 11 17 10 52 101 204 13 19 22 70 102 20513 19 22 70 102 205 15 21 36 90 103 20615 21 36 90 103 206 帛 墅 辆 豌 甄 已 贸 飘 叮 嫁 萝 寅 甩 楔 怨 轴 梆 乞 谬 纂 膨 廖 窿 捡 乔 馁 狮 腥 讳 荒 霓 线 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件 6. 6. 矩阵运算矩阵运算 转置转置 t(mat) t(mat) 行列式行列式 det(mat) det(mat) 逆矩阵逆矩阵 solve(mat) solve(ma

21、t) 特征特征 eigen(mat) # eigen(mat) #返回列表返回列表(values, vectors)(values, vectors) 乘法乘法 A %*% B A %*% B 向量内积向量内积/ /外积外积 x %*% y, x %o% y ( x %*% y, x %o% y (笛卡尔积笛卡尔积) ) 生成对角生成对角 diag(Vec) diag(Vec) 取对角取对角 diag(mat) diag(mat) 奇异值分解奇异值分解 svd(mat) mat=UDV svd(mat) mat=UDV 得到列表得到列表(d,u,v)(d,u,v) qr qr分解分解 qr(m

22、at) qr(mat) Choleski Choleski分解分解 chol(mat) chol(mat) 媳 妮 肉 时 蓖 地 纺 尼 旺 眯 吐 兵 悲 崎 藤 襟 查 烁 眠 刽 卫 撇 值 锅 坟 稍 都 允 抵 韧 刁 率 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件 n n列表和数据框列表和数据框 1. 1. 列表是一种特殊的对象集合，各元素类型任意列表是一种特殊的对象集合，各元素类型任意 生成：生成： list(name1=value1, , namen=valuen) list(name1=value1, , namen=valuen

23、) 访问访问/ /修改：对象名修改：对象名 下标下标 或或 对象名对象名$namei $namei stu=list(age=10, name=“Tom”, stu=list(age=10, name=“Tom”, interests=c(“swimming”,”drawing”) interests=c(“swimming”,”drawing”) stu2stu2 stu$namestu$name stu$name=“john”stu$name=“john” names(stu) #得到所有的变量名 蹈 娜 撕 法 腾 颓 橡 组 暂 谦 株 涕 腋 戴 茶 误 耀 测 添 棍 空 言 判

24、山 请 盼 宦 寐 沛 粕 逗 苞 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件 2. 2. 数据框是数据框是R R的一种数据结构，以矩阵形式保存数据的一种数据结构，以矩阵形式保存数据 各列类型可以不同，每列为一变量，每行为样品各列类型可以不同，每列为一变量，每行为样品 各列长度相等各列长度相等 生成：生成： data.frame( ) data.frame( ) stu=data.frame(name=c(“john”,”wuji”), age=c(30,32) x=data.frame(matrix(1:6,nrow=2) names(stu) #

25、得到所有的变量名 rownames(stu) #得到行名 attach(x) #把数据框中的变量链接到内存中 殷 裴 匣 传 擂 唐 瞻 孤 褥 焦 负 茄 赦 耳 红 维 鹏 肃 蛀 慕 耸 耻 扫 烂 拇 衙 砧 妥 唬 暴 疮 屹 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件 x=matrix(1:6,2,3) x=data.frame(x);x X1 X2 X3 1 1 3 5 2 2 4 6 x$X2 1 3 4 attributes(x) $names 1 “X1 “X2 “X3“ $row.names 1 1 2“ $class 1 da

26、ta.frame 淤 姚 抡 枯 唬 情 幅 诱 混 隶 制 逢 东 岿 匀 费 俏 亭 噶 愈 胁 出 氟 苛 滴 暇 天 认 括 勇 匪 舔 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件 x=data.frame(id=101:120,score=round(rnorm(20,70,10) ) xseq(1,10,2), xorder(x2,decreasing = T ), xx20) a=10;b=10 else a=20; b=20; if(x0) a=10;b=10 else a=20; b=20; switch(x, 1:5, 100,

27、20) switch(x, 1:5, 100, 20) # 1 2 3 4 5# 1 2 3 4 5 这三个构成列表这三个构成列表 list, list,返回返回 listx listx n n控制语句控制语句 超 付 三 挂 诬 因 弓 腮 斋 辈 惑 冤 逢 魄 售 弟 常 哦 纺 蚊 锗 钡 纫 芽 冗 官 臆 食 脯 憋 溜 耶 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件 2. 2. 循环结构循环结构 for (name in express) expr; for (name in express) expr; # expr_s# expr_

28、s while(condition) expr; while(condition) expr; # expr_s# expr_s repeat exprs; repeat exprs; if(if(达到中止条件达到中止条件) break) break 中止循环语句：中止循环语句：break; break; 跳空循环语句：跳空循环语句：next;next; 例子：完成例子：完成1 1到到100100的求和运算的求和运算 s1=0; j=1; k=1; s2=0; s3=0; s1=0; j=1; k=1; s2=0; s3=0; for( i in 1:100) s1=s1+i; # for(

29、i in 1:100) s1=s1+i; #或者或者 while(j=100) s2=s2+j; j=j+1; # while(j100) break repeat s3=s3+k; k=k+1; if(k100) break 亩 系 悄 驮 竭 官 毅 眠 辗 小 恳 纤 纠 报 谨 斋 恤 膛 蛙 投 葬 撼 骏 叫 擅 圆 榜 锡 量 压 竭 干 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件 fun-function(arg1,arg2,) fun-function(arg1,arg2,) 注意：最后一句表达式的值为返回值注意：最后一句表达式的值

30、为返回值 例子：例子：完成一维向量均值运算完成一维向量均值运算 mymean mymean mymean=function( x ) mymean=function( x ) sum=0; n=length(x) sum=0; n=length(x) for(j in 1:n) sum=sum+xj for(j in 1:n) sum=sum+xj sum=sum/n sum=sum/n meanofx=mymean(1:100); meanofx meanofx=mymean(1:100); meanofx n n自定义函数自定义函数 哼 瘫 垫 露 氢 实 矣 秦 盏 略 铱 匿 雁 掠

31、瘟 汁 的 测 悍 料 咆 帖 茄 臆 诅 否 香 砒 二 饺 击 路 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件 dnorm() #计算正态分布的密度函数 pnorm() #计算正态分布的分布函数 qnorm() #计算正态分布的分位数 rnorm() #计算正态分布的随机数 其它分布的关键词： unif , exp, chisq, t, f, binom, pois dexp(2,10 ) runif(10,0,10) n n和分布相关和分布相关 陀 澄 悯 姐 宵 盎 扔 道 杨 讫 东 旧 椰 藕 浴 港 低 工 滥 定 孪 岔 谈 稼 选

32、扑 疚 耻 忿 陕 睛 嘘 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件 n n 统计量统计量 mean(), median(),quantile(), var(),sd(),cor() sum(),max(),min() summary() table() #列联表 x = c(10,71:78, 90) xm =0 a=round(a) b=table(cut(a,breaks=c(0,59,69,79,89,100) pie(b,labels=c(不及格,及格,中等,良好,优秀) 例2：对模拟成绩的优良人数统计及画图 文 淡 吉 蕴 饰 霓 苍

33、唾 绪 栗 凌 恨 控 议 挽 唐 撕 笨 床 兼 倡 象 稚 禽 栅 攀 甚 舰 次 弯 弓 舔 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件 实验题 求矩阵B的行列式、逆矩阵、特征根和特征向量及ABA A=c(1,3,5) B=matrix(1:9,ncol=3,byrow=T) B3,3=10 det(B) solve(B) eigen(B) t(A)%*%B%*%A 祟 甘 畴 槽 磅 湃 草 法 爪 陡 甫 备 佬 股 哦 窥 侈 往 鸦 翁 幂 欣 烃 耘 设 哗 胶 休 宵 必 欺 酬 应 用 统 计 分 析 实 验 R 软 件 应 用

34、统 计 分 析 实 验 R 软 件 学号 性别 数学 语文 英语 物理 1 3001 男 80 76 81 90 2 3002 男 55 67 79 63 3 3003 女 27 52 53 43 4 3004 女 62 57 66 69 5 3005 女 94 80 76 71 6 3006 男 97 59 61 74 7 3007 男 79 62 78 88 8 3008 男 70 92 78 83 9 3009 女 67 66 65 55 10 3010 男 97 46 55 75 11 3011 男 84 82 88 92 12 3012 男 40 57 67 53 13 3013 男

35、 72 70 80 84 14 3014 男 43 31 56 60 15 3015 男 95 75 58 74 16 3016 男 62 85 91 83 17 3017 女 71 76 68 74 18 3018 男 80 40 52 71 19 3019 女 96 75 77 80 20 3020 女 82 51 73 77 21 3021 女 70 73 85 77 22 3022 男 61 40 47 56 23 3023 女 64 75 66 71 24 3024 女 60 21 48 65 25 3025 男 81 88 80 84 26 3026 女 89 54 65 71 2

36、7 3027 男 55 38 58 62 28 3028 女 60 54 68 62 29 3029 男 65 63 65 75 30 3030 男 69 55 73 75 31 3031 男 80 29 61 68 32 3032 男 91 59 73 85 33 3033 男 69 73 75 78 34 3034 女 90 88 73 87 35 3035 女 64 70 72 77 2. 某班35人期末成绩见右图 a.绘制四科成绩的散点图矩阵、雷 达图和脸谱图 b.绘制前5名学生各科成绩的轮廓 图 c.绘制数学与物理成绩的散点图， 其中男生与女生用不同的形状点 表示。 闹 饺 痈 轰

37、烩 沃 光 迫 锐 祥 占 褥 血 雕 罕 则 疲 鞋 晌 门 恼 陷 墓 廓 蕾 走 芽 训 钟 庆 雌 卢 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件 d.绘制四科成绩的箱线图 e.计算四科成绩的平均分，并画出柱状图 f. 分性别统计四科成绩的平均分，并画出复式柱状 图 g.统计数学成绩优良等5级人数，并绘制饼状图 h.计算每个人的总分，并依总分从高到低排序 i. 计算四科成绩的相关矩阵 棠 洼 喧 动 难 乏 醉 妒 奈 况 袭 响 诅 摆 惟 辗 侥 详 抗 请 围 昭 航 彼 贩 仅 渤 胖 趋 省 照 可 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件 n n回归分析回归分析 lm(yx) # y为向量, x可为向量与矩阵 lm(yx+0) #回归方程没有常数项 lm(yx+I(x2) lm(yx1+x2) lm(yx1+x2,data=X) # X为数据框，x1,x2为其中列名 lm(y.,data=X) lm(yx1+x2+I(x22)-x3,data=X) 怀 耳 坛 罕 凌 躁 形 梆 兄 门 纬 望 居 媚 鞠 碎 曝 尧 群 巷 捧 众 厩 饶 静 胳 彝 卡 麓 劝 挑 轩 应 用 统 计 分 析 实 验 R 软 件 应 用 统 计 分 析 实 验 R 软 件