## First fitting model without random effects
model1 <- multdrc(SLOPE~DOSE, CURVE, data=PestSci,
collapse=data.frame(HERBICIDE, 1, 1, HERBICIDE))
## Then fitting model with random effects (assigned to d)
model2<-mixdrc(model1, random="d~1|CURVE", data=PestSci)
compParm(model2, "e", "-")
rm(model1, model2)
Run the code above in your browser using DataLab