# (1) load of data
#
data(endosulfan)
# (2) plot and description of data for non Australian fish in logarithm
#
log10ATV <-log10(subset(endosulfan,(Australian == "no") & (group == "Fish"))$ATV)
plotdist(log10ATV)
descdist(log10ATV,boot=1000)
# (3) fit of a lognormal and a logistic distribution to data
# (classical distributions used for SSD)
# and visual comparison of the fits using a cdf plot
#
fln <- fitdist(log10ATV,"norm")
summary(fln)
fll <- fitdist(log10ATV,"logis")
summary(fll)
cdfcomp(list(fln,fll),legendtext=c("normal","logistic"))
# (4) fit of a logistic distribution
# by minimizing the modified Anderson-Darling AD2L distance
# cf. ?mgedist for definition of this distance
#
fllAD2L <- fitdist(log10ATV,"logis",method="mge",gof="AD2L")
summary(fllAD2L)
plot(fllAD2L)
Run the code above in your browser using DataLab