r语言常用统计方法实现

r语言常用统计方法实现

ID:24848962

大小:150.00 KB

页数:26页

时间:2018-11-16

r语言常用统计方法实现_第1页
r语言常用统计方法实现_第2页
r语言常用统计方法实现_第3页
r语言常用统计方法实现_第4页
r语言常用统计方法实现_第5页
资源描述:

《r语言常用统计方法实现》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库

1、常用统计方法用R实现描述性统计位置的度量:均值、顺序统计量、中位数、百分位数。均值计算:若x是向量、矩阵,则mean(x)返回其全部元素均值。若要返回数组某一维的均值:apply(x,dim,mean);dim=1计算行均值,dim=2计算列均值。若x是数据框,则mean(x)返回各列的均值Mean的一般用法:mean(x,trim=0,na.rm=FALSE)trim指定去掉x两端数的比例;na.rm=TRUE允许有缺失值。类似有sum(x)函数可求x的和。顺序统计量将n个数据(观测值)按从小到大的顺序排列后,称其为顺序统计量.函数sort(x)给出了样本x的顺序统计量ord

2、er()给出排序后的下标rank()给出了样本x的秩次统计量x<-c(75,64,47.4,66.9,62.2,62.2,58.7,63.5)sort(x)order(x)中位数中位数描述数据中心位置的数字特征.大体上比中位数大或小的数据个数为整个数据的一半.对于对称分布的数据,均值与中位数比较接近;对于偏态分布的数据,均值与中位数不同.中位数的又一显著特点是不受异常值的影响,具有稳健性,因此它是数据分析中相当重要的统计量.在R软件中,函数median()给观测量的中位数.如x<-c(75,64,47.4,66.9,62.2,62.2,58.7,63.5)median(x)me

3、dian(x,na.rm=TRUE)#若数据中有缺失值百分位数百分位数(percentile)是中位数的推广.将数据按从小到大的排列后,0

4、差、变异系数和标准误等.在R软件中,用var()和sd()计算方差、标准差:var(x,na.rm=FALSE,)sd(x,na.rm=FALSE)变异系数、平方和对于变异系数、校正平方和、未校正平方和等指标,需要编写简单的程序.变异系数CV计算:cv<-100*sd(x)/mean(x);cv校正平方和CSS:css<-sum((x-mean(x))^2);css未校正平方和USS:uss<-sum(x^2);uss极差与标准误样本极差(记为R)的计算:R=max(x)-min(x)样本上、下四分位数之差称为四分位差(或半极差),记为R1.它也是度量样本分散性的重要数字特征,

5、特别对于具有异常值的数据,它作为分散性具有稳健性,因此在稳健性数据分析中具有重要作用.半极差计算:R1=quantile(x,0.75)-quantile(x,0.25)样本标准误(记为sm)定义为s/sqrt(n)样本标准误计算:sm=sd(x)/sqrt(length(x))分布形状的度量偏度系数Kurtosis是刻划数据的对称性指标.关于均值对称的数据其偏度系数为0.右侧更分散的数据偏度系数为正,左侧更分散的数据偏度系数为负.当数据的总体分布为正态分布时,峰度系数Skewness近似为0;当峰度系数为正时,两侧极端数据较多;当峰度系数为负时,两侧极端数据较少.偏度系数Sk

6、ewness样本峰度系数sk计算程序n<-length(x)m<-mean(x)s<-sd(x)sk<-n/((n-1)*(n-2))*sum((x-m)^3)/s^3计算公式峰度系数Kurtosis计算样本峰度系数ku计算程序n<-length(xm<-mean(x)s<-sd(x)ku<-((n*(n+1))/((n-1)*(n-2)*(n-3))*sum((x-m)^4)/s^4-(3*(n-1)^2)/((n-2)*(n-3)))计算公式相关分析R软件采用用cov()函数计算协方差或协方差阵,用cor()函数计算相关矩阵(相关系数)。函数cov()和cor()的使用格式

7、为:cov(x,y=NULL,use="all.obs“,method=c("pearson","kendall","spearman"))cor(x,y=NULL,use="all.obs“,method=c("pearson","kendall","spearman"))其中x是数值型向量、矩阵或数据框.y是空值(NULL,缺省值)、向量、矩阵或数据框,但需要与x的维数相一致.与cov和cor有关的函数还有:cov.wt----计算加权协方差(加权协方差矩阵);cor.test---计算

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。