diff options
| -rw-r--r-- | src/lossdistrib.c | 2 |
1 files changed, 2 insertions, 0 deletions
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)); } |