[R-sig-ME] error: model is nearly unidentifiable

Ligia Pizzatto do Prado ligia_oceanica at hotmail.com
Mon Apr 9 08:14:09 CEST 2018


Hi there, I'm new to mixed models but have ran a few with success before. Now, while trying to analyse this new experiment I am  having an error that I quite don't understand...


The experiment is a two choice habitat ("choice": poor [0] vs rich [1]) for frogs under two-state treatments, lets say F and C. Then I have as potential variables frog size ("size"), air temperature ("temp"), humidity ("hum") and date of experiment (recorded as continuos variable starting at day 1...). This is a repeated measure design as frogs were tested both in F and C trials (thus id is my random effect). I want to know if the choice is affected by treat, but also considering size, temp, humidity, and date in my model.


First I did:


data$treat<- factor(data$treat)

data$id<- factor(data$id)

data$choice<- factor(data$choice)


summary(data)

treat        id     choice

C:24   1      : 2   0:24

F:24   2      : 2   1:24

           3      : 2

           4      : 2

           5      : 2

           6      : 2

      (Other):36



size                     temp                      hum                   date

 Min.   :35.70      Min.   :24.80       Min.   :53.00       Min.   : 1.00

 1st Qu.:38.25     1st Qu.:26.30     1st Qu.:58.50     1st Qu.: 4.00

 Median :42.40   Median :27.30   Median :63.00   Median : 9.00

 Mean   :42.02    Mean   :26.74     Mean   :63.65    Mean   :10.75

 3rd Qu.:44.27    3rd Qu.:27.40     3rd Qu.:70.50    3rd Qu.:16.75

 Max.   :51.00      Max.   :27.90     Max.   :76.00       Max.   :24.00

m1<- glmer(Fchoice ~ treat + SUL + temp + hum + date + (1|id), data = data, family = binomial)

Warning message:
In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv,  :
  Model is nearly unidentifiable: large eigenvalue ratio
 - Rescale variables?


I keep getting this message in all models except when I exclude both temp and hum, but I also get the message when tried null model: null<- glmer(choice ~ 1 + (1|id), data = data, family = binomial)


I tried to transform/re-scale all continuous variable (temp, hum, date) and nothing changed, and I quite don't understand why the error also appears in the null model, given id is a factor... If its a scale problem wouldn't this only appear in the continuous variables?


Can anyone provide some guidance here, please?


TIA,


Ligia


	[[alternative HTML version deleted]]



More information about the R-sig-mixed-models mailing list