主要内容

anova1

单向方差分析

描述

例子

p= anova1 (y)执行单向方差分析对于样例数据y并返回p价值。anova1将每一列的y作为一个单独的组。功能测试的假设样本的列y都来自相同的人口意味着对备择假设人口意味着并非都是相同的。函数也显示了箱线图每组的y和标准方差分析表(资源描述)。

p= anova1 (y,集团)对样本数据进行单向方差分析y分组,集团

例子

p= anova1 (y,集团,displayopt)使方差分析表和盒型图显示时displayopt“上”(默认)和抑制时显示displayopt“关闭”

例子

(p,资源描述)= anova1 (___)回报的方差分析表(包括列和行标签)单元阵列资源描述使用任何输入参数组合在前面的语法。复制文本到剪贴板版本的方差分析表,选择编辑>复制文本从方差分析表图。

例子

(p,资源描述,统计数据)= anova1 (___)返回一个结构,统计数据,你可以使用它来执行多重比较检验。多重比较测试使您能够确定哪些双组意味着是明显不同的。执行这个测试,使用multcompare,提供了统计数据结构作为输入参数。

例子

全部折叠

创建示例数据矩阵y列是常数,加上随机干扰正常的意思是0和标准偏差1。

y = meshgrid (1:5);rng默认的;%的再现性y = y + normrnd (0、1、5、5)
y =5×51.5377 0.6923 1.6501 3.7950 5.6715 2.8339 1.5664 6.0349 3.8759 3.7925 -1.2588 2.3426 3.7254 5.4897 5.7172 1.8622 5.5784 2.9369 5.4090 6.6302 1.3188 4.7694 3.7147 5.4172 5.4889

执行单向方差分析。

p = anova1 (y)

图单向方差分析包含类型的对象uicontrol。

图包含一个坐标轴对象。坐标轴对象包含35行类型的对象。一个或多个行显示的值只使用标记

p = 0.0023

方差分析表显示了组间变异()和群体内部的变化(错误)。党卫军是平方和,df的自由度。总自由度总数的观察- 1,也就是25 - 1 = 24。组织的团体之间的自由度是- 1,5 - 1 = 4。群体内部自由度的总自由度-团体之间的自由度,24 - 4 = 20。

女士是均方误差,它是什么SS / df为每一个变异的来源。的F统计量的比值均方误差(13.4309/2.2204)。的p测试统计值的概率是可以把一个值大于计算检验统计量的值,即:P (F > 6.05)。小p价值0.0023表明差异列意味着意义重大。

输入样例数据。

强度= [82 86 79 83 84 85 86 87 74 8278 75 76 77 79 79 77 78 82 79];合金= {“圣”,“圣”,“圣”,“圣”,“圣”,“圣”,“圣”,“圣”,“到”,“到”,“到”,“到”,“到”,“到”,“al2”,“al2”,“al2”,“al2”,“al2”,“al2”};

一项研究的数据结构梁的强度霍格(1987)。光束的矢量强度措施变位四舍五入英寸3000磅的力量。向量合金标识每个梁钢(“圣”1)、合金(“到”),或者合金2 (“al2”)。虽然在这个例子中,合金进行排序分组变量不需要排序。

测试零假设的钢梁相等强度的光束由两个昂贵的合金。关掉图显示单元阵列和回报方差分析结果。

(p(资源)= anova1(力量、合金“关闭”)
p = 1.5264 e-04
台=4×6单元阵列{'源'}{“党卫军”}{“df”}{‘女士’}{' F '}{'不利于> F '}{“团体”}{[184.8000]}{[2]}{[92.4000]}{[15.4000]}{[1.5264 e-04]}{‘错误’}{[102.0000]}{[17]}{[6.0000]}{0 x0双}{0 x0双}{“总”}{[286.8000]}{[19]}{0 x0双}{0 x0双}{0 x0双}

总自由度总数的观察- 1,这是 2 0 - - - - - - 1 = 1 9 。组织的团体之间的自由度是- 1,这是 3 - - - - - - 1 = 2 。群体内部自由度的总自由度-团体之间的自由度,这是 1 9 - - - - - - 2 = 1 7

女士是均方误差,它是什么SS / df为每一个变异的来源。的F统计是均方误差的比值。的p测试统计值的概率是可以把一个值大于或等于检验统计量的值。的p价值1.5264 e-04表明拒绝零假设。

您可以检索值的方差分析表索引单元阵列。保存F统计值和p新变量值函数pvalue

函数=台{2、5}
函数= 15.4000
pvalue =台{2,6}
pvalue = 1.5264 e-04

输入样例数据。

强度= [82 86 79 83 84 85 86 87 74 8278 75 76 77 79 79 77 78 82 79];合金= {“圣”,“圣”,“圣”,“圣”,“圣”,“圣”,“圣”,“圣”,“到”,“到”,“到”,“到”,“到”,“到”,“al2”,“al2”,“al2”,“al2”,“al2”,“al2”};

一项研究的数据结构梁的强度霍格(1987)。光束的矢量强度措施变位四舍五入英寸3000磅的力量。向量合金标识每个梁钢(1)、合金(到此),或者合金2 (al2)。虽然在这个例子中,合金进行排序分组变量不需要排序。

执行单向方差分析使用anova1。返回的结构统计数据,其中包含的数据multcompare需要执行多重比较

[~,~,统计]= anova1(强度、合金);

图单向方差分析包含类型的对象uicontrol。

图包含一个坐标轴对象。坐标轴对象包含21线类型的对象。一个或多个行显示的值只使用标记

p价值0.0002表明,梁的强度是不一样的。

执行一个多重比较平均强度的光束。

[c ~ ~, gnames] = multcompare(统计);

图意味着包含一个坐标轴对象的多重比较。坐标轴对象与标题点击您想要测试,包含两组明显意味着不同于圣包含7线类型的对象。一个或多个行显示的值只使用标记

在图中,蓝色的酒吧代表比较区间意味着对钢铁材料强度。红酒吧代表比较间隔指的是材料强度的合金合金1和2。红酒吧都与蓝色的酒吧,这表明,意味着对钢铁材料强度明显不同于合金合金1和2。您可以通过单击确认显著差异代表合金1和2的酒吧。

显示多重比较结果和表中相应的组名。

台= array2table (c,“VariableNames”,(“组”,“B组”,“下限”,“a - b”,“上限”,“假定值”]);资源描述。“组”)= gnames(资源描述。“组”));资源描述。“B组”)= gnames(资源描述。“B组”))
台=3×6表A组B组A - B上限下限假定值_________……___ ___________ _____{‘圣’}{'到'}3.6064 7 10.394 - 0.00016831{‘圣’}{‘al2} 1.6064 5 8.3936 - 0.0040464{'到'}{‘al2} 0.35601 1.628 -5.628 - 2

前两列显示,两组进行比较。第四列显示了估计组意味着之间的区别。第三和第五列显示的上下极限的95%置信区间的真正区别的意思。第六列显示了p价值的假设意味着相应的团体的真正差别等于零。

前两行显示两个比较涉及第一组(钢)置信区间,不包括零。因为相应的p值(分别为1.6831 e-04和0.0040)很小,这些差异是显著的。

第三行显示,两种合金的强度差异不显著。95%置信区间的差异[-5.6,1.6],所以你不能拒绝的假设真正的区别是零。相应的p第六列中的值0.3560证实了这个结果。

输入参数

全部折叠

样本数据,指定为一个向量或矩阵。

  • 如果y你是一个矢量,必须指定吗集团输入参数。中的每个元素集团代表一组对应的元素的名称y。的anova1函数将y值对应于相同的值集团同一组的一部分。组时使用这个设计有不同数量的元素(不平衡方差分析)。

    示例的样本数据输入参数y在g组输入参数g。每个元素代表一个组名称的对应元素y。

  • 如果y你是一个矩阵,不指定吗集团,然后anova1将每一列的y作为一个单独的组。在这个设计中,功能评估列的总体均值是否相等。使用这个设计时每组有相同数量的元素(平衡方差分析)。

    示例的样本数据输入参数Y矩阵形式,说明anova1如何对待每一列(Y)作为一个单独的组

  • 如果y是一个矩阵和指定吗集团,然后每个元素集团代表一个组名中相应的列y。的anova1函数将列有相同的组名是同一组的一部分。

    示例的样本数据输入参数Y矩阵形式和组输入参数组。组中的每个元素表示一组对应的列的名称。

请注意

anova1忽略任何一个y。同样,如果集团包含空的或值,anova1忽略了相应的观测y。的anova1函数执行平衡方差分析每组是否有相同数量的观察后,无视空或者函数值。否则,anova1执行不平衡的方差分析。

数据类型:|

分组变量包含组名称指定为一个数值向量,逻辑向量分类向量,字符数组,字符串数组或单元阵列特征向量。

  • 如果y每个元素是一个矢量,然后在吗集团代表一组对应的元素的名称y。的anova1函数将y值对应于相同的值集团同一组的一部分。

    示例的样本数据输入参数y在g组输入参数g。每个元素代表一个组名称的对应元素y。

    N是观察的总数。

  • 如果y是一个矩阵,然后每个元素集团代表一个组名中相应的列y。的anova1函数对待的列y这有相同的组名是同一组的一部分。

    示例的样本数据输入参数Y矩阵形式和组输入参数组。组中的每个元素表示一组对应的列的名称。

    如果你不想为矩阵示例数据指定组名y,输入一个空数组([])或省略这个论点。在这种情况下,anova1将每一列的y作为一个单独的组。

如果集团包含空的或值,anova1忽略了相应的观测y

关于分组变量的更多信息,请参阅分组变量

例子:“集团”,[1、2、1,3,1,…,3,1]y是一个向量与观测分为团体1,2,3

例子:“集团”,{‘白’,‘红’,‘白’,‘黑’,‘红色’}y是一个矩阵与5列分为组织红色,白色,和黑色的吗

数据类型:||逻辑|分类|字符|字符串|细胞

指示器显示方差分析表和箱线图,指定为“上”“关闭”。当displayopt“关闭”,anova1返回输出参数,。它不显示标准方差分析表和箱线图。

例子:p =方差分析(x,集团“了”)

输出参数

全部折叠

p价值的F以及,作为标量值返回。p值的概率F统计可以取一个值大于计算检验统计量的值。anova1测试零假设,所有组意味着相等对备择假设,至少有一组指的是不同于其他人。函数获取p从提供的价值F分布。

一个p值,小于的显著性水平,表明至少有一个样本均值明显不同于其他人。常见的重要性水平是0.05或0.01。

方差分析表,作为细胞数组返回。资源描述有六列。

定义
变化的来源。
党卫军 由于每个源的平方和。
df 与每个源相关的自由度。假设N观察和总数吗k是团体的数量。然后,N- - - - - -k是群体内部自由度(错误),k- 1是团体之间的自由度(),N- 1是总自由度。N- 1 = (N- - - - - -k)+ (k- 1)
女士 平均每个源广场,这是比例SS / df
F F统计,意思是正方形的比例。
概率F > p值,也就是概率F统计可以取一个值大于计算检验统计量的值。anova1这个概率来自的运作F分布。

方差分析表的行显示变化的数据除以源。

定义
变化造成的差异意味着(可变性之间的组)
错误 变化将在每组数据之间的差异和集团(可变性组)
总变化

统计数据多重比较测试,作为一个结构返回这个表中描述的领域。

字段名 定义
gnames 组的名称
n 观察每组的数量
的来源统计数据输出
意味着 估计价值的手段
df 错误(群体内部)自由度(N- - - - - -k,在那里N观察和总数吗k是组)的数量
年代 均方误差的平方根

更多关于

全部折叠

箱线图

anova1返回一个观察每组的箱线图y。箱形图提供了一个视觉比较组的位置参数。

在每个箱子上中央标明是中位数(2分位数,2)和盒子的边缘是第25和第75百分位数(第一和第三分位数,13分别)。胡须延伸到最极端的数据点,不被认为是离群值。离群值绘制单独使用“+”的象征。胡须对应的极端3+ 1.5×(3- - - - - -1)1- 1.5×(3- - - - - -1)

箱形图包括为中间值的比较级。两个中位数在5%显著性水平如果他们明显不同的间隔,代表等级,不重叠。这个测试是不同的F以及方差分析执行;然而,大框对应的中心线的差异很大F统计值和相应较小p价值。级距对应的极端2- 1.57 (3- - - - - -1)/√n),2+ 1.57 (3- - - - - -1)/√n),n没有任何是观测的数量值。在某些情况下,等级可以扩展外框。

箱线图的例子与标签汇总统计

箱形图的更多信息,请参阅“晶须”“缺口”箱线图

引用

[1]豪格,r . V。,和J. Ledolter.工程数据。纽约:麦克米兰,1987。

版本历史

之前介绍过的R2006a