R 如何在线性混合模型中求F和p值

我仍然不太熟悉R，我想使用

lmer

函数构建混合效果的线性模型。
我们有3个自变量，每次有2种模式：

• 组（G3和G5）
• 频率（低频和高频）
• 条件（NR和PH）

和一个因变量：

• 反应时间（RT）

是我的数据样本，数据提取

manipfreq

：

subject group   item    RT  frequency   condition
502 G3  ARGENT  677 HF  NR
502 G3  ACTION  814 HF  NR
502 G3  AFFAIRE 1278    HF  NR
502 G3  ARMOIRE 716 LF  PH
502 G3  AIGLE   761 LF  PH
503 G3  AFFAIRE 841 HF  PH
503 G3  ARGENT  919 HF  PH
503 G3  ACTION  1071    HF  PH
503 G3  AIGLE   1253    LF  NR
505 G3  ARGENT  785 HF  PH
505 G3  AFFAIRE 1028    HF  PH
505 G3  ACTION  1190    HF  PH
505 G3  AIGLE   1248    LF  NR
506 G3  ACTION  898 HF  NR
506 G3  AFFAIRE 1566    HF  NR
506 G3  ARGENT  2037    HF  NR
506 G3  AIGLE   707 LF  PH
506 G3  ARMOIRE 839 LF  PH
1530    G5  ACTION  705 HF  NR
1530    G5  ARGENT  755 HF  NR
1530    G5  AFFAIRE 789 HF  NR
1530    G5  AIGLE   481 LF  PH
1530    G5  ARMOIRE 495 LF  PH
1532    G5  ACTION  532 HF  NR
1532    G5  ARGENT  679 HF  NR
1532    G5  AFFAIRE 1142    HF  NR
1532    G5  ARMOIRE 423 LF  PH
1532    G5  AIGLE   484 LF  PH
1533    G5  ACTION  690 HF  PH
1533    G5  ARGENT  708 HF  PH
1533    G5  AFFAIRE 775 HF  PH
1533    G5  AIGLE   1027    LF  NR
1533    G5  ARMOIRE 1299    LF  NR
1539    G5  ARGENT  747 HF  NR
1539    G5  ACTION  771 HF  NR
1539    G5  AFFAIRE 827 HF  NR
1539    G5  ARMOIRE 518 LF  PH
1539    G5  AIGLE   566 LF  PH

我特别感兴趣的是：

• “组”、“频率”、“条件”的主要影响
• 双向交互：“条件*组”和“条件*频率”
• 三方互动：“组*频率*条件”

以及仅在NR条件下的双向交互“组*频率”

我写了一个模型，该模型考虑了我的IV与截距和坡度上与我的主题和项目相关的可变性之间的相互作用。
这是我的模型（如果我错了，请纠正我）：

我需要找到F和p的值

---

我将非常感谢您的帮助。

请查看

lmerTest

包。请阅读

help（“pvalues”）

。您不应在模型中包含

（1 |主题）+（1 |项）

；它们与

（1+条件|主题）+（1+条件|项目）

model=lmer(RT~frequency*condition*group+(1|subject)+(1|item)+(1+condition|subject)+(1+condition|item),data=manipfreq)