R语言与统计分析大作业

第五章

5.1 设总体x是用无线电测距仪测量距离的误差，它服从（α，β）上的均匀分布，在200次测量中，误差为xi的次数有ni次： Xi：3 5 7 9 11 13 15 17 19 21 Ni：21 16 15 26 22 14 21 22 18 25 求α，β的矩法估计值

α=u- 3s β=u+ 3s 程序代码： x=seq(3,21,by=2)

y=c(21,16,15,26,22,14,21,22,18,25) u=rep(x,y) u1=mean(u) s=var(u) s1=sqrt(s) a=u1-sqrt(3)*s1

b=u1+sqrt(3)*s1b=u1+sqrt(3)*s1 得出结果： a= 2.217379 b= 22.40262

5.2为检验某自来水消毒设备的效果，现从消毒后的水中随机抽取50L，化验每升水中大肠杆菌的个数（假设1L水中大肠杆菌的个数服

从泊松分布），其化验结果如下表所示：试问平均每升水中大肠杆菌个数为多少时，才能使上述情况的概率达到最大 大肠杆菌数/L：0 1 2 3 4 5 6 水的升数：17 20 10 2 1 0 0 γ=u是最大似然估计 程序代码： a=seq(0,6,by=1) b=c(17

中国虐学c提 2 1,72)1714丕 0 1 (7:2—3 2C ie eAg iutr l ce c l t h n s rc l a in eBu l i u S en

基于 R语言的 4 4’‘ 1’ 3肥效试验的统计分析张祯勇，高明文，肖启银，杨开俊，罗孝贵，陈代全(四川省甘孜藏族自治州农业科学研究所，四川康定 6 6 0 ) 2 0 0

摘要：了介绍R语言在农业试验分析中的应用，用R语言对“44肥料效应试验具体实例的部分为使 3 1”实施方案和完全实施方案结果进行统计分析，拟合了二元二次、三元二次肥料效应方程，并根据肥效方

程模拟寻优。使用R语言对“4 4肥料效应试验分析结果表明： 31”通过三元二次肥料效应函数获得的玉米最佳施肥指标为：以每公顷施 N 7 . gPO 8.2 gK2 53 g期望产量为 1741 gh, 252 k, 66, O4 .2, 6 k k 06 .3/m2产 k投比达到9 3为最佳方案。这与使用全国农技中心提供的“4 4 ., 7 3 1田间试验设计与数据分析管理系统”进行计算是一致的；考虑到田间的可操作性，如建议每公顷施 N 2 2 8 27 9k,: 2 2 9 . , 6. - 8 . g P 8

第一讲作业只要在计算机完成即可，不要写在作业本。

第二讲 作业

1、欲研究不同收入的群体对某种特定的商品是否有相同的购买习惯，市场研究人员调查了四个收入的组的消费者共527人，购买习惯分为：经常购买、不购买、有时购买。调查结果如下： 购买习惯 低收入组 偏低收入组 偏高收入组 高收入组 经常购买 25 40 47 46 不购买 69 51 74 57 有时购买 36 26 19 37 试分别绘制简单条图、分组条图和堆积图来分析不同收入组的购买习惯。 2、对第一讲作业中学生消费调查的数据文件，利用描述统计过程分析这31名学生消费的分布情况。

（1）计算平均月生活费、月平均衣物消费的方差、标准差、全距、最大值、最小值、中位数、众数、偏度与峰度等指标； （2）列出频数分布表；

（3）画出带有正态曲线的频数直方图；

（4）将平均月生活费范围从400以下、400到600、600到800、800以上进行分组，再做频数分布。

3、案例分析：亚太地区25所知名商学院有关情况 一、案例背景

寻求工商管理专业较高的学历已是一种世界趋势。有调查表明，越来越多的亚洲人选择就读工商管理硕士（MBA）学位，将它当作通往成功的道路。在近一年来，亚太地区的学院里申请

统计分析与SPSS的应用

原始数据：10级市场营销2班成绩

分析一：综测成绩四分位数

Statistics 综测成绩 N Valid Missing Minimum Maximum Percentiles 25 50 75

上表表明：综测成绩的最小值为68.61分，最大值为89.15分。其中25%的学生综测成绩为74.4100分，50%的学生综测成绩为80.3740分，75%的学生综测成绩为85.2200分。四分位数差从侧面证实了学生综测成绩呈一定左偏分布。

31 0 68.61 89.15 74.4100 80.3740 85.2200 分析二：综测成绩直方图

上图表明：该班学生的综测成绩均分为80.07分，标准差为5.62。从图中可以看出，综测成绩呈左偏性分布，在85分左右的学生人数最多，70分左右的学生人数最少。

分析三：综测成绩的基本统计量分析 Descriptive Statistics N Range Minimum Maximum Mean Std. Deviation Skewness Std. Kurtosis Std. Statistic Statistic Statistic Statistic Stat

统计分析与SPSS的应用

摘要：为对统计分析与spss应用分析所学知识进行巩固和检验，特运用所学知识进行简单的

统计分析应用，下文以某校学生学期成绩进行模拟分析。

一：原始数据：10级市场营销2班成绩

分析一：综测成绩四分位数

Statistics 综测成绩 N Valid Missing Minimum Maximum Percentiles 25 50 75

上表表明：综测成绩的最小值为68.61分，最大值为89.15分。其中25%的学生综测成绩为74.4100分，50%的学生综测成绩为80.3740分，75%的学生综测成绩为85.2200分。四分位数差从侧面证实了学生综测成绩呈一定左偏分布。

31 0 68.61 89.15 74.4100 80.3740 85.2200 分析二：综测成绩直方图

上图表明：该班学生的综测成绩均分为80.07分，标准差为5.62。从图中可以看出，综测成绩呈左偏性分布，在85分左右的学生人数最多，70分左右的学生人数最少。

分析三：综测成绩的基本统计量分析 Descriptive Statistics N Range Minimum Maximum Mean Std. Deviation Sk

多元统计分析

实验报告

实验课程名称多元统计分析

实验项目名称多元统计理论的计算机实现 年级 2013 专业应用统计学 学生姓名侯杰 成绩

理学院

实验时间：2015 年05月07日

学生所在学院：理学院专业：应用统计学班级：9131137001

姓名 实验时间 实验项目名称 侯杰 16.05.07 学号 指导教师 913113700124 李建军 实验组 多元统计理论的计算机实现 实验目的及要求： 目的：熟悉R（或SPSS）软件，掌握多元统计分析中多元正态分布均值向量和协差阵的检验，判别方法，聚类分析，主成分分析，因子分析，相应分析内容。 要求：程序要有注释，尽量体现多元统计分析多元正态分布均值向量和协差阵的检验，判别方法，聚类分析，主成分分析，因子分析，相应分析内容内容的基本原理。 实验硬件及软件平台： 计算机、R、网络 实验内容（包括实验具体内容、算法分析、源代码等等）： 指导教师意见： 签名：年月日

1

代码及运行结果分析

1、均值检验

问题重述：某医生观察了16名正常人的24小时动态心电图，分析出早晨3小时各小时的低

频心电频谱值（LF）、高频心电频谱值（HF），数据见压缩

多元统计分析

实验报告

实验课程名称多元统计分析

实验项目名称多元统计理论的计算机实现 年级 2013 专业应用统计学 学生姓名侯杰 成绩

理学院

实验时间：2015 年05月07日

学生所在学院：理学院专业：应用统计学班级：9131137001

姓名 实验时间 实验项目名称 侯杰 16.05.07 学号 指导教师 913113700124 李建军 实验组 多元统计理论的计算机实现 实验目的及要求： 目的：熟悉R（或SPSS）软件，掌握多元统计分析中多元正态分布均值向量和协差阵的检验，判别方法，聚类分析，主成分分析，因子分析，相应分析内容。 要求：程序要有注释，尽量体现多元统计分析多元正态分布均值向量和协差阵的检验，判别方法，聚类分析，主成分分析，因子分析，相应分析内容内容的基本原理。 实验硬件及软件平台： 计算机、R、网络 实验内容（包括实验具体内容、算法分析、源代码等等）： 指导教师意见： 签名：年月日

1

代码及运行结果分析

1、均值检验

问题重述：某医生观察了16名正常人的24小时动态心电图，分析出早晨3小时各小时的低

频心电频谱值（LF）、高频心电频谱值（HF），数据见压缩

多元统计分析

实验报告

实验课程名称多元统计分析

实验项目名称多元统计理论的计算机实现 年级 2013 专业应用统计学 学生姓名侯杰 成绩

理学院

实验时间：2015 年05月07日

学生所在学院：理学院专业：应用统计学班级：9131137001

姓名 实验时间 实验项目名称 侯杰 16.05.07 学号 指导教师 913113700124 李建军 实验组 多元统计理论的计算机实现 实验目的及要求： 目的：熟悉R（或SPSS）软件，掌握多元统计分析中多元正态分布均值向量和协差阵的检验，判别方法，聚类分析，主成分分析，因子分析，相应分析内容。 要求：程序要有注释，尽量体现多元统计分析多元正态分布均值向量和协差阵的检验，判别方法，聚类分析，主成分分析，因子分析，相应分析内容内容的基本原理。 实验硬件及软件平台： 计算机、R、网络 实验内容（包括实验具体内容、算法分析、源代码等等）： 指导教师意见： 签名：年月日

1

代码及运行结果分析

1、均值检验

问题重述：某医生观察了16名正常人的24小时动态心电图，分析出早晨3小时各小时的低

频心电频谱值（LF）、高频心电频谱值（HF），数据见压缩

多元统计分析

实验报告

实验课程名称多元统计分析

实验项目名称多元统计理论的计算机实现 年级 2013 专业应用统计学 学生姓名侯杰 成绩

理学院

实验时间：2015 年05月07日

学生所在学院：理学院专业：应用统计学班级：9131137001

姓名 实验时间 实验项目名称 侯杰 16.05.07 学号 指导教师 913113700124 李建军 实验组 多元统计理论的计算机实现 实验目的及要求： 目的：熟悉R（或SPSS）软件，掌握多元统计分析中多元正态分布均值向量和协差阵的检验，判别方法，聚类分析，主成分分析，因子分析，相应分析内容。 要求：程序要有注释，尽量体现多元统计分析多元正态分布均值向量和协差阵的检验，判别方法，聚类分析，主成分分析，因子分析，相应分析内容内容的基本原理。 实验硬件及软件平台： 计算机、R、网络 实验内容（包括实验具体内容、算法分析、源代码等等）： 指导教师意见： 签名：年月日

1

代码及运行结果分析

1、均值检验

问题重述：某医生观察了16名正常人的24小时动态心电图，分析出早晨3小时各小时的低

频心电频谱值（LF）、高频心电频谱值（HF），数据见压缩

多元统计分析

现实中的统计对象经常用多个指标来表示，比如人口普查，就可以有姓名、性别、出生年月日、籍贯、婚姻状况、民族、政治面貌、地区等，企业调查，可以有净资产、负债、盈利、职工人数、还贷情况等等。多个指标（变量）可以分别进行分析，但是，我们往往希望综合使用这些指标，这时，有主成分分析、因子分析等方法可以把数据的维数降低，同时又尽量不损失数据中的信息。 一、主成分分析

1、主成分分析的原理

多变量的主成分分析是在不损失或很少损失原有信息(指方差)的前提下，将原来多个彼此相关的指标转换为新的少数几个彼此独立的综合指标的一种统计分析方法。

多变量的主成分分析在教育评估中可用以寻找反映或影响评估对象的综合指标。如描述教师能力的指标很多：对教育对象的控制能力，对教育影响的控制能力，表达能力、教学思维能力、创新能力、组织协调能力等等。这些评估教师能力的指标个数可能很多，且指标之间彼此相关，多变量的主成分分析就是要综合这些指标，从而找出反映教师能力的少数几个彼此独立的指标，以便综合出教师能力的重要信息。

主成分分析的目的是从原始的多个变量取若干线性组合，能尽可能多地保留原始变量中

?的信息。从原始变量到新变量是一个正交变换（坐标变换）。设有