从glm中提取唯一系数估计值

时间:2015-10-23 15:26:22

标签: r glm lme4

我有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个人中的每个人得到它。我确定我错过了一些简单的事情......

1 个答案:

答案 0 :(得分:0)

这更像我在寻找的东西:

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