1

I am unable to include more than 2 fixed effects in my LME without getting an error.

When I include a variable called DD in my fixed effects. I get an error. R recognizes my previous model - which is a simplified model without DD.

I know DD is the problem because when I take it off, the model works perfectly. But I know that for my clinical data DD is an important fixed effect that I need to account for. This is my model:

  mod1.lme <- lme(Variable1 ~ Age + Sex + DD + Category, 
random= ~1|Patient, data = AFN1, na.action = na.omit)
    summary(mod1.lme)

Error in MEEM(object, conLin, control$niterEM) :Singularity in backsolve at level 0, block 1

Data:

Patient Variable1  Age  Sex  DD  Category

   1a       12      30    0    346   A
   1b       22      30    0    346   F
   2a       29      67    0    700   A
   2b       12      67    0    700   A
   3a       50      19    1    36    N
   3b       54      19    1    36    N
   4a       44      36    1    55    F
   4b       22      36    1    55    A

Could this happen because Disease duration and another variable (maybe Age?) are colinear in my model? Could be LME warning me about the confounded fixed effects... I am not sure if this statement would mean that I can confidently drop DD from my model or would mean something a bit more serious.

I would really appreciate your help on this.

Best regards, Lili

NelsonGon
  • 13,015
  • 7
  • 27
  • 57
Lili
  • 547
  • 6
  • 19

0 Answers0