我正在尝试比较完整模型和简化模型,但是由于我上次这样做已经有一段时间了,所以我遇到了一些问题。
任务: 我有一个关于实验A和B分别或一起介绍的实验的读数(设置单人与双人)。我需要比较两个模型:完全模型,具有处理A和B之间的相互作用,以及简化模型-没有相互作用。
现在输入数据如下:
readout treatA treatB setup
85.29 1 0 single
91.74 1 0 single
94.08 1 0 single
70.76 1 0 single
97.37 1 0 single
98.23 0 1 single
61.08 0 1 single
54.18 0 1 single
53.18 0 1 single
27.93 0 1 single
55.71 1 1 double
24.83 1 1 double
43.3 1 1 double
4.52 1 1 double
36.06 1 1 double
我正在尝试使用方差分析(这是我尝试过的方法之一)进行比较:
full_model <- lm(readout ~ treatA + treatA:setup, data = data)
red_model <- lm(readout ~ treatA, data = data)
但是完整的模型给我treatA0:setup
的不适用
在安排解释变量方面,我觉得我对完整模型的输入不正确。谁能帮助我找出正确的做法?否则还有什么问题...