У меня 57 человек (с ~ 60 000 повторов) с надписью «ID_Year», используемой для создания модели выбора среды обитания GLME. Как извлечь уникальные оценки коэффициентов для каждого «ID_Year»?Извлечь уникальные оценки коэффициентов от glm
Вот модель:
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 физические лица. Я уверен, что я что-то просто не хватает ...
Может использовать 'coef'? – Roland
@Roland Я уже пробовал это раньше, и хотя он возвращает оценки для всех 57 человек, единственным коэффициентом, который изменяется индивидуальным, является перехват. Остальные коэффициенты остаются статическими ... – Charlie
Да, конечно. В конце концов, у вас есть только случайный перехват. Склоны должны быть одинаковыми для всех людей. Если вы хотите случайные склоны, вам нужно указать это. (Я заметил комментарий об ошибке конвергенции. Вы должны изучить это.) – Roland