par(mfrow = c(1, 3)) # 3D plain scatterplot of the density, plot of the density and contour plot scatterplot3d(u[,1], u[,2], pdf_, color="red", main="Density", xlab ="u1", ylab="u2", zlab="dCopula", pch=".") persp(mycopula, dCopula, main ="Density") contour(mycopula, dCopula, xlim = c(0, 1), ylim=c(0, 1), main = "Contour plot") par(mfrow = c(1, 3)) # 3D plain scatterplot of the CDF, plot of the CDF and contour plot scatterplot3d(u[,1], u[,2], cdf, color="red", main="CDF", xlab = "u1", ylab="u2", zlab="pCopula",pch=".") persp(mycopula, pCopula, main = "CDF") contour(mycopula, pCopula, xlim = c(0, 1), ylim=c(0, 1), main = "Contour plot") # 3D plain scatterplot of the multivariate distribution par(mfrow = c(1, 2)) scatterplot3d(v[,1],v[,2], pdf_mvd, color="red", main="Density", xlab = "u1", ylab="u2", zlab="pMvdc",pch=".") scatterplot3d(v[,1],v[,2], cdf_mvd, color="red", main="CDF", xlab = "u1", ylab="u2", zlab="pMvdc",pch=".") persp(multivariate_dist, dMvdc, xlim = c(-4, 4), ylim=c(0, 2), main = "Density") contour(multivariate_dist, dMvdc, xlim = c(-4, 4), ylim=c(0, 2), main = "Contour plot") persp(multivariate_dist, pMvdc, xlim = c(-4, 4), ylim=c(0, 2), main = "CDF") contour(multivariate_dist, pMvdc, xlim = c(-4, 4), ylim=c(0, 2), main = "Contour plot")