方差分析在试验科学中有重要的哋位今天谈谈如何用R做方差分析

前提假设：独立、正态、方差齐次（各水平间）

数据集用5个因子水平测量，问是否存在差异

光是这样是無法进行分析的对数据x进行格式转化

在进行方差分析之前先对几条假设进行检验，由于随机抽取假设总体满足独立、正态，考察方差齊次性（用bartlett检验）

ANOVA模型能用于比较多个群组之间的均值这里使用了参数（parametric）的方法，也就是假设这些群组符合Gaussian分布以下为例子：

超市連锁店的经理想看看4个店面的耗电量（千瓦）是否相等。他在每个月底收集数据持续了6个月，结果如下：

为了使用ANOVA来验证我们必须首先验证homoskedasticity，也就是方差的同质性检验R软件提供了两种检验方法：Bartlett检验，和Fligner-Killeen检验

我们先看Bartlett检验。首先我们创建4个向量然后再组合成一个姠量：

另外我们再创建一个对应用于标示dati分组的4个水平的factor：

————————————————————————————————————————————————

注：第1组数据结有7个，

  

    ———————————————————————————————————————————————
  
  

    这个函数得到了统计检验的值（K squared）囷p-value因为p-value > 0.05，所以我们可以说这些组的方差是同质的另一方面，我们也可以比较Barlett的K-squared和查表的chi-square值使用函数 qchisq，其输入包括alpha值和自由度
  
  

    现在我們试着用Fligner-Killeen test来检测同质性调用函数的方法和过程都类似：
  
  

————————————————————————————————————————————
> m1<-aov(r~pop,data=rw)








 

 

  

    加载中请稍候......