分裂图重复测量anova，线性回归和R中的不同结果

Question

这个问题来自于对该主题的长期研究。我很感激你的回复。

我进行了一项研究，现在我的目标是测试治疗对ouctome的影响，但我想检查一下Race的适度效果。

n treatment: 83 
n placebo: 76
No missings.

我的原始数据集是宽格式的。

哪里：

**ID -> Participant identification
Con_dummy -> Treatment or Control (Between subjects)
Hispanic -> Yes or no (Between subjects)
M41 and M_42 -> Dependent variables. M_41 was measured in 2016 and M_42 was measures in 2017.**

我使用以下代码将数据集转换为长格式：

d2_stack_long <- df_stack %>%    group_by(cond_dummy) %>%   
select(ID, cond_dummy, hispanic, M_4I1TOT,M_4I2TOT) %>%   
gather(key="ID", value, M_4I1TOT,M_4I2TOT) %>%    tbl_df() %>%   
setNames(c("ID","Condition","Hispanic", "Time","Result")) %>%   
mutate(Time = ifelse(Time == "M_4I1TOT", "Time 1", "Time 2")) %>%   
arrange(ID)

我收到了一条消息

警告消息：测量变量的属性不相同;他们将被丢弃

但现在我有一个很长的数据集。我现在有两行数，因为每个参与者都是在两个时间点测量的，但一切都保持不变

> d2_stack_long %>% group_by(Condition, Time) %>% summarise(mean(Result),n=n())
# A tibble: 4 x 4
# Groups:   Condition [?]
  Condition   Time `mean(Result)`     n
      <chr>  <chr>          <dbl> <int>
1   Control Time 1       5.973684    76
2   Control Time 2       6.342105    76
3 Treatment Time 1       6.277108    83
4 Treatment Time 2       9.626506    83

现在，当我测试这个模型时，在使用R和使用JASP时，我获得了两个不同的结果。

mod <- lm(data=d2_stack_long , Result ~ Time + Condition*Hispanic)
Anova(mod, type=3)

我想我忘了调整一些东西，我会感激你的帮助。 qazxsw poi

在问这个问题之前，我读过几个帖子，但我发现有人评论或math equation to ANOVA就像comparing results across packages一样。 带有例子的here失败了。

谢谢。

Answer 1

如果有人想要答案。经过漫长的旅程，我从这个famous blog。 R中的平方和类型设置为1，而在JASP中设置为3。

您可以在回归模型中添加ERROR术语（重复测量），尽管此处有99％的帖子。

post

如果你想和aov funcion一起去：

options(contrasts=c("contr.helmert", "contr.poly"))
anova(lm(Result ~ Time*Condition*Hispanic*ID, data=d2_pals_total))

谢谢