model1 <- multdrc(rgr~dose, data=TerMet)
diagnostics(model1)
model2 <- multdrc(rgr~dose, data=TerMet, startVal=c(1,0.01,0.38,130))
diagnostics(model2)
## the start values given yield swifter convergence
rm(model1, model2)
Run the code above in your browser using DataLab