From 2d4c32070bfe6572272d99caeb95edb7c24d6085 Mon Sep 17 00:00:00 2001 From: Boyuan Ning <73546421+ningboyuan@users.noreply.github.com> Date: Sat, 12 Mar 2022 23:01:26 +0100 Subject: [PATCH] Update MSMpdfasymmetric.R --- MSMpdfasymmetric/MSMpdfasymmetric.R | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/MSMpdfasymmetric/MSMpdfasymmetric.R b/MSMpdfasymmetric/MSMpdfasymmetric.R index 4a43a51..89fce43 100644 --- a/MSMpdfasymmetric/MSMpdfasymmetric.R +++ b/MSMpdfasymmetric/MSMpdfasymmetric.R @@ -17,12 +17,12 @@ x3 = seq(-1, 3, 0.001) plot(x1, dsn(x1, xi = 0, omega = 1, alpha = -5, log = FALSE), type = "l", lwd = 2.5, ylab = "PDF", xlab = "X", col = "chocolate3", ylim = c(0, 0.9)) -abline(v = median(rsn(x1, xi = 0, omega = 1, alpha = -5))) +abline(v = median(rsn(n=1, xi = 0, omega = 1, alpha = -5))) plot(x2, dsn(x2, xi = 0, omega = 0.55, alpha = 0, log = FALSE), type = "l", lwd = 2.5, ylab = "", xlab = "X", col = "chartreuse4", ylim = c(0, 0.9)) -abline(v = median(rsn(x2, xi = 0, omega = 0.55, alpha = 0))) +abline(v = median(rsn(n=1, xi = 0, omega = 0.55, alpha = 0))) plot(x3, dsn(x3, xi = 0, omega = 1, alpha = 5, log = FALSE), type = "l", lwd = 2.5, ylab = "", xlab = "X", col = "blue3", ylim = c(0, 0.9)) -abline(v = median(rsn(x3, xi = 0, omega = 1, alpha = 5))) +abline(v = median(rsn(n=1, xi = 0, omega = 1, alpha = 5)))