统计软件实验教学要件 实验项目十 R软件中数据的描述.docVIP

实验项目十 R软件中数据的描述 一、实验目的 运用R中的barplot( )函数画条形图，hist( ) 函数画直方图，pie( )函数画饼图，boxplot( )函数画箱线图，plot( ) 函数画线图，stem( )函数画茎叶图。 掌握运用R软件中的函数对数据进行描述性分析。 二、实验环境 系统软件Windows2000或WindowsXP或Windows7； 统计软件R2.13.2或更高版本。 三、实验内容 1.常见统计图的绘制 (1) 分类数据和顺序数据的频数分布表 (2) 分类数据和顺序数据的条形图及饼图 (3) 数值型数据的频数分布表 (4) 数值型数据的直方图、箱线图、线图、茎叶图 2. 描述性统计分析 第10章 数据的描述 10.1 常见统计图的绘制 10.1.1 分类数据和顺序数据的频数分布表 10.1.2 分类数据和顺序数据的条形图及饼图 10.1.3 数值型数据的频数分布表 10.1.4 数值型数据的直方图、箱线图、线图、茎叶图 10.2 描述性统计分析 10.1 常见统计图的绘制 10.1.1 分类数据和顺序数据的频数分布表 【例10.1】(数据文件为li10.1.txt)对某高校经济系30名教师性别及职称登记结果，如表10.1所示，试用R分别编制教师职称的频数分布表 。 【统计理论】 对分类数据和顺序数据，计算出每一类别出现的频数或频率后，可通过频数分布表 来展示。 【软件操作】 采用table( )函数编制频数分布表 setwd(“D:/R-Statistics/data/chap-10”) #设定工作路径 dat=read.table(li10.1.txt,header=T) #从li10.1.txt中读入数据，记为dat fre-table(x$ZHCH) # 生成教师职称变量的分组频数，如果绘制教师性别的频数分布表只需要把命令改为fre- table(x$ZHCH)即可 library(sca) # 加载扩展包sca per-percent(fre/sum(fre),d=1) # 使用扩展包sca中的percent函数，可以用带%符号的表示方式显示结果，参数d来指定小数点后的表示位数。在不指定的情况下，小数点后的数字不表示。 rbind(fre,per) # 合并表格命令，生成教师职称分组频数分布表 运行结果 副教授 讲师 教授 助教 fre 12 8 6 4 per 40.0 % 26.7 % 20.0 % 13.3 % 10.1.2 分类数据和顺序数据的条形图及饼图 【例10.2】 (数据文件为li10.1.txt) 根据表10.1资料, 用R软件绘制条形图。 （数据来源：费宇等,《统计学》第2章,高等教育出版社,2010） 【统计理论】 对分类数据和顺序数据，计算出每一类别出现的频数或频率后，还可通过统计图来展示，例如条形图 。 【软件操作】 采用barplot()函数画条形图 setwd(“D:/R-Statistics/data/chap-10”) #设定工作路径 dat=read.table(li10.1.txt,header=T) #从li10.1.txt中读入数据，记为dat ZHCH.table-table(x$ZHCH) # 生成教师职称变量的列联表 barplot(sort(ZHCH.table,decreasing=TRUE),ylab=人数) # 其中sort命令是对列联表进行排序，decreasing=TRUE按升序排，参数ylab为y轴的标签 10.1.3 数值型数据的频数分布表 【例10.4】 (数据文件为li10.4.txt)2007年我国各地区农村居民家庭人均纯收入资料如表10.2所示，试用R软件制作频数分布表。 【统计理论】 数值型数据经过整理后，首先可用频数分布表(变量数列) 将其展示出来。 【软件操作】 x=read.table(li10.4.txt,header=T) #从li10.4.txt中读入样本数据x fre- table(x[,4]) #生成教师职称变量的分组频数 library(sca) # 加载扩展包sca per-percent(fre/sum(fre),d=1) # 使用扩展包sca中的percent函数，可以用带%符号的表示方式显示结果，参数d来指定小数点后的表示位数。在不指定的情况下，小数点后的数字不表示。 rbind(fre,