分裂图重复测量anova，线性回归和R

Question

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

我进行了一项研究，现在我的目标是测试治疗对 ouctome 的影响，但我想检查审核效果种族即可。

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)

我想我正在努力调整一些事情，我会感激你的帮助。

在提出此问题之前，我已阅读了几篇帖子，但我发现有人对math equation to ANOVA或comparing results across packages评论here。
带有例子的famous blog已经失效。

感谢。

Answer 1

如果有人想要答案。经过一段漫长的旅程后，我从post开始。 R中的平方和类型设置为1，而在JASP中设置为3。

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

--extra_vars

如果你想和aov funcion一起去：

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

由于