diff --git a/README.md b/README.md
index 8522ef2..ca7bd2a 100644
--- a/README.md
+++ b/README.md
@@ -3,7 +3,7 @@
 Simulated data examples for an article about multivariate extreme quantile
 region estimation. Simulations produce figures of the following form.
 
-![image](figures/fig-n_5000-k_400-p_low-gamma_1.00.png)
+![image](figures/fig-n_5000-k_400-p_medium-gamma_1.00.jpg)
 
 In above figure $Q_p$ denotes the true $(1-p)$-quantile region, $\bar Q_p$
 denotes estimate based on sample quantile and $\hat Q_p$ denotes estimate based
diff --git a/figures/fig-n_5000-k_400-p_high-gamma_1.00.jpg b/figures/fig-n_5000-k_400-p_high-gamma_1.00.jpg
new file mode 100644
index 0000000..04c75ed
Binary files /dev/null and b/figures/fig-n_5000-k_400-p_high-gamma_1.00.jpg differ
diff --git a/figures/fig-n_5000-k_400-p_high-gamma_1.00.png b/figures/fig-n_5000-k_400-p_high-gamma_1.00.png
deleted file mode 100644
index 9112502..0000000
Binary files a/figures/fig-n_5000-k_400-p_high-gamma_1.00.png and /dev/null differ
diff --git a/figures/fig-n_5000-k_400-p_low-gamma_1.00.jpg b/figures/fig-n_5000-k_400-p_low-gamma_1.00.jpg
new file mode 100644
index 0000000..9b6d920
Binary files /dev/null and b/figures/fig-n_5000-k_400-p_low-gamma_1.00.jpg differ
diff --git a/figures/fig-n_5000-k_400-p_low-gamma_1.00.png b/figures/fig-n_5000-k_400-p_low-gamma_1.00.png
deleted file mode 100644
index 690b663..0000000
Binary files a/figures/fig-n_5000-k_400-p_low-gamma_1.00.png and /dev/null differ
diff --git a/figures/fig-n_5000-k_400-p_medium-gamma_1.00.jpg b/figures/fig-n_5000-k_400-p_medium-gamma_1.00.jpg
new file mode 100644
index 0000000..5e967bf
Binary files /dev/null and b/figures/fig-n_5000-k_400-p_medium-gamma_1.00.jpg differ
diff --git a/figures/fig-n_5000-k_400-p_medium-gamma_1.00.png b/figures/fig-n_5000-k_400-p_medium-gamma_1.00.png
deleted file mode 100644
index ae6907a..0000000
Binary files a/figures/fig-n_5000-k_400-p_medium-gamma_1.00.png and /dev/null differ
diff --git a/plot.R b/plot.R
index cbf9db4..15eaf7b 100644
--- a/plot.R
+++ b/plot.R
@@ -71,7 +71,7 @@ plot_ellipse <- function(tsample, q, qbar, qhat) {
 
   # Plotting
   dir.create("figures", showWarnings = FALSE)
-  f <- sprintf("%s/fig-n_%d-k_%d-p_%s-gamma_%.2f.png", "figures", opt$n,
+  f <- sprintf("%s/fig-n_%d-k_%d-p_%s-gamma_%.2f.jpg", "figures", opt$n,
                opt$k, opt$p, opt$gamma)
   ggplot(data) +
     geom_ellipse(aes(x0 = x, y0 = y, a = a, b = b, angle = theta,
@@ -92,7 +92,7 @@ plot_ellipse <- function(tsample, q, qbar, qhat) {
                           labels = c(TeX("$Q_p$"),
                                      TeX("$\\hat{Q}_p$"),
                                      TeX("$\\bar{Q}_p$")))
-  ggsave(f, width = 7, height = 7)
+  ggsave(f, width = 7, height = 7, dpi = 1000)
 }