《应用统计分析实验软件.ppt》由会员分享,可在线阅读,更多相关《应用统计分析实验软件.ppt(54页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、应用统计分析实验软件现在学习的是第1页,共54页SPSS:这是一个很受欢迎的统计软件这是一个很受欢迎的统计软件u容易操作,容易操作,u输出漂亮,输出漂亮,u功能齐全,功能齐全,u价格合理。价格合理。u对于非统计工作者是很好的选择。对于非统计工作者是很好的选择。现在学习的是第2页,共54页ExcelExcel:严格说来不是统计软件,但作为数据表:严格说来不是统计软件,但作为数据表格软件,必然有一定统计计算功能。格软件,必然有一定统计计算功能。u多数装多数装OfficeOffice时没有装数据分析的功能,画图时没有装数据分析的功能,画图功能都具备功能都具备(虽然不好看虽然不好看)。u对于简单分析,
2、对于简单分析,ExcelExcel还算方便,但随着问题还算方便,但随着问题的深入,就不那么的深入,就不那么“傻瓜傻瓜”,需要很麻烦地使,需要很麻烦地使用函数,甚至根本没有相应的方法了。用函数,甚至根本没有相应的方法了。现在学习的是第3页,共54页SAS:这是功能这是功能非常非常齐全的软件;齐全的软件;u美国政府政策倾斜美国政府政策倾斜(“权威性权威性”)u许多美国公司使用。许多美国公司使用。u价格不菲价格不菲,每年交费每年交费.即使赠送即使赠送,条件苛刻条件苛刻u尽管现在已经尽量尽管现在已经尽量“傻瓜化傻瓜化”,仍然需要一定的,仍然需要一定的训练才可以进入。训练才可以进入。现在学习的是第4页,
3、共54页S-plus:这是统计学家喜爱的软件。:这是统计学家喜爱的软件。u功能齐全,图形漂亮功能齐全,图形漂亮u有不断加入的各个方向统计学家编写的统计软件有不断加入的各个方向统计学家编写的统计软件包。也可以自己加入算法包。也可以自己加入算法.u强大而又方便的编程功能,使得研究人员可以编强大而又方便的编程功能,使得研究人员可以编制自己的程序来实现自己的理论和方法。制自己的程序来实现自己的理论和方法。u它也在进行它也在进行“傻瓜化傻瓜化”以争取顾客。但主要以其以争取顾客。但主要以其方便的编程为顾客所青睐。但是对于不会编程者方便的编程为顾客所青睐。但是对于不会编程者,不那么不那么“傻瓜傻瓜”现在学习
4、的是第5页,共54页R软件软件:免费的:免费的,志愿者管理的软件。志愿者管理的软件。u编程方便,语言灵活,图形功能强大编程方便,语言灵活,图形功能强大u有不断加入的各个方向统计学家编写的统计软件有不断加入的各个方向统计学家编写的统计软件包。也可以自己加入自己算法的软件包包。也可以自己加入自己算法的软件包.u这是发展最快的软件,受到世界上统计师生的欢这是发展最快的软件,受到世界上统计师生的欢迎。是用户量增加最快的统计软件。迎。是用户量增加最快的统计软件。u对于一般非统计工作者来说,主要问题是它没有对于一般非统计工作者来说,主要问题是它没有“傻瓜化傻瓜化”。现在学习的是第6页,共54页Minita
5、b:这个软件是很方便的功能强大而又齐:这个软件是很方便的功能强大而又齐全的软件,也已经全的软件,也已经“傻瓜化傻瓜化”,在我国用的不如,在我国用的不如SPSS与与SAS那么普遍。那么普遍。Eviews:这是一个主要处理回归和时间序列的:这是一个主要处理回归和时间序列的软件。软件。GAUSS:这是一个很好用的统计软件,许多搞:这是一个很好用的统计软件,许多搞经济的喜欢它。主要也是编程功能强大。目前在经济的喜欢它。主要也是编程功能强大。目前在我国使用的人不多。我国使用的人不多。MATLAB:这也是应用于各个领域的以编程为主这也是应用于各个领域的以编程为主的软件,在工程上应用广泛。但是统计方法不多。
6、的软件,在工程上应用广泛。但是统计方法不多。现在学习的是第7页,共54页使用傻瓜软件的问题使用傻瓜软件的问题.无法任意取出计算过程中产生的任何中间结果无法任意取出计算过程中产生的任何中间结果;只能输出软件规定的输出只能输出软件规定的输出.无法在中间插入任何算法无法在中间插入任何算法.无法实现软件所没有的计算无法实现软件所没有的计算.无法实现任何方法无法实现任何方法或计算方面的创新或计算方面的创新.是输入输出皆有限制的黑盒子是输入输出皆有限制的黑盒子.用语句的任何计算用语句的任何计算(即使即使1+1=2)1+1=2)都需类似八股都需类似八股文的文的“花架子花架子”.现在学习的是第8页,共54页搞
7、统计是否不需要学习编程语言,可以不学习搞统计是否不需要学习编程语言,可以不学习,如如果你果你搞纯粹数学推导的搞纯粹数学推导的,只搞只搞“理论理论”,不面对数据不面对数据,不用计算机不用计算机觉得岁数太大觉得岁数太大,学不会学不会“非傻瓜非傻瓜”的计算的计算不想有创新不想有创新,仅使用现成方法套用仅使用现成方法套用现在学习的是第9页,共54页但如果你要创新但如果你要创新你需要反复试验你的新方法你需要反复试验你的新方法需要用各种数据来检验你的新方法需要用各种数据来检验你的新方法需要把你的方法和老方法进行比较需要把你的方法和老方法进行比较需要介绍自己的新方法需要介绍自己的新方法使用编程语言不可避免使
8、用编程语言不可避免现在学习的是第10页,共54页R软件现在学习的是第11页,共54页R免费免费 资源公开资源公开(不是黑盒子不是黑盒子,也不是吝啬鬼也不是吝啬鬼)R可以在可以在UNIX,Windows和和Macintosh运行运行.R 有优秀的有优秀的内在帮助系统内在帮助系统.R有优秀的有优秀的画图功能画图功能学生能够轻松地学生能够轻松地转转到商业支持的到商业支持的 S-Plus程序程序(如如果需要使用商业软件果需要使用商业软件)R语言有一个语言有一个强大的强大的,容易学习容易学习的语法的语法,有许多内有许多内在的统计函数在的统计函数.现在学习的是第12页,共54页通过用户自编程序通过用户自编
9、程序,R语言很容易延伸和扩大语言很容易延伸和扩大.它它就是这样成长的就是这样成长的.R 是计算机编程语言是计算机编程语言.类似于类似于UNIX语言语言,C语言语言,Pascal,Gauss语言等语言等.对于熟练的编程者对于熟练的编程者,它将觉得该语言比其他语言它将觉得该语言比其他语言更熟悉更熟悉.而对计算机初学者而对计算机初学者,学习学习R语言使得学习下一步语言使得学习下一步的其他编程不那么困难的其他编程不那么困难.那些傻瓜软件那些傻瓜软件(SAS,SPSS等等)语言的语法则完全语言的语法则完全不同不同.现在学习的是第13页,共54页R的历史的历史S语言在语言在1980年代后期在年代后期在AT
10、&T实验室开发实验室开发.R 项目由项目由Auckland 大学统计系的大学统计系的Robert Gentleman和和Ross Ihaka于于1995年开始的年开始的.它很快得到广泛用户的欢迎它很快得到广泛用户的欢迎.目前它是由目前它是由R核心核心发展团队维持发展团队维持;它是一个由志愿者组成的工作努它是一个由志愿者组成的工作努力的国际团队力的国际团队现在学习的是第14页,共54页下载下载R软件软件http:/www.r-project.org学习网站http:/www.biosino.org/pages/newhtm/r/schtml/现在学习的是第15页,共54页1R软件的使用软件的使用
11、1.基本语法基本语法2.向量、矩阵向量、矩阵 3.list与与data.frame4.读写数据文件读写数据文件5.控制语句与自定义函数控制语句与自定义函数二二.数据描述性分析数据描述性分析1.分布分布2.统计量统计量3.一维数据的统计图形一维数据的统计图形4.多维数据的图形表示多维数据的图形表示R软件现在学习的是第16页,共54页三三.回归分析回归分析四四.判别分析判别分析五五.聚类分析聚类分析六六.主成分分析主成分分析现在学习的是第17页,共54页n 基本语法基本语法1.1.变量使用即定义变量使用即定义变量使用即定义变量使用即定义,变量名区分大小写变量名区分大小写变量名区分大小写变量名区分大
12、小写,可用中文命名可用中文命名可用中文命名可用中文命名 变量赋值可采用变量赋值可采用变量赋值可采用变量赋值可采用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)中国
13、中国中国中国=“=“中华人民共和国中华人民共和国中华人民共和国中华人民共和国”#生成字符串变量生成字符串变量生成字符串变量生成字符串变量assign(“assign(“中国中国中国中国”,“,“中华人民共和国中华人民共和国中华人民共和国中华人民共和国”)a=10;A=10;a;A现在学习的是第18页,共54页2.算术运算符算术运算符:+,-,*,/,(乘方乘方),%(模模),%/%(整除整除)3.常用的数学函数有常用的数学函数有:abs,sign,log,log2,log10,sqrt,exp,sin,cos,tan,acos,asin,atan,cosh,sinh,tanh4.查看帮助,例如
14、:help(round),?abs现在学习的是第19页,共54页n 向量、矩阵向量、矩阵 1.1.向量向量向量向量 (一维数组一维数组一维数组一维数组,下标从下标从下标从下标从1 1开始开始开始开始)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=)重复函数:重复函数:重复函数:重复函
15、数:rep(vec,times)rep(vec,times)rep(vec,times,len=,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个服从标准正态分布的随机数个服从标准正态分布的随机数现在学习的是第20页,共54页z=1:5z7=8;z1 1 2 3 4 5
16、 NA 8#缺失数据缺失数据 NAzc(1,3,5)=1:3;z1 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.现在学习的是第26页,共54页5.5.基本矩阵运算基本矩阵运算基本矩阵运算基本矩阵运算 矩阵间四则运算:矩阵间四则运算:矩阵间四则运算:矩阵间四则运算:+,-,*,/+,-,*,/分别是对应元素的四则运算
17、分别是对应元素的四则运算分别是对应元素的四则运算分别是对应元素的四则运算 向量矩阵间:向量按列匹配与矩阵运算向量矩阵间:向量按列匹配与矩阵运算向量矩阵间:向量按列匹配与矩阵运算向量矩阵间:向量按列匹配与矩阵运算例如:例如:例如:例如: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+CA+B A*B A+C11 17 10 52 101 204 13 19 22 70 102 20515 21
18、36 90 103 206现在学习的是第27页,共54页6.6.矩阵运算矩阵运算矩阵运算矩阵运算 转置转置转置转置 t(mat)t(mat)行列式行列式行列式行列式 det(mat)det(mat)逆矩阵逆矩阵逆矩阵逆矩阵 solve(mat)solve(mat)特征特征特征特征 eigen(mat)#eigen(mat)#返回列表返回列表返回列表返回列表(values,vectors)(values,vectors)乘法乘法乘法乘法 A%*%BA%*%B 向量内积向量内积向量内积向量内积/外积外积外积外积 x%*%y,x%o%y (x%*%y,x%o%y (笛卡尔积笛卡尔积笛卡尔积笛卡尔积)
19、生成对角生成对角生成对角生成对角 diag(Vec)diag(Vec)取对角取对角取对角取对角 diag(mat)diag(mat)奇异值分解奇异值分解奇异值分解奇异值分解 svd(mat)mat=UDVsvd(mat)mat=UDV 得到列表得到列表得到列表得到列表(d,u,v)(d,u,v)qr qr分解分解分解分解 qr(mat)qr(mat)Choleski Choleski分解分解分解分解 chol(mat)chol(mat)现在学习的是第28页,共54页n列表和数据框列表和数据框1.1.列表是一种特殊的对象集合,各元素类型任意列表是一种特殊的对象集合,各元素类型任意列表是一种特殊的
20、对象集合,各元素类型任意列表是一种特殊的对象集合,各元素类型任意生成:生成:生成:生成:list(name1=value1,namen=valuen)list(name1=value1,namen=valuen)访问访问访问访问/修改:对象名修改:对象名修改:对象名修改:对象名 下标下标下标下标 或或或或 对象名对象名对象名对象名$namei$namei stu=list(age=10,name=“Tom”,stu=list(age=10,name=“Tom”,interests=c(“swimming”,”drawing”)interests=c(“swimming”,”drawing”)s
21、tu2stu2stu$namestu$namestu$name=“john”stu$name=“john”names(stu)#得到所有的变量名得到所有的变量名现在学习的是第29页,共54页2.2.数据框是数据框是数据框是数据框是R R的一种数据结构,以矩阵形式保存数据各列类的一种数据结构,以矩阵形式保存数据各列类的一种数据结构,以矩阵形式保存数据各列类的一种数据结构,以矩阵形式保存数据各列类型可以不同,每列为一变量,每行为样品型可以不同,每列为一变量,每行为样品型可以不同,每列为一变量,每行为样品型可以不同,每列为一变量,每行为样品 各列长度相等各列长度相等各列长度相等各列长度相等生成:生成
22、:生成:生成: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)#得到所有的变量名得到所有的变量名 rownames(stu)#得到行名得到行名 attach(x)#把数据框中的变量链接到内存中把数据框中的变量链接到内存中 现在学习的是第30页,共54页x=matrix(1:6,2,3)x=data.frame(x);x X1 X2 X31 1 3 52 2 4 6x$X21 3 4attributes(x)$name
23、s1“X1“X2“X3“$row.names1 1 2“$class1 data.frame现在学习的是第31页,共54页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,20)switch(x,1:5,100,20)#1 2 3 4 5#1 2 3 4 5 这三个构成列表这三个构成列表这三个构成列表这三个构成列表
24、list,list,返回返回返回返回 listxlistxn控制语句控制语句现在学习的是第35页,共54页2.2.循环结构循环结构循环结构循环结构 for(name in express)expr;for(name in express)expr;#expr_s#expr_s while(condition)expr;while(condition)expr;#expr_s#expr_s repeat exprs;repeat exprs;if(if(达到中止条件达到中止条件达到中止条件达到中止条件)break)break 中止循环语句:中止循环语句:中止循环语句:中止循环语句:break;b
25、reak;跳空循环语句:跳空循环语句:跳空循环语句:跳空循环语句: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(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 现在学习的是第36页,共54页 fun-f
26、unction(arg1,arg2,)fun-function(arg1,arg2,)注意:最后一句表达式的值为返回值注意:最后一句表达式的值为返回值注意:最后一句表达式的值为返回值注意:最后一句表达式的值为返回值 例子:例子:例子:例子:完成一维向量均值运算完成一维向量均值运算完成一维向量均值运算完成一维向量均值运算 mymeanmymean 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/
27、n sum=sum/n meanofx=mymean(1:100);meanofx meanofx=mymean(1:100);meanofxn自定义函数自定义函数现在学习的是第37页,共54页dnorm()#计算正态分布的密度函数计算正态分布的密度函数pnorm()#计算正态分布的分布函数计算正态分布的分布函数qnorm()#计算正态分布的分位数计算正态分布的分位数rnorm()#计算正态分布的随机数计算正态分布的随机数其它分布的关键词:其它分布的关键词:unif ,exp,chisq,t,f,binom,poisdexp(2,10)runif(10,0,10)n和分布相关和分布相关现在学习
28、的是第38页,共54页n 统计量统计量mean(),median(),quantile(),var(),sd(),cor()sum(),max(),min()summary()table()#列联表列联表x=c(10,71:78,90)xm-mean(x)c(xm,mean(x,trim=0.10)A=cbind(c(90,60,70),c(85,66,95)apply(A,1,mean)#对每一行取平均apply(A,2,mean)#对每一列取平均var(A)cor(A)现在学习的是第39页,共54页n数据的图形表示数据的图形表示x=seq(-10,10,2);y=rnorm(11)plot
29、(x,y)plot(x,y,main=散点图散点图,xlab=横坐标横坐标x,ylab=纵坐标纵坐标y)text(x,y,1:length(x),-1)#写入点序号写入点序号,第三个参数默认如此第三个参数默认如此plot(x,y,type=o,col=red,pch=2,lty=3)1.散点图散点图 plot(x,y,)现在学习的是第40页,共54页type=“p”#散点图,默认散点图,默认type=“l#绘实线绘实线type=“o”#实线通过所有的点实线通过所有的点type=“n”#不画点不画点参数参数lty表示线的类型表示线的类型0=blank,1=solid(default),2=das
30、hed,3=dotted,4=dotdash,5=longdash,6=twodash 现在学习的是第41页,共54页低水平作图函数还有低水平作图函数还有points()lines()ablines()text()legend()现在学习的是第42页,共54页x=seq(-10,10,2);y=rnorm(11)par(mfrow=c(3,1)#准备画准备画3个图个图plot(x,y)plot(x,y,main=散点图散点图,xlab=横坐标横坐标x,ylab=纵坐标纵坐标y)text(x,y,1:length(x),-1)#写入点序号写入点序号plot(x,y,type=o,col=red,
31、pch=2,lty=3)现在学习的是第43页,共54页2.散点图矩阵散点图矩阵 pairs(x,)A=matrix(rnorm(12,70,20),nrow=4)pairs(A)现在学习的是第44页,共54页3.多组散点图多组散点图matplot(x,y,)x=seq(0,2*pi,0.2)matplot(x,cbind(sin(x),cos(x),type=c(“o”,”o”),col=2:3)现在学习的是第45页,共54页4.星座图星座图 5.脸谱图脸谱图library(aplpack)faces(x,nrow.plot,ncol.plot)stars(x,full=TRUE,draw.s
32、egments=FALSE)现在学习的是第46页,共54页例例1.在同一图中绘制标准正态分布和自由度为在同一图中绘制标准正态分布和自由度为5的的t分布的概率密度图分布的概率密度图.plot(x,dnorm(x),type=l,ylab=f(x),main=概率密度图)lines(x,dt(x,5),col=red,lty=2)legend(2,0.4,c(N(0,1),t(5),lty=1:2,col=c(black,red)matplot(x,cbind(dnorm(x),dt(x,5),type=l,ylab=f(x),main=概率密度图)legend(2,0.4,c(N(0,1),t(
33、5),lty=1:2,col=1:2)现在学习的是第47页,共54页lwl.andrews=function(X)t=seq(-pi,pi,0.2)n=nrow(X)f=matrix(0,nrow=length(t),ncol=n)for(i in 1:n)f,i=Xi,1/sqrt(2)for(j in 2:ncol(X)if(j%2=0)f,i=f,i+Xi,j*sin(j/2*t)else f,i=f,i+Xi,j*cos(j%/%2*t)matplot(t,f,col=1:n,type=o,main=调和曲线图调和曲线图)例2:绘制调和曲线图现在学习的是第48页,共54页hist(x,
34、freq)#hist(x,freq)#直方图,直方图,freq=counts,(freq=counts,(默认默认)为频数直方图为频数直方图,freq=TRUE ,freq=TRUE ,为频率直方图,为频率直方图,freq=FALSE,freq=FALSE,为密度为密度直方图直方图boxplot(x)#boxplot(x)#箱线图箱线图stem(x)#stem(x)#茎叶图茎叶图n其它图形其它图形a=rnorm(100,70,10);hist(a)b=rbinom(100,5,0.5)nums=table(b)barplot(nums)pie()#pie()#饼状图饼状图barplot(x,b
35、eside=F)#barplot(x,beside=F)#柱状图柱状图现在学习的是第49页,共54页例例1:31省市地区的消费支出数据,进行如下统计省市地区的消费支出数据,进行如下统计1.计算各地区的人均消费支出总和,并画出柱状图计算各地区的人均消费支出总和,并画出柱状图2.画出前十个地区的关于衣着和居住两项的复式柱状图。画出前十个地区的关于衣着和居住两项的复式柱状图。3.计算各消费类型的均值,并画出饼状图和柱状图计算各消费类型的均值,并画出饼状图和柱状图3.画出各消费类型的箱线图画出各消费类型的箱线图4.画出各地区的雷达图和脸谱图画出各地区的雷达图和脸谱图 barplot(apply(X,1
36、,sum),col=rainbow(31),cex.names=0.5)pie(apply(X,2,mean)barplot(apply(X,2,mean),names=NA,legend=names(X),col=1:ncol(X)barplot(t(X1:10,c(2,7),beside=T,legend=T)现在学习的是第50页,共54页实验题求矩阵B的行列式、逆矩阵、特征根和特征向量及ABA现在学习的是第51页,共54页学号 性别 数学 语文 英语 物理1 3001 男 80 76 81 902 3002 男 55 67 79 633 3003 女 27 52 53 434 3004
37、女 62 57 66 695 3005 女 94 80 76 716 3006 男 97 59 61 747 3007 男 79 62 78 888 3008 男 70 92 78 839 3009 女 67 66 65 5510 3010 男 97 46 55 7511 3011 男 84 82 88 9212 3012 男 40 57 67 5313 3013 男 72 70 80 8414 3014 男 43 31 56 6015 3015 男 95 75 58 7416 3016 男 62 85 91 8317 3017 女 71 76 68 7418 3018 男 80 40 52
38、7119 3019 女 96 75 77 8020 3020 女 82 51 73 7721 3021 女 70 73 85 7722 3022 男 61 40 47 5623 3023 女 64 75 66 7124 3024 女 60 21 48 6525 3025 男 81 88 80 8426 3026 女 89 54 65 7127 3027 男 55 38 58 6228 3028 女 60 54 68 6229 3029 男 65 63 65 7530 3030 男 69 55 73 7531 3031 男 80 29 61 6832 3032 男 91 59 73 8533 3
39、033 男 69 73 75 7834 3034 女 90 88 73 8735 3035 女 64 70 72 772.某班某班35人期末成绩见右图人期末成绩见右图a.绘制四科成绩的散点图矩阵、雷达图绘制四科成绩的散点图矩阵、雷达图和脸谱图和脸谱图b.绘制前绘制前5名学生各科成绩的轮廓图名学生各科成绩的轮廓图c.绘制数学与物理成绩的散点图,其中绘制数学与物理成绩的散点图,其中男生与女生用不同的形状点表示。男生与女生用不同的形状点表示。现在学习的是第52页,共54页d.绘制四科成绩的箱线图绘制四科成绩的箱线图e.计算四科成绩的平均分,并画出柱状图计算四科成绩的平均分,并画出柱状图f.分性别统计四科成绩的平均分,并画出复式柱状图分性别统计四科成绩的平均分,并画出复式柱状图g.在数据阵中增加总分列,并依总分从高到低排序在数据阵中增加总分列,并依总分从高到低排序h.计算四科成绩的相关矩阵计算四科成绩的相关矩阵现在学习的是第53页,共54页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)现在学习的是第54页,共54页
限制150内