data(models_dawson)
plot(models_dawson[,1])
lines(models_dawson[,2])
diagnostic_dawson(modelled=models_dawson[,2], measured=models_dawson[,1], use_qualV=TRUE)
do.call(rbind,lapply(models_dawson,FUN=diagnostic_dawson, measured=models_dawson[,1]))
Run the code above in your browser using DataLab