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)))