R语言统计分析简介.pptVIP

残差检验 正态性检验 程序命令： shapiro.test(lm1$res) t.test(lm1$res) 残差序列满足εt ~ N(0, σ2) 序列相关检验（D.W.） 程序命令： library(car) durbin.watson(lm1) 残差检验 自相关系数（ACF） 程序命令： acf(lm1$res) 偏自相关系数（PACF） 程序命令： pacf(lm1$res) 残差序列满足0均值、同方差、无自相关，即白噪声序列 聚类分析 根据事物本身的特征研究个体分类的方法 函数名：hclust(d, method =ward) 数据文件名：beer.txt 分析方法：欧氏距离、ward法、系统聚类 程序命令： beer=read.table(beer.txt, header=T) s(beer)=beer$beername hclust(dist(beer[, -1]), method=ward)$merge plot(hclust(dist(beer[, -1]), method=ward)) 聚类的凝聚过程 聚类分析R 判别分析 根据观察或测量到若干变量值，判断研究对象属于哪一类 函数名：lda(x, grouping) 数据文件名：disc.sav 分析方法：距离判别法 library(MASS) library(foreign) disc=read.spss(disc.sav, to.data.frame=T) lda1=lda(disc[,-1], disc$GROUP) table(predict(lda1, disc[,-1])$class, disc$GROUP) plot(predict(lda1, disc[,-1])$x, col=as.numeric(disc$GROUP)) 判别分析R与SPSS对比 R输出结果 SPSS输出结果 主成分分析 利用降维的思想，把多个指标转换成较少的几个互不相关的综合指标的统计方法 函数名：prcomp(x, scale=F) 数据文件名：company.txt 分析方法：从相关矩阵出发 程序命令： company=read.table(company.txt, header=T) prc=prcomp(company[,-1], scale=T) summary(prc) prc$rotation prc$x[,1:2] 主成分分析R与SPSS比较 R输出结果 SPSS输出结果 因子分析 探讨存在相关关系的变量之间，是否存在不能直接观察到但对可观测变量变化起支配作用的潜在因子的分析方法 函数名：factanal(x, factors, rotation=varimax) 数据文件名：student.sav 分析方法：极大似然法、方差最大正交旋转、回归法 程序命令： library(foreign) student=read.spss(student.sav, to.data.frame=T) factanal(student, 2, rotation=varimax , scores=regression) 因子分析R与SPSS对比 R使用mle方法估计参数，默认因子旋转方法为方差最大化 R输出结果 SPSS输出结果 对应分析 在一个低维度空间中对列联表里的两个名义变量间的关系进行描述 函数名：corresp(xtabs) 数据文件：caith（在MASS程序包中） 程序命令： library(MASS) caith plot(corresp(caith, nf=2)) 对应分析R caith数据文件格式： 对应分析行和列得分 R语言统计分析 易 楠 张青青 R语言简介 集数据分析与图形显示于一体的统计编程软件。 由Auckland大学统计系的Robert Gentleman和Ross Ihaka于1995年创立的。 免费下载地址 其它统计软件 SAS：世界上最权威的统计软件，价格昂贵，支持编程，统计学专业专用软件 SPSS：操作容易、输出漂亮、功能齐全、价格合理，非统计学专业首选软件 Eviews：具有强大的多元回归和时间序列分析功能，计量经济学专业专用软件 Matlab：功能强大的编程软件，统计分析功能较少，工程专业首选软件 Excel：具有简单的统计分析功能，商务办公首选软件 通病： “黑匣子” “傻的可爱” R语言特点 程序命令完全公开（非黑匣子） 计算过程保留在对象中，可随时调用 有不断加入的各个方向统计学家编写的统计软件包 编程语言极其简单（归纳为“三个一”） 一小内时掌握基本运算和绘图语句 一天内学会编写函数和利用软件包 一年内