Displaying 1 result from an estimated 1 matches for "treatslop".
Did you mean:
treatslope
2007 Feb 26
0
LD50 contrasts with lmer/lme4
...40%
with B at dose xx" has less meaning than giving "equivalent effective
doses".
Dieter
----- Simulated data -----
library(lme4)
animal = data.frame(ID = as.factor(1:20), da = rnorm(1:20))
treat = data.frame(treat=c('A','B','C'), treatoff=c(1,2,1.5),
treatslope = c(0.5,0.6,0.7))
gr = expand.grid(animal=animal$ID,treat=treat$treat,logdose=c(-4:4))
gr$resp = as.integer(treat$treatoff[gr$treat]+
treat$treatslope[gr$treat]*gr$logdose+
animal$da[gr$animal] + rnorm(nrow(gr),0,2) >0)
gr.lmer = lmer(resp ~ treat*logdose+(1|animal),data=gr,family=binomi...