应用统计分析实验r软件

应用统计分析实验 R 软件SPSS： 这是一个很受欢迎的统计软件u容易操作，u输出漂亮，u功能齐全，u价格合理。u对于非统计工作者是很好的选择。Excel： 严格说来不是统计软件，但作为数据表格软件，必然有一定统计计算功能。u多数装 Office时没有装数据分析的功能，画图功能都具备 (虽然不好看 )。u对于简单分析， Excel还算方便，但随着问题的深入，就不那么 “ 傻瓜 ” ，需要很麻烦地使用函数，甚至根本没有相应的方法了。SAS： 这是功能 非常 齐全的软件；u美国政府政策倾斜 (“权威性 ”)u许多美国公司使用。u价格不菲 ,每年交费 .即使赠送 ,条件苛刻u尽管现在已经尽量 “傻瓜化 ”，仍然需要一定的训练才可以进入。S-plus： 这是统计学家喜爱的软件。u功能齐全，图形漂亮u有不断加入的各个方向统计学家编写的统计软件包。也可以自己加入算法 .u强大而又方便的编程功能，使得研究人员可以编制自己的程序来实现自己的理论和方法。u它也在进行 “傻瓜化 ”以争取顾客。但主要以其方便的编程为顾客所青睐。但是对于不会编程者 ,不那么 “傻瓜 ”R软件 ：免费的 ,志愿者管理的软件。u编程方便，语言灵活，图形功能强大u有不断加入的各个方向统计学家编写的统计软件包。也可以自己加入自己算法的软件包 .u这是发展最快的软件，受到世界上统计师生的欢迎。是用户量增加最快的统计软件。u对于一般非统计工作者来说，主要问题是它没有“傻瓜化 ”。 Minitab： 这个软件是很方便的功能强大而又齐全的软件，也已经 “傻瓜化 ”，在我国用的不如SPSS与 SAS那么普遍。 Eviews： 这是一个主要处理回归和时间序列的软件。 GAUSS： 这是一个很好用的统计软件，许多搞经济的喜欢它。主要也是编程功能强大。目前在我国使用的人不多。 MATLAB： 这也是应用于各个领域的以编程为主的软件，在工程上应用广泛。但是统计方法不多。使用傻瓜软件的问题 无法任意取出计算过程中产生的任何中间结果 ;只能输出软件规定的输出 . 无法在中间插入任何算法 . 无法实现软件所没有的计算 .无法实现任何方法或计算方面的创新 . 是输入输出皆有限制的黑盒子 . 用语句的任何计算 (即使 1+1=2) 都需类似八股文的 “ 花架子 ” .搞统计是否不需要学习编程语言，可以不学习 ,如果你 搞纯粹数学推导的 ,只搞 “ 理论 ” ,不面对数据 , 不用计算机 觉得岁数太大 ,学不会 “ 非傻瓜 ” 的计算 不想有创新 ,仅使用现成方法套用但 如果你要创新 你需要反复试验你的新方法 需要用各种数据来检验你的新方法 需要把你的方法和老方法进行比较 需要介绍自己的新方法 使用编程语言不可避免R软件软件 R免费 资源公开 (不是黑盒子 ,也不是吝啬鬼 ) R可以在 UNIX, Windows和 Macintosh运行 . R 有 优秀的 内在帮助系统 . R有优秀的 画图功能 学生能够轻松地 转 到商业支持的 S-Plus程序 (如果需要使用商业软件 ) R语言有一个 强大的 ,容易学习 的语法 ,有许多内在的统计函数 . 通过用户自编程序 , R语言很容易延伸和扩大 . 它就是这样成长的 . R 是 计算机编程语言 . 类似于 UNIX语言 ,C语言,Pascal,Gauss语言等 . 对于熟练的编程者 , 它将觉得该语言比其他语言更熟悉 . 而对计算机初学者 , 学习 R语言使得学习下一步的其他编程不那么困难 . 那些傻瓜软件 (SAS,SPSS等 )语言的语法则完全不同 .R的历史 S语言在 1980年代后期在 AT它是一个由志愿者组成的工作努力的国际团队下载 R软件学习网站/pages/newhtm/r/schtml/1 R软件的使用1. 基本语法2. 向量、矩阵 3. list与 data.frame4. 读写数据文件5. 控制语句与自定义函数二 . 数据描述性分析1.分布2.统计量3.一维数据的统计图形4.多维数据的图形表示R软件三 . 回归分析四 . 判别分析五 . 聚类分析六 . 主成分分析n 基本语法基本语法1. 变量使用即定义变量使用即定义 ,变量名区分大小写变量名区分大小写 , 可用中文命名可用中文命名变量赋值可采用变量赋值可采用 4种形式：种形式： =， , assign()变量类型自动由变量赋值确定。变量类型自动由变量赋值确定。# 注释符号注释符号 , 分号分号 ; 语句连接符语句连接符例子：例子：a=10 aa assign(“a”,10)A=10 AA assgin(ab,200)中国中国 =“中华人民共和国中华人民共和国 ” #生成字符串变量生成字符串变量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开始开始 ) a=c(d1,d2,d3,)间隔为间隔为 1的等差序列的等差序列 : a:b指定间隔的等差序列指定间隔的等差序列 : seq(from,to,by) seq(length=, from=, by=) 重复函数：重复函数： rep(vec, times) rep(vec,times,len=,each=) a=c(3,5,8,10)a=1:10; b=seq(1,10,2); c=rep(a, 2,each=3)a=seq(-pi,pi, 0.2)随机向量 rnorm(10) #10个服从标准正态分布的随机数z=1:5z7=8;z1 1 2 3 4 5 NA 8 # 缺失数据 NAzc(1,3,5)=1:3; z1 1 2 2 4 3 NA 8 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.5. 基本矩阵运算基本矩阵运算矩阵间四则运算：矩阵间四则运算： +,-,*,/ 分别是对应元素的四则运算分别是对应元素的四则运算向量矩阵间：向量按列匹配与矩阵运算向量矩阵间：向量按列匹配与矩阵运算例如：例如：A=matrix(1:6, nrow=3)B=matrix(10:15,nrow=3)C=c(100,200)则：则： A+B A*B A+C11 17 10 52 101 204 13 19 22 70 102 20515 21 36 90 103 2066. 矩阵运算矩阵运算转置转置 t(mat)行列式行列式 det(mat)逆矩阵逆矩阵 solve(mat)特征特征 eigen(mat) #返回列表返回列表 (values, vectors)乘法乘法 A %*% B向量内积向量内积 /外积外积 x %*% y, x %o% y (笛卡尔积笛卡尔积 )生成对角生成对角 diag(Vec) 取对角取对角 diag(mat)奇异值分解奇异值分解 svd(mat) mat=UDV 得到列表得到列表 (d,u,v)qr分解分解 qr(mat)Choleski分解分解 chol(mat)n列表和数据框列表和数据框1. 列表是一种特殊的对象集合，各元素类型任意列表是一种特殊的对象集合，各元素类型任意生成：生成： list(name1=value1, , namen=valuen)访问访问 /修改：对象名修改：对象名 下标下标 或或 对象名对象名 $namei stu=list(age=10, name=“Tom”, interests=c(“swimming”,”drawing”)stu2stu$namestu$name=“john”names(stu) #得到所有的变量名2. 数据框是数据框是 R的一种数据结构，以矩阵形式保存数据的一种数据结构，以矩阵形式保存数据各列类型可以不同，每列为一变量，每行为样品各列类型可以不同，每列为一变