应用多元统计分析作业.pdfVIP

古之立大事者，不惟有超世之才，亦必有坚忍不拔之志。——苏轼

多元统计分析

实验报告

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

实验项目名称多元统计理论的计算机实现

年级2013

专业应用统计学

学生姓名侯杰

成绩

古之立大事者，不惟有超世之才，亦必有坚忍不拔之志。——苏轼

理学院

实验时间：2015年05月07日

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

姓名侯杰学号913113700124实验组

实验时间16.05.07指导教师李建军

实验项目名称多元统计理论的计算机实现

实验目的及要求：

目的：熟悉R（或SPSS）软件，掌握多元统计分析中多元正态分布均值向量和协差阵的检验，判别方法，聚类

分析，主成分分析，因子分析，相应分析内容。

要求：程序要有注释，尽量体现多元统计分析多元正态分布均值向量和协差阵的检验，判别方法，聚类分析，

主成分分析，因子分析，相应分析内容内容的基本原理。

实验硬件及软件平台：

计算机、R、网络

实验内容（包括实验具体内容、算法分析、源代码等等）：

指导教师意见：

签名：年月日

1

古之立大事者，不惟有超世之才，亦必有坚忍不拔之志。——苏轼

代码及运行结果分析

1、均值检验

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

频心电频谱值（LF）、高频心电频谱值（HF），数据见压缩包，试分析这两个指标的各次重复测定

均值向量是否有显著差异。

代码如下：

Tsq.test-function(data,alpha=0.05){

data-#读取数据

xdat-data[,2:4];

xbar-apply(xdat,2,mean);#计算LF指标的均值

ydat-data[,5:7];

ybar-apply(ydat,2,mean);#计算HF指标数据

xcov-cov(xdat);#计算LF样本协差阵

ycov-cov(ydat);#计算HF样本协差阵

sinv-solve(xcov+ycov);#求逆矩阵

Tsq-(16+16-2)*t(sqrt(16*16/(16+16)*(xbar-ybar)))%*%sinv%*%sqrt(16*16/(16+16)*(xbar-ybar));

#计算T统计量