从 glm 中提取唯一系数估计
Extract unique coefficient estimates from 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 个人中的每一个人都得到它。我确定我遗漏了一些简单的东西...
我有 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 个人中的每一个人都得到它。我确定我遗漏了一些简单的东西...