R_语言-统计.ppt

Stat.R R 语言与统计分析 目录 R 简介 R 语法与数据结构 程序控制语句 矩阵运算 统计函数 输入与输出 图形 1.3 统计分析软件包 CRAN提供了许多便于统计分析的宏 /src/contrib/PACKAGES.html VaR – 风险值分析 tseries – 时间序列分析 matrix – 矩阵运算 cinterface – C与R的接口 foreign – 读写由S, Minitab, SAS, SPSS, Stata等软件的数据 normix – 混合正态分布分析 nortest – 正态分布的Anderson-Darling检验 MCMCpack – 基于Gibbs抽样的MCMC抽样方法 还有很多…… 2. R语法与数据结构 2.1 语法 符号 命令或运算提示符 + 续行符 基本算术运算 + 加号 - 减号 * 乘号 / 除号 ^ 乘方 赋值符 = 或 - 建立向量的方法 seq() 向量(序列)具有较为简单的规律 rep() 向量(序列)具有较为复杂的规律 c() 向量(序列)没有什么规律 例： 1:10 seq(1,10,by=0.5) 或者 seq(from=1,to=10,by=0.5) 或者 seq(1,10,length=21) rep(2:5,2) 重复第一个自变量(2:5)若干次 rep(2:5,rep(2,4)) 2 2 3 3 4 4 5 5 x=c(42,7,64,9) length(x) 注意向量运算中的循环法则(recycling rule) 1:2+1:4 1:4+1:7 向量的下标(index)与向量子集(元素)的提取 正的下标 提取向量中对应的元素 负的下标 去掉向量中对应的元素 逻辑运算 提出向量中元素的值满足条件的元素 注：R中向量的下标从1开始，这与通常的统计或数学软件 一致而象C语言等计算机高级语言的向量下标则从0开始！ 例子： x = c(42,7,64,9) x[1] x[-2] x扣掉第2个位置后的向量 x[c(1,4)] x[which(x9)] 2.3 因子(factor) 统计中常处理的一类数据：分类数据(categorical data)；涉及的变量称作：名义(nominal)变量或分类变量。 R中用factor来表示分类变量；对于factor类型的数据，常用的函数有： table(); ftable(); 产生列联表(contingency table) chisq.test(); 对列联表做卡方检验 2.4 数据框(data frame) 一个数据框就是将许多向量组合起来的一个对象，它是二维的，通常其列表示变量，其行表示观测 数据框的用途 数据框的主要用途是保存统计建模需要的数据。 数据框的生成 例子： d - data.frame(name=c(李明, 张聪, 王建), age=c(30, 35, 28), height=c(180, 162, 175)) 数据框的读取 若数据本身保存在一个文件中，则可以使用 read.table() 仅接受带有分界符的ASCII数据 如果数据是电子报表的形式，则采用下面的两种变型 read.csv() 先将数据另存为带逗号的数据(Comma Seperated values) read.delim() 先将数据另存为用tab作为分界符的数据 2.5 列表(list) 有时需要生成包含不同类型的对象 R的列表(list) 是一种特别的对象集合，它的元素也由下标区分，但可以是任意类型。元素本身也可以是其它数据类型。 列表元素的引用： 列表名[[下标]]； 列表名[[“元素名”]] 或者 列表名$元素名 注意和 “列表名[下标]” 的区别：它取得的是子列表，而不是元素。 例子: foo = list(x = 1:6, y = matrix(1:4, nrow = 2)) foo $x [1] 1 2 3 4 5 6 $y [,1] [,2] [1,] 1 3 [2,] 2 4 列表子集的提取 提取一个子对象如foo的x, foo$xfoo[1] 3.3 函数 函数是一系列语句的组合，在R中可以写出自己的函数 形式: 变量名 = function( 变量列表 ) {函数体} 函数引用: 变量名(变量的值) 函数可以递归引用，但不提倡！ 例子 – 使用gamma函数求n!factorial = function(n)