[R-sig-ME] completely nested glmer
Troels Ring
tring at gvdnet.dk
Wed Aug 30 10:15:33 CEST 2017
Dear friends - I have a question which is basic and probably simple but
disturbing to me. I am on Windows 7 or 10, R version 3.4.1 (2017-06-30)
-- "Single Candle"
I have 2013 measurements of capillary flow speed obtained from 204
glomeruli originating from 29 rats, 10 of whom are controls, 11 made
diabetic, and 8 made hyperglycemic otherwise. Altogether I have 621
capillaries from controls, 964 from diabetics and 428 from hyperglycemic.
If I make a direct aov
summary(z2 <- aov(Speed~TRT,SCAN))
TukeyHSD(z2)
#$TRT
# diff lwr upr p adj
#Diab-Ctrl -0.4266708 -0.65092625 -0.20241526 0.0000255
#Hyper-Ctrl -0.2206938 -0.49449343 0.05310581 0.1416500
#Hyper-Diab 0.2059769 -0.04716959 0.45912348 0.1365243
and a similar result is obtained with kruskal.test
I happen to know the flow is approximately gamma distributed and
summary(z1a <- glm(Speed~TRT,SCAN,family=Gamma)) where TRT is the three
level factor, leaves me a significant effect of diabetes
However, as it is, glomeruli are nested within rats, and rats are by
design nested within one of the three groups
Hence I make
summary(z3a <- glmer(Speed~TRT + (1|RAT)+(1|ind) + (1|RAT:ind)
,data=SCAN,family=Gamma)) where ind is an indicator for each of the 204
glomeruli. And the fixed effect TRT is not significant.
But taking into account that RAT is nested within TRT leaves me
something like
summary(z77 <- glmer(Speed~ 1 +(1|TRT/RAT/ind) ,data=SCAN,family=Gamma))
which fits without protests and gives me
Random effects:
Groups Name Variance Std.Dev.
ind:(RAT:TRT) (Intercept) 0.003839 0.06196
RAT:TRT (Intercept) 0.001479 0.03846
TRT (Intercept) 0.000000 0.00000
Residual 0.237689 0.48753
Number of obs: 2013, groups: ind:(RAT:TRT), 204; RAT:TRT, 29; TRT, 3
which hardly indicates strong effect of TRT levels?
The results of the aov or first glm would fit expectations - but ignores
the nesting. So is the mixed model wrong? (it looks reasonable from
diagnostic plots)
Best wishes
Troels Ring
Aalborg, Denmark
More information about the R-sig-mixed-models
mailing list