Question

我正在尝试运行glmer模型，但收到一条错误消息，指出缺少起始值。我试图找到一个解决方案，但没有理解大部分提示。请问任何人都可以通过一种易于理解的方式解释什么是起始值以及解决方案是什么？

学习说明

我测量了具有不同空间量=治疗的动物之间的攻击性（计数数据）。这个想法是使用治疗作为连续变量。我测量了三次攻击性，其中在第一个年龄时所有的笔具有相同的空间（covar），然后应用处理（年龄）。在每支笔（位置）内的两个位置测量侵略。

示例数据

agg <- structure(list(age = c(21L, 21L, 39L, 39L), pen = c(1L, 1L, 1L, 1L), treatment = c(4L, 4L, 4L, 4L), position = structure(c(1L, 2L, 1L, 2L), .Label = c("AKB", "BG"), class = "factor"), jost_feed = c(NA, NA, NA, 3L), ovar_jostF_19 = c(1L, 1L, 1L, 1L), aggression = c(NA, NA, NA, 13L), covar_aggr_19 = c(12L, 15L, 12L, 15L)), .Names = c("age", "pen", "treatment", "position", "jost_feed", "covar_jostF_19", "aggression", "covar_aggr_19"), row.names = c(NA, 4L), class = "data.frame") 
#

尝试过的代码

在初步探索数据后，我决定运行高斯族并使用日志链接：

agg1<- glmer (y~ treatment + age 
            +treatment:age
            + (1|pen/covar)+(1|pen/position), 
            data= agg, 
            family = gaussian(link = "log"),
            na.action=na.omit)
summary(agg1)

glmer的起始值

0 个答案: