条件下R中的一种方差分析和TUKEY

时间:2019-07-03 22:16:38

标签: r statistics anova tukey

我正在尝试查找变量stim_ending_t之间的平均差异,该变量包含以下6个因素:1、1.5、2、2.5、3、3.5

您可以访问df Here

stim_ending_t visbility soundvolume Opening_text               m    sd coefVar
           <dbl>     <dbl>       <dbl> <chr>                  <dbl> <dbl>   <dbl>
 1           1           0           0 Now focus on the Image  1.70 1.14    0.670
 2           1           0           0 Now focus on the Sound  1.57 0.794   0.504
 3           1           0           1 Now focus on the Image  1.55 1.09    0.701
 4           1           0           1 Now focus on the Sound  1.77 0.953   0.540
 5           1           1           0 Now focus on the Image  1.38 0.859   0.621
 6           1           1           0 Now focus on the Sound  1.59 0.706   0.444
 7           1.5         0           0 Now focus on the Image  1.86 0.718   0.387
 8           1.5         0           0 Now focus on the Sound  2.04 0.713   0.350
 9           1.5         0           1 Now focus on the Image  1.93 1.00    0.520
10           1.5         0           1 Now focus on the Sound  2.14 0.901   0.422

这是我的数据的直观表示 Here is a visual representation of my data

问:如何使用“ Opening_test”(包含“现在关注图像”和“现在关注声音”)比较均值的条件进行方差分析。

问:另外,我也想进行事后测试。

这是我尝试过的方法,但显然不是正确的方法!

# Compute one-way ANOVA test

res.aov <- aov(m ~ stim_ending_t, data = clean_test_master2)
summary(res.aov)

              Df Sum Sq Mean Sq F value Pr(>F)    
stim_ending_t  1  7.589   7.589   418.8 <2e-16 ***
Residuals     34  0.616   0.018                   
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

我认为aov结果有问题! stim_ending_t有6个因素,因此上表的Fredom(Df)的Dgree应该= 5而不是!= 1。

# post hoc test 
TukeyHSD(res.aov, conf.level = 0.99)

Here is the message I got

Error in TukeyHSD.aov(res.aov, conf.level = 0.99) : 
  no factors in the fitted model
In addition: Warning message:
In replications(paste("~", xx), data = mf) :
  non-factors ignored: stim_ending_t

注意:参与者在一个会话中以任一条件-Opening_text开始,并随机完成另一个,从而完成了实验。

1 个答案:

答案 0 :(得分:1)

  1. 以下6个因素:1、1.5、2、2.5、3、3.5

    不!如果您将其视为要素,它将是具有六个要素的一个要素。您将其用作定量变量,请参阅ANOVA表中的Df。应该改为5。1。在as.factor()之前尝试aov函数。

  2. m是因变量吗?如果是,visbilitysoundvolume是什么?如果它们也是因素,则独立性的假设是错误的。在这种情况下,您应该引入这些因素进行建模。