如何将协变量作为glmer的固定因素?

时间:2014-10-14 10:57:51

标签: r

我是R的初学者,所以如果我问一个愚蠢的问题,请原谅我。

我在语法结构(OS.SO.sentence;主语优先和对象优先结构)上进行了一个关于语速减慢(speech.rate;三种语速)效果的实验,看看语音是否会减慢便于理解,特别是对象优先结构(精度测量,二进制代码)。 我测试了听力受损的不同年龄组的非听力障碍者。众所周知,与老年人一样,听力受损的成年人也会从较慢的语音中获益,因此在OS结构中受损更多。

在SPSS中(以及如何在SPSS仍然花哨的时候错过这些时间......)我本来能够将hear.state(二进制,2个值)和age.group(二进制,4个值)作为协变量放入推广估计方程。或者,我可以将听力和年龄视为连续变量...在R中,我无法弄清楚如何将这两个变量拟合为协变量。我不认为我能用与语速和句法结构相同的力量来对待它们。

到目前为止,我已经尝试了几种模型,将所有模型定义为固定因素:

gm1 = glmer(Accuracy ~ OS.SOSentenceType + (1|Subject) + (1|ItemID), data=dataC, family='binomial', REML = FALSE)

gm1a = glmer(Accuracy ~ OS.SOSentenceType + speech.rate + (1|Subject) + (1|ItemID), data=dataC, family='binomial', REML = FALSE)

gm1b = glmer(Accuracy ~ OS.SOSentenceType  +  speech.rate + hearing.state  + (1|Subject) + (1|ItemID), data=dataC, family='binomial', REML = FALSE)

gm1c = glmer(Accuracy ~ OS.SOSentenceType  +  speech.rate + hearing.state  + age.group + (1|Subject) + (1|ItemID), data=dataC, family='binomial', REML = FALSE)

gm2 = glmer(Accuracy ~ OS.SOSentenceType*speech.rate + hearing.state  + age.group + (1|Subject) + (1|ItemID), data=dataC, family='binomial', REML = FALSE)

gm2a = glmer(Accuracy ~ OS.SOSentenceType*speech.rate*hearing.state  + age.group + (1|Subject) + (1|ItemID), data=dataC, family='binomial', REML = FALSE)

gm2b = glmer(Accuracy ~ OS.SOSentenceType*speech.rate *hearing.state*age.group + (1|Subject) + (1|ItemID), data=dataC, family='binomial', REML = FALSE)

模型越大,实际解释的越少...... Atm,我觉得在黑暗中钓鱼甚至没有钓鱼竿。 非常感谢帮助!

干杯,格拉

0 个答案:

没有答案