最新应用统计分析实验R软件PPT课件.ppt

应用统计分析实验应用统计分析实验R R软件软件SPSS：这是一个很受欢迎的统计软件这是一个很受欢迎的统计软件u容易操作，容易操作，u输出漂亮，输出漂亮，u功能齐全，功能齐全，u价格合理。价格合理。u对于非统计工作者是很好的选择。对于非统计工作者是很好的选择。搞统计是否不需要学习编程语言，可以不学习搞统计是否不需要学习编程语言，可以不学习,如果你如果你搞纯粹数学推导的搞纯粹数学推导的,只搞只搞“理论理论”,”,不面对数据不面对数据,不用计算机不用计算机觉得岁数太大觉得岁数太大,学不会学不会“非傻瓜非傻瓜”的计算的计算不想有创新不想有创新,仅使用现成方法套用仅使用现成方法套用但如果你要创新但如果你要创新你需要反复试验你的新方法你需要反复试验你的新方法需要用各种数据来检验你的新方法需要用各种数据来检验你的新方法需要把你的方法和老方法进行比较需要把你的方法和老方法进行比较需要介绍自己的新方法需要介绍自己的新方法使用编程语言不可避免使用编程语言不可避免R软件R免费免费资源公开资源公开(不是黑盒子不是黑盒子,也不是吝啬鬼也不是吝啬鬼)R可以在可以在UNIX,Windows和和Macintosh运行运行.R有优秀的有优秀的内在帮助系统内在帮助系统.R有优秀的有优秀的画图功能画图功能学生能够轻松地学生能够轻松地转转到商业支持的到商业支持的S-Plus程序程序(如如果需要使用商业软件果需要使用商业软件)R语言有一个语言有一个强大的强大的,容易学习容易学习的语法的语法,有许多内有许多内在的统计函数在的统计函数.通过用户自编程序通过用户自编程序,R语言很容易延伸和扩大语言很容易延伸和扩大.它它就是这样成长的就是这样成长的.R是计算机编程语言是计算机编程语言.类似于类似于UNIX语言语言,C语言语言,Pascal,Gauss语言等语言等.对于熟练的编程者对于熟练的编程者,它将觉得该语言比其他语言它将觉得该语言比其他语言更熟悉更熟悉.而对计算机初学者而对计算机初学者,学习学习R语言使得学习下一步语言使得学习下一步的其他编程不那么困难的其他编程不那么困难.那些傻瓜软件那些傻瓜软件(SAS,SPSS等等)语言的语法则完全语言的语法则完全不同不同.R的历史的历史S语言在语言在1980年代后期在年代后期在AT&T实验室开发实验室开发.R项目由项目由Auckland大学统计系的大学统计系的RobertGentleman和和RossIhaka于于1995年开始的年开始的.它很快得到广泛用户的欢迎它很快得到广泛用户的欢迎.目前它是由目前它是由R核心核心发展团队维持发展团队维持;它是一个由志愿者组成的工作努它是一个由志愿者组成的工作努力的国际团队力的国际团队下载下载R软件软件http:/www.r-project.org学习网站http:/www.biosino.org/pages/newhtm/r/schtml/1R软件的使用软件的使用1.基本语法基本语法2.向量、矩阵向量、矩阵3.list与与data.frame4.读写数据文件读写数据文件5.控制语句与自定义函数控制语句与自定义函数二二.数据描述性分析数据描述性分析1.分布分布2.统计量统计量3.一维数据的统计图形一维数据的统计图形4.多维数据的图形表示多维数据的图形表示R软件三三.回归分析回归分析四四.判别分析判别分析五五.聚类分析聚类分析六六.主成分分析主成分分析n基本语法基本语法1.1.变量使用即定义变量使用即定义变量使用即定义变量使用即定义,变量名区分大小写变量名区分大小写变量名区分大小写变量名区分大小写,可用中文命名可用中文命名可用中文命名可用中文命名变量赋值可采用变量赋值可采用变量赋值可采用变量赋值可采用4 4种形式：种形式：种形式：种形式：=，,assign(),assign()变量类型自动由变量赋值确定。变量类型自动由变量赋值确定。变量类型自动由变量赋值确定。变量类型自动由变量赋值确定。#注释符号注释符号注释符号注释符号,分号分号分号分号;语句连接符语句连接符语句连接符语句连接符例子：例子：例子：例子：a=10aaassign(“a”,10)a=10aaassign(“a”,10)A=10AAassgin(ab,200)A=10AAassgin(ab,200)中国中国中国中国=“=“中华人民共和国中华人民共和国中华人民共和国中华人民共和国”#生成字符串变量生成字符串变量生成字符串变量生成字符串变量assign(“assign(“中国中国中国中国”,“”,“中华人民共和国中华人民共和国中华人民共和国中华人民共和国”)a=10;A=10;a;A2.算术运算符算术运算符:+,-,*,/,(乘方乘方),%(模模),%/%(整整除除)3.常用的数学函数有常用的数学函数有:abs,sign,log,log2,log10,sqrt,exp,sin,cos,tan,acos,asin,atan,cosh,sinh，tanh4.查看帮助，例如：help(round),?absn向量、矩阵向量、矩阵1.1.向量向量向量向量(一维数组一维数组一维数组一维数组,下标从下标从下标从下标从1 1开始开始开始开始)2.2.a=c(d1,d2,d3,)a=c(d1,d2,d3,)间隔为间隔为间隔为间隔为1 1的等差序列的等差序列的等差序列的等差序列:a:ba:b指定间隔的等差序列指定间隔的等差序列指定间隔的等差序列指定间隔的等差序列:seq(from,to,byseq(from,to,by)seq(lengthseq(length=,from=,by=)=,from=,by=)重复函数：重复函数：重复函数：重复函数：rep(vecrep(vec,times),times)rep(vec,times,lenrep(vec,times,len=,each=)=,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个服从标准正态分布的随机数个服从标准正态分布的随机数z=1:5z7=8;z112345NA8#缺失数据缺失数据NAzc(1,3,5)=1:3;z112243NA8zis.na(z)=0zis.na(z)=0#函数函数is.na()判断数据是否缺失判断数据是否缺失,将其为真的赋为将其为真的赋为0zz0,1#第第1列大于列大于0的元素的元素x,-c(1,3)#没有第没有第1、3列的列的x.x-2,-c(1,3)#没有第没有第2行、第行、第1、3列的列的x.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+BA*BA+CA+BA*BA+C1117105210120413192270102205152136901032066.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(笛卡尔积笛卡尔积笛卡尔积笛卡尔积)生成对角生成对角生成对角生成对角diag(Vec)diag(Vec)取对角取对角取对角取对角diag(mat)diag(mat)奇异值分解奇异值分解奇异值分解奇异值分解svd(mat)mat=UDVsvd(mat)mat=UDV 得到列表得到列表得到列表得到列表(d,u,v)(d,u,v)qrqr分解分解分解分解qr(mat)qr(mat)CholeskiCholeski分解分解分解分解chol(mat)chol(mat)n列表和数据框列表和数据框1.1.列表是一种特殊的对象集合，各元素类型任意列表是一种特殊的对象集合，各元素类型任意列表是一种特殊的对象集合，各元素类型任意列表是一种特殊的对象集合，各元素类型任意生成：生成：生成：生成：list(name1=value1,namen=valuen)list(name1=value1,namen=valuen)访问访问访问访问/修改：对象名修改：对象名修改：对象名修改：对象名 下标下标下标下标或或或或 对象名对象名对象名对象名$namei$nameistu=list(age=10,name=“Tom”,stu=list(age=10,name=“Tom”,interests=c(“swimming”,”drawing”)interests=c(“swimming”,”drawing”)stu2stu2stu$namestu$namestu$name=“john”stu$name=“john”names(stu)#得到所有的变量名得到所有的变量名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)#得到所有的变量名得到所有的变量名rownames(stu)#得到行名得到行名attach(x)#把数据框中的变量链接到内存中把数据框中的变量链接到内存中x=matrix(1:6,2,3)x=data.frame(x);xX1X2X311352246x$X2134attributes(x)$names1“X1“X2“X3“$row.names112“$class1data.framex=data.frame(id=101:120,score=round(rnorm(20,70,10)xseq(1,10,2),xorder(x2,decreasing=T),xx20)a=10;b=10elsea=20;b=20;if(x0)a=10;b=10elsea=20;b=20;switch(x,1:5,100,20)switch(x,1:5,100,20)#12345#12345这三个构成列表这三个构成列表这三个构成列表这三个构成列表list,list,返回返回返回返回listxlistxn控制语句控制语句2.2.循环结构循环结构循环结构循环结构 for(nameinexpress)expr;for(nameinexpress)expr;#expr_s#expr_swhile(condition)expr;while(condition)expr;#expr_s#expr_srepeatexprs;repeatexprs;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(iin1:100)s1=s1+i;#for(iin1:100)s1=s1+i;#或者或者或者或者 while(j=100)s2=s2+j;j=j+1;#while(j100)breakrepeats3=s3+k;k=k+1;if(k100)breakfun-function(arg1,arg2,)fun-function(arg1,arg2,)注意：最后一句表达式的值为返回值注意：最后一句表达式的值为返回值注意：最后一句表达式的值为返回值注意：最后一句表达式的值为返回值 例子：例子：例子：例子：完成一维向量均值运算完成一维向量均值运算完成一维向量均值运算完成一维向量均值运算mymeanmymeanmymean=function(x)mymean=function(x)sum=0;n=length(x)sum=0;n=length(x)for(jin1:n)sum=sum+xjfor(jin1:n)sum=sum+xjsum=sum/nsum=sum/n meanofx=mymean(1:100);meanofxmeanofx=mymean(1:100);meanofxn自定义函数自定义函数dnorm()#计算正态分布的密度函数计算正态分布的密度函数pnorm()#计算正态分布的分布函数计算正态分布的分布函数qnorm()#计算正态分布的分位数计算正态分布的分位数rnorm()#计算正态分布的随机数计算正态分布的随机数其它分布的关键词：其它分布的关键词：unif,exp,chisq,t,f,binom,poisdexp(2,10)runif(10,0,10)n和分布相关和分布相关n统计量统计量mean(),median(),quantile(),var(),sd(),cor()sum(),max(),min()summary()table()#列联表列联表x=c(10,71:78,90)xm=0&a=0&a=100;a=round(a)b=table(cut(a,breaks=c(0,59,69,79,89,100)pie(b,labels=c(不及格,及格,中等,良好,优秀)例例2：对模拟成绩的优良人数统计及画图：对模拟成绩的优良人数统计及画图实验题求矩阵B的行列式、逆矩阵、特征根和特征向量及ABAA=c(1,3,5)B=matrix(1:9,ncol=3,byrow=T)B3,3=10det(B)solve(B)eigen(B)t(A)%*%B%*%A学号 性别 数学 语文 英语 物理1 3001 男 80 76 81 902 3002 男 55 67 79 633 3003 女 27 52 53 434 3004 女 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 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 3033 男 69 73 75 7834 3034 女 90 88 73 8735 3035 女 64 70 72 772.某班某班35人期末成绩见右图人期末成绩见右图a.绘制四科成绩的散点图矩阵、雷绘制四科成绩的散点图矩阵、雷达图和脸谱图达图和脸谱图b.绘制前绘制前5名学生各科成绩的轮廓名学生各科成绩的轮廓图图c.绘制数学与物理成绩的散点图，绘制数学与物理成绩的散点图，其中男生与女生用不同的形状点其中男生与女生用不同的形状点表示。表示。d.绘制四科成绩的箱线图绘制四科成绩的箱线图e.计算四科成绩的平均分，并画出柱状图计算四科成绩的平均分，并画出柱状图f.分性别统计四科成绩的平均分，并画出复式柱状分性别统计四科成绩的平均分，并画出复式柱状图图g.统计数学成绩优良等统计数学成绩优良等5级人数，并绘制饼状图级人数，并绘制饼状图h.计算每个人的总分，并依总分从高到低排序计算每个人的总分，并依总分从高到低排序i.计算四科成绩的相关矩阵计算四科成绩的相关矩阵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)结束语结束语谢谢大家聆听！谢谢大家聆听！53