Извлечение оценок уникальных коэффициентов из glm

У меня есть 57 особей (с ~ 60 000 повторений), помеченных «ID_Year», которые используются для создания модели выбора среды обитания GLME. Как извлечь уникальные оценки коэффициентов для каждого "ID_Year"?

Вот модель:

Cand.mod[[25]] <-glmer(Used ~ WBEMesicHa + Treatments + WBERoads45 + WBERoads25 + WBERoadsUn + Powerlines + Agricultur + UrbanEucli+ NonUrbanEu + ClippedCan + (1|ID_Year),family=binomial(logit),data=data)#failed to converge

Я попробовал fixef(Cand.mod[[25]]) и получил один вывод оценок модели:

 (Intercept)   WBEMesicHa   Treatments   WBERoads45   WBERoads25   WBERoadsUn   Powerlines   Agricultur   UrbanEucli   NonUrbanEu  ClippedCan2  ClippedCan3 
-4.195335406  0.045302307 -0.111918732 -0.036906751  0.025336140  0.240330206  0.209806088  0.073736439 -0.004829703 -0.217925175 -1.070915726 -1.866897616 
 ClippedCan4  ClippedCan5 
-2.833117637 -5.124756532

Но я просто не могу получить его для каждого из 57 человек. Я уверен, что упускаю что-то простое...


r lme4 glm
person Charlie    schedule 23.10.2015    source источник
comment
@Roland Я пробовал это раньше, и хотя он возвращает оценки для всех 57 человек, единственный коэффициент, который меняется в зависимости от человека, - это перехват. Остальные коэффициенты остаются неизменными...   -  person Charlie    schedule 23.10.2015
comment
Да, конечно. В конце концов, у вас есть только случайный перехват. Наклоны должны быть одинаковыми для всех особей. Если вам нужны случайные наклоны, вам нужно указать это. (Я заметил комментарий о сбое сходимости. Вы должны изучить это.)   -  person Roland    schedule 23.10.2015
comment
@ Роланд, я хотел удалить этот комментарий в конце. Я исследовал это с предложениями Бена Болкера - пока это не проблема. Как указать случайные наклоны?   -  person Charlie    schedule 23.10.2015
comment
stats.stackexchange.com/a/13173/11849   -  person Roland    schedule 23.10.2015

Ответы (1)


arrow_upward
0
arrow_downward

Это больше похоже на то, что я искал:

http://www.inside-r.org/packages/cran/lme4/docs/isREML

person Charlie    schedule 28.10.2015