錯誤在glmer：在EVAL錯誤（表達式，ENVIR，enclos）：無法找到有效的開始>值：請指定一些

我試圖運行glmer模型：錯誤在glmer：在EVAL錯誤（表達式，ENVIR，enclos）：無法找到有效的開始>值：請指定一些

這是我裝盤運行模型

md$count<-as.integer(md$count) 
Model1 <- glmer (count ~ distance_nest_tag + ep_nest + (1|tag_ID), 
       family = gaussian (link="log"),   
       data=md)

count是從0到19的整數，其指示的訪問次數。
distance_nest_tag是以米爲單位的連續可變
ep_nest是二進制變量（1,0）

和隨機效果方面都只是巢和個人號碼。

我也試過這種

Model1 <- glmer (count ~ distance_nest_tag + ep_nest + (1|tag_ID), 
       family = gaussian (link="log"), 
       start = coef (lm(md$count ~ md$distance_nest_tag + md$ep_nest)),     
       data=md)

但我總是得到這個

> Error in eval(expr, envir, enclos) : cannot find valid starting 
> values: please specify some

有什麼想法？？？

dput 



     Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod'] 
Family: binomial (logit) 
Formula: EP0 ~ counts + dist_mean_hour + nest_reproductive_period + X.vistors + 
    (1 | nest_ID) + (1 | tag_ring_ID) + (1 | distance_nesttag_nest) +  (1 | date) 
    Data: only_visits_encounternet 
    AIC  BIC logLik deviance df.resid 
34.7643 81.1238 -7.3822 14.7643  752 
Random effects: 
Groups    Name  Std.Dev. 
date     (Intercept) 8.507 
distance_nesttag_nest (Intercept) 68.095 
tag_ring_ID   (Intercept) 64.491 
nest_ID    (Intercept) 2.251 
Number of obs: 762, groups: date, 96; distance_nesttag_nest, 35; tag_ring_ID, 20; nest_ID, 19 
Fixed Effects: 
       (Intercept)      counts    dist_mean_hour 
       -13.08952     -1.80786     -0.04528 
nest_reproductive_periodI nest_reproductive_periodP     X.vistors 
       -4.04944     -1.56264     0.06601 
convergence code 0; 1 optimizer warnings; 0 lme4 warnings 
> dput(R) 
structure(list(nest_ID = c(6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L), nest_reproductive_period = structure(c(3L, 
1L, 2L, 3L, 1L, 2L, 3L, 1L), .Label = c("F", "I", "P"), class = "factor"), 
    tag_ID = c(23L, 23L, 23L, 4L, 4L, 4L, 5L, 5L), ring_tag_id = structure(c(3L, 
    3L, 3L, 1L, 1L, 1L, 2L, 2L), .Label = c("AU34180", "AU70442", 
    "BD12273"), class = "factor"), age_tag = c(5L, 5L, 5L, 4L, 
    4L, 4L, NA, NA), p3_tag = c(590L, 590L, 590L, 585L, 585L, 
    585L, 590L, 590L), weight_tag = c(1770L, 1770L, 1770L, 1843L, 
    1843L, 1843L, 1856L, 1856L), next_box_tag = c(42L, 42L, 42L, 
    56L, 56L, 56L, 0L, 0L), dist_mean = c(NA, NA, NA, NA, NA, 
    NA, NA, NA), count = c(3L, 2L, 1L, 0L, 0L, 0L, 0L, 0L), comb = structure(c(2L, 
    2L, 2L, 3L, 3L, 3L, 1L, 1L), .Label = c("6_0", "6_42", "6_56" 
    ), class = "factor"), distance_nest_tag = c(0.189813484, 
    0.189813484, 0.189813484, 0.649465717, 0.649465717, 0.649465717, 
    NA, NA), epp_male_comb = c(NA, NA, NA, NA, NA, NA, NA, NA 
    ), epp_male = c(NA, NA, NA, NA, NA, NA, NA, NA), nest_epmale = c(NA, 
    NA, NA, NA, NA, NA, NA, NA), neighbour = c(0L, 0L, 0L, 0L, 
    0L, 0L, 0L, 0L), ep_nest = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 
    0L), roost_nextbox_tag = c(8L, 8L, 8L, 335L, 335L, 335L, 
    332L, 332L), nestroost = structure(c(3L, 3L, 3L, 2L, 2L, 
    2L, 1L, 1L), .Label = c("6_332", "6_335", "6_8"), class = "factor"), 
    neighbours_roost = c(1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L)), .Names = c("nest_ID", 
"nest_reproductive_period", "tag_ID", "ring_tag_id", "age_tag", 
"p3_tag", "weight_tag", "next_box_tag", "dist_mean", "count", 
"comb", "distance_nest_tag", "epp_male_comb", "epp_male", "nest_epmale", 
"neighbour", "ep_nest", "roost_nextbox_tag", "nestroost", "neighbours_roost" 
), class = "data.frame", row.names = c(NA, -8L))

編輯

這裏是包含0在響應變量的數據集和一個分組變量的小再現的例子。

dat = structure(list(y = c(0, 13.0988072077744, 1.53920266020577, 12.1207857178524, 
33.9470080593601, 0, 0, 3.46572339150589, 1.05917038733605, 14.295924854232, 
11.29307367, 8.1866351980716, 0.0106620794860646, 0.731016675010324, 
5.99548577982932, 0, 0, 3.15404516097624, 7.62275500199758, 0.604545763926581, 
0, 2.37143378704786, 2.39386320579797, 0.800569675164297), x = c(0, 
0, 0, 0, 0, 0, 2.43684896701947, 0.808418724797666, 0.672910050840583, 
0, 0.763154948372394, 1.44573423753027, 0.96113385074772, 0, 
0.498556550480425, 2.43977373047965, 1.55665618954226, 0.88557694968069, 
4.12758995011915, 3.16827587767271, 2.55354765986558, 3.99502024875255, 
1.71053826174466, 3.59887218330055), tag = c(1, 1, 1, 1, 1, 1, 
2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 5, 5, 5, 5, 5, 5)), .Names = c("y", 
"x", "tag"), row.names = c(NA, -24L), class = "data.frame") 

mod2 = glmer(y ~ x + (1|tag), 
      data = dat, family = gaussian(link = log), 
      start = list(theta = 1, coef(lm(y ~ x, data = dat))))

來源

2016-07-13 Marita Ruiz Miñano

'count'是否有零值？ –

是的，它從0 ... –

請提供數據，使其成爲一個可重複的例子。 http://stackoverflow.com/help/mcve感謝您剛剛做出的更新，但數據應採用可重複的格式，例如'dput'，以便在我們的R會話中輕鬆導入，以幫助我們爲您提供幫助。 –

TL;博士爲了使這項工作，你需要設置mustart參數的東西明智的;據報道，與mgcv類似的問題報告爲here。在glmer當前版本，你只需要包括

mustart=pmax(dd$count,1e-3)

在 glmer通話

。

解釋：通常用於GLM（迭代重新加權最小二乘）的迭代算法需要設置一個起點，不僅僅是參數，而且是預測的響應值。高斯模型的默認設置是將預測值設置爲等於觀察值;如果觀測值對於特定鏈路函數是不可能的（例如，當對數變換時零值變爲負 - 無限），則這將破壞東西（即使觀測值是合法的（即，，我們可以觀察到零值甚至當平均約束爲正）因此，gaussian(link="log")家庭項目的$initialize成分是：

if (is.null(etastart) && is.null(start) && is.null(mustart) && 
    ((family$link == "inverse" && any(y == 0)) || (family$link == 
     "log" && any(y <= 0)))) 
    stop("cannot find valid starting values: please specify some") 
mustart <- y

所以，我們真正需要做的是建立mustart所以所有的這個值是正數

然後就是這是否是正確的統計方法的問題，它不是瘋狂的;我t也可能是合理的嘗試泊松模型，或基於log(1+x)的簡單線性模型或序數模型......取決於您對推動訪問的過程以及您觀察到的有關分佈/條件分佈的觀點數量...

來源

2016-07-15 19:39:13

錯誤在glmer：在EVAL錯誤（表達式，ENVIR，enclos）：無法找到有效的開始>值：請指定一些

回答

相關問題