2016-08-22 2 views
0

lme4 패키지를 사용하여 R에서 최대 glmer 모델을 실행하려고하는데 행운이 없이는 24 시간 이상이 걸렸습니다.R에서 처리 모델의 속도를 높이려면 어떻게해야합니까?

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 length 

787 TQ22 TQ 여성 26 I II 긴 788 TQ22 TQ 여성 26 I II 작성한 작성 불일치? u는 U 기입 UU 790 TQ22 TQ 여성 26 긴 기록 불일치 UU 긴 789 TQ22 TQ 여성 26 일치 IIS를 불일치가 길다 791 TQ22 TQ 여성 26 uuu 작성이 길다. 792 TQ22 TQ 여성 26 u uu 작성 시간이 길다. context nature freq. 실제 787 Q 실제 6.40 788 인두 실제 0.54 789 일반 실제 0.25 790 강조 실제 0.53 791 Q 실제 0.13 792 인두 0.32 '

모델하거나 단순화를 처리하는 속도를 높일 수있는 방법이 있나요?

감사합니다.

+0

데이터를 제공하십시오 (데이터 집합이 큰 경우 작은 예제 데이터 집합을 사용하여 재현 할 수있는 예제를 만들려고합니다). 'dput'을 사용하여 그것을 공유하거나'data()'로 내장 된 데이터 셋을 호출 할 수 있습니다. –

+0

죄송합니다. 귀하의 모델이 너무 복잡해서 무의미한 모습으로 보일 수 있습니다. 이것은 대단히 적합 할 가능성이 큽니다. – Roland

+0

감사합니다. Hack-R. 도움이된다면 데이터의 일부를 내 게시물에 추가했습니다. – Shad

답변

0

문제를 해결하여 모델을 수렴 할 수있었습니다. 각 무작위 효과에 대해 하나의 기울기 (주 실험 관심 변수 만 포함)를 포함시켜 무작위 효과 구조를 단순화했습니다. (Listgp | stimulus) + (vowel.quality | listener)입니다.

관련 문제