## Not run: par(mfrow = c(1, 3))
# ymat <- rbilogis(n = 2000, loc1 = 5, loc2 = 7, scale2 = exp(1))
# myxlim <- c(-2, 15); myylim <- c(-10, 30)
# plot(ymat, xlim = myxlim, ylim = myylim)
#
# N <- 100
# x1 <- seq(myxlim[1], myxlim[2], len = N)
# x2 <- seq(myylim[1], myylim[2], len = N)
# ox <- expand.grid(x1, x2)
# z <- dbilogis(ox[,1], ox[,2], loc1 = 5, loc2 = 7, scale2 = exp(1))
# contour(x1, x2, matrix(z, N, N), main = "density")
# z <- pbilogis(ox[,1], ox[,2], loc1 = 5, loc2 = 7, scale2 = exp(1))
# contour(x1, x2, matrix(z, N, N), main = "cdf") ## End(Not run)
Run the code above in your browser using DataLab