Error in the calculation of Lauricella for the case when nulambda[1]<1 and nulambda[p]>1

R/rdstudent.R

@@ -94,7 +94,7 @@ rdstudent <- function(nu1, Sigma1, nu2, Sigma2, bet, eps = 1e-06) {
   } else {
     lauric <- lauricella(delta2, c(rep(0.5, p-1), delta1 + delta2 - p/2), delta1 + delta2,
                          c(1 - lambda[1:(p-1)]/lambda[p], 1/nulambda[p]), eps = eps)
-    lauric <- 1/nulambda[p] * prod(sqrt(nulambda)) * lauric
+    lauric <- (nulambda[p])^(-delta2) * prod(sqrt(nulambda)) * lauric
   }
 
   gamma1p <- gamma((nu1 + p)/2); gamma2p <- gamma((nu2 + p)/2)