par(mfrow = c(1, 3))
n = 2000
ymat = rbilogis4(n, 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 = dbilogis4(ox[,1], ox[,2], loc1 = 5, loc2 = 7, scale2 = exp(1))
contour(x1, x2, matrix(z, N, N), main = "density")
z = pbilogis4(ox[,1], ox[,2], loc1 = 5, loc2 = 7, scale2 = exp(1))
contour(x1, x2, matrix(z, N, N), main = "cdf")
Run the code above in your browser using DataLab