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