From 83124e719f6bd7cf33a24786e248a7d677065484 Mon Sep 17 00:00:00 2001 From: Filipe Rodrigues Date: Sat, 17 Jun 2023 15:31:31 +0100 Subject: [PATCH] Fixed 8. --- 8/script.R | 7 ++++--- 1 file changed, 4 insertions(+), 3 deletions(-) diff --git a/8/script.R b/8/script.R index 17af68d..22a6fb2 100755 --- a/8/script.R +++ b/8/script.R @@ -15,13 +15,13 @@ cauchy_gen <- sort(cauchy_gen) quantile_idxs <- 1:107 / (107 + 1) cauchy_quantiles <- sapply(quantile_idxs, function(p) { - cauchy_location + cauchy_scale * tan(pi * (p - 0.5)) + qcauchy(p, location = cauchy_location, scale = cauchy_scale) }) norm_mean <- -2.8 norm_var <- 1.4 norm_quantiles <- sapply(quantile_idxs, function(p) { - norm_mean + sqrt(2 * norm_var) * erfinv(2 * p - 1) + qnorm(p, mean = norm_mean, sd = sqrt(norm_var)) }) df <- data.frame( @@ -36,9 +36,10 @@ plot <- ggplot(df, aes(x = quantile_idxs)) + geom_line(aes(y = cauchy_gen, color = "Generated")) + geom_line(aes(y = cauchy_quantiles, color = "Cauchy Quantiles")) + geom_line(aes(y = norm_quantiles, color = "Normal Quantiles")) + + geom_line(aes(y = quantile_idxs, color = "Line bisector of odd quadrants")) + xlab("Quantiles") + ylab("") + - scale_color_manual(values = c("red", "blue", "green")) + + scale_color_manual(values = c("red", "blue", "black", "green")) + theme(legend.title = element_blank()) ggsave(plot = plot, "output.png")