model1 <- multdrc(FA)
coef(model1)
model2 <- multdrc(FA, boxcox=TRUE)
summary(model2) # not a satisfactory fit
## Usinf estimates from 'model1' as start values
model3 <- multdrc(FA, boxcox=TRUE, startVal=coef(model1))
summary(model3) # a satisfactory fit
rm(model1, model2, model3)
Run the code above in your browser using DataLab