From 69d12d4599b693e6a96d0c46c9dc0845f0a6ac32 Mon Sep 17 00:00:00 2001 From: Guillaume Horel Date: Fri, 16 Nov 2018 16:41:59 -0500 Subject: one more special case --- src/lossdistrib.c | 2 ++ 1 file changed, 2 insertions(+) diff --git a/src/lossdistrib.c b/src/lossdistrib.c index 08f7a89..5c70e44 100644 --- a/src/lossdistrib.c +++ b/src/lossdistrib.c @@ -403,6 +403,8 @@ double shockprob(double p, double rho, double Z, int give_log){ return p; } else if (p == 1.) { return give_log ? 0. : 1.; + } else if (p == 0.) { + return give_log ? -INFINITY : 0.; } else { return( pnorm( (qnorm(p, 0, 1, 1, 0) - sqrt(rho) * Z)/sqrt(1 - rho), 0, 1, 1, give_log)); } -- cgit v1.2.3-70-g09d2