|
1 | 1 | require(fitdistrplus)
|
2 | 2 | data(salinity)
|
3 | 3 | log10LC50 <-log10(salinity)
|
4 |
| -fit <- fitdistcens(log10LC50, "norm") |
| 4 | +fitn <- fitdistcens(log10LC50, "norm") |
| 5 | +fitc <- fitdistcens(log10LC50, "cauchy") |
| 6 | +cdfcompcens(list(fitc, fitn), xlim = c(1, 1.8), main = "", |
| 7 | + fitcol = c("blue", "orange"), |
| 8 | + fitlty = c(1, 1), fitlwd = c(2, 2)) |
| 9 | + |
| 10 | +data(salinity) |
| 11 | +fitln <- fitdistcens(salinity, "lnorm") |
| 12 | +fitw <- fitdistcens(salinity, "weibull") |
| 13 | +fitg <- fitdistcens(salinity, "gamma") |
| 14 | +cdfcompcens(list(fitln, fitw, fitg), main = "", |
| 15 | + fitcol = c("blue", "orange", "green"), |
| 16 | + fitlty = c(1, 1, 1), fitlwd = c(2, 2, 2)) |
| 17 | +cdfcompcens(list(fitln, fitw, fitg), main = "", |
| 18 | + fitcol = c("blue", "orange", "green"), |
| 19 | + fitlty = c(1, 1, 1), fitlwd = c(2, 2, 2), |
| 20 | + xlogscale = TRUE) |
0 commit comments