重复测量方差分析

重复测量方差分析

<script type="text/x-mathjax-config">MathJax.Hub.Config({tex2jax:{inlineMath:[['','']]}});</script>
<script type="text/javascript" src="http://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML"></script>

需要满足的前提条件

独立。各组数据独立,互不相关。这一点一般都可以满足,所以一般不做检验。

  1. 正态性。各组数据服从正态分布。
  2. 方差齐性。各组方差相等。
  3. 球对称假设。对于自变量的各取值水平组合而言(对于被试内因素的各个水平组合而言),因变量的协方差矩阵相等。
    前面1,2是方差分析都要满足的基本条件,第3条球对称假设是重复测量方差分析需要特别满足的一个额外条件。
    球形检验(Mauchly’s test of sphericity),适用于重复测量时,检验不同测量之间的差值的方差是否相等,用于三次以及三次以上水平的重复测量(因为两次重复测量根本就没有办法比较差值的方差,只有一个方差)。
    image

以上三种条件spss操作:

  • 检查正态性。
    参数检验要求数据服从正态分布。非参数检验虽然不要求总体正态分布,但能效性不高。
    方法有好以下几种:
    1. 峰度(kurtosis)和偏度(skewness)。n<100
      峰度:反映了频数分布曲线顶端的尖峭或扁平程度。当样本数据的分布符合正态分布时,峰度系数K为0,不高不低;系数K大于0时,为高狭峰(频数过于集中);系数K小于0时,为低阔峰(频数过于分散)。系数K值越大,偏离正态的情况越明显。
      偏度:描述了分布偏离对称性的程度。当样本数据的分布符合正态分布时,偏度系数S为0,不左不右;系数S大于0时,为正偏态(极端值偏向右侧高值端);系数S小于0时,为负偏态(极端值偏向左侧低值端)。系数S值越大,偏离正态的情况越明显。
      spss操作:
      分析-描述统计-描述-选变量-然后点击选项-峰度和偏度勾选
      在使用峰度和偏度来判断数据是否满足正态分布时,不一定要K和S完全等于0,很少有数据可以满足。只要K和S的值与0的差异不满足统计学意义,就仍然可以认为是正态分布。计算公式如下:
      Z = (峰度值K-零)/峰度的标准误
      Z = (偏度值S-零)/偏度的标准误
      Z_{峰度}= \frac{K-0}{σ}
      Z_{偏度} = \frac{K-0}{σ}
      将得到的Z值,也就是所谓的峰度和偏度系数与1.96进行比较,如果峰度(偏度)系数的绝对值大于或者等于1.96,就认为峰度(偏度)显著不为0,意味着变量不服从正态分布。反之,必须峰度和偏度系数的绝对值都小于1.96,则认为满足正态分布
      因为标准误会受到样本量的影响,很敏感,所以一般在样本量小于100的时候比较常用峰度和偏度来检验正态分布。
      样本量大于100,建议使用k-s正态性检验。
      image
      image
      image
    2. k-s检验和s-w检验。n>100


      image
      image
      image

      通过查看p值,如果p>0.05,则认为满足正态分布。这种方法也可以输出峰度和偏度值。

    3. p-p图。n>100
      分析-描述统计-p-p图-选变量-然后检验分布选择正态-确定
      散点如果全部落在了45度角的对角线上,就可以认为满足正态分布,只要不偏离对角线太多就可以近似认为是正态分布。
    4. q-q图
    5. 直方图
      绘图法往往没有直接检验p值准确。所以通常在查看p值后,发现是非正态时,在画图来看下数据偏态是不是很严重。如果不严重,还可以近似正态分布来做参数检验。
  • 方差齐性。
    重复测量-选项-同质性检验(Homogeneity test)
  • 球形检验。
    通过spss做重复测量的时候自动输出。


    image

    如上图,因为此数据是每个因素只有两个水平,不需要进行球形检验,所以没有输出p值。

重复测量方差分析spss操作

  • 分析-一般线性模型-重复测量。
    image
  • 填入被试内因子-添加-定义。
    image
  • 主体内变量-因子列表。
    image
  • 模型。
    image

    全因子模型是指所有的主效应和交互作用都包含在内的模型。一般默认这个不用改动。
    定制模型是自己可以定制只做某一因素的主效应,或者只做交互作用等。

  • 绘图。
  • 事后多重比较。
    也叫事后检验/多重比较。主要用于,当自变量有3个以上的水平时,并且自变量主效应显著。为了区分究竟是自变量的哪两个水平之间差异显著,做事后比较
    可以通过点击 事后多重比较来操作。
    image
    image

    也可以通过点击 选项-估计边际平均值来操作。
    image
    image

    两者得到的结果是一致的。只不过事后多重比较只能对组间因子进行操作,而选项那里可以选择所有的,组间,组内或者交互都可以选择。
  • 选项。
    可以选择 描述统计,功效估计,同质性检验等。按需选择。
  • 简单效应。
    交互作用显著时,需要做简单效应分析成对比较,两两比较)。
    简单效应分析在spss中没有专门的按钮,需要通过写代码来实现,代码如下:
    /EMMEANS=TABLES(a*b) COMPARE(a) ADJ(SIDAK)
    在b的某一个水平上,检验a变量不同水平差异的显著性。如,在b1上看,a1与a2的差异是否显著。
    /EMMEANS=TABLES(a*b) COMPARE(b) ADJ(SIDAK)
    点击 粘贴-插入代码,具体可以自己修改变量名称-然后点击运行-全部。
    image
    image
    image

输出结果

需要看球形检验方差齐性检验

image
。如果满足方差齐性,则在看下面的。如果不满足,没办法做方差分析。如果满足了方差齐性,球形检验不显著,则在看主体内效应的检验输出结果时,需要看第一行假设为球形。如果满足方差齐性,不满足球形假设,建议看第二行Greenhouse-Geisser的结果。
image

主体内效应的检验主体间效应的检验
image
。来判断主效应和交互作用是否显著。
如果主效应显著,看事后检验的结果。
如果交互作用显著,看简单效应的结果。

事后检验和简单效应检验的结果怎么报告

image

上图中的红圈所示部分,就是简单效应的报告。里面的F值和P值就是简单效应结果里面的单变量检验。单变量检验相当于把每个因素的每个水平分别做方差分析,看在一个自变量在另一个自变量的不同水平上是不是对因变量有不同的影响。上图的没有报告完整,还应该加上,到底是那个水平显著高于哪个(p值)。这个时候需要看成对比较的结果。

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 199,902评论 5 468
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 84,037评论 2 377
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 146,978评论 0 332
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 53,867评论 1 272
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 62,763评论 5 360
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,104评论 1 277
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,565评论 3 390
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,236评论 0 254
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,379评论 1 294
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,313评论 2 317
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,363评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,034评论 3 315
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,637评论 3 303
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,719评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 30,952评论 1 255
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 42,371评论 2 346
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 41,948评论 2 341

推荐阅读更多精彩内容

  • Chapter1 什么是统计学(statistics)?统计学是描述一系列可用于描述/整理/解释资料或数据的统计工...
    芒果芭乐阅读 5,016评论 0 18
  • 1. 简述相关分析和回归分析的区别和联系。 回归分析和相关分析都是研究两个或两个以上变量之间关系的方法。 广义上说...
    安也也阅读 8,634评论 0 3
  • 数据分析方法分为四大类: 1、单纯的数据加工方法 a.描述性统计分析(集中、离中趋势分析和数据分布) b.相关性分...
    重生之魂阅读 2,792评论 0 3
  • 当才到中午的时候,你就觉得已经精疲力竭了。你挣扎着保持清醒可能没有用,你不知道怎样才能阻止自己在错误的时间里犯困。...
    杂技团暖男阅读 3,194评论 0 3