我正在尝试使用lme4包在R中运行下面的最大glmer模型,但到目前为止,它已经花费了超过24小时,没有任何运气。
library(lme4)
maxwvq <- glmer(match~Listgp+context+length+age+ I(age^2) + gender+ freq.+
I(freq.^2) + vowel.quality +
Listgp:context+ Listgp:length+
Listgp:age+Listgp:gender+Listgp:freq.+
Listgp:vowel.quality + context:length+context:age+
context:gender+context:freq.+
context:vowel.quality + length:age+length:gender+
length:freq.+length:vowel.quality +
age:gender+age:freq.+age:vowel.quality +
gender:freq.+ gender:vowel.quality +
freq.:vowel.quality+
(Listgp+age+gender|stimulus) +
(length+context+freq.+vowel.quality|listener),
data = SBwritten, family = "binomial",
control=glmerControl(optCtrl=list(maxfun=2e5)), nAGQ =1)下面是数据的摘录。
'> tail(SBwritten)
listener Listgp gender age vowel.quality st.vowel stimulus.presentation match length787 TQ22 TQ母接头26 i ii写入不匹配长788 TQ22 TQ母接头26 i ii写入?iis匹配长789 TQ22 TQ母接头26 u uu写入不匹配长790 TQ22 TQ母接头26 uuu写入不匹配长791 TQ22 TQ母26 u uu写入不匹配长792 TQ22 TQ母26 u uu写入不匹配长上下文自然频率。787 q实数6.40 788咽音实数0.54 789平实0.25 790强音实数0.53 791 q实数0.13 792咽实0.32‘
有没有一种方法可以加快模型的处理速度或简化模型?
感谢你的帮助。
发布于 2016-08-25 15:17:03
我解决了这个问题,并能够使模型收敛。我简化了随机效果结构,为每个随机效果只包含一个斜率(仅包含主要实验兴趣的变量);即(Listgp|刺激)+ (vowel.quality|listener)用于上述模型。
https://stackoverflow.com/questions/39082758
复制相似问题