From d895d6279956861b273adcbeb3bde31d57dc400f Mon Sep 17 00:00:00 2001
From: dmalizia <48389764+dmalizia@users.noreply.github.com>
Date: Fri, 15 May 2020 18:12:03 +0100
Subject: [PATCH] Update sn-funct.R

---
 R/sn-funct.R | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/R/sn-funct.R b/R/sn-funct.R
index 306aa43..0674b3d 100644
--- a/R/sn-funct.R
+++ b/R/sn-funct.R
@@ -258,7 +258,7 @@ pmsn <- function(x, xi=rep(0,length(alpha)), Omega, alpha, tau=0,
   Obig <- matrix(rbind(c(1,-delta), cbind(-delta,Ocor)), d+1, d+1)
   x <- if (is.vector(x)) matrix(x, 1, d) else data.matrix(x)
   if (is.vector(xi)) xi <- outer(rep(1, nrow(x)), as.vector(matrix(xi,1,d)))
-  z0 <- cbind(tau, t(t(x - xi))/omega) 
+  z0 <- cbind(tau, t(t(x - xi)/omega))
   mnormt::pmnorm(z0, mean=rep(0,d+1), varcov=Obig, ...)/pnorm(tau) 
 }