diff options
| author | Guillaume Horel <guillaume.horel@gmail.com> | 2018-11-16 16:41:59 -0500 |
|---|---|---|
| committer | Guillaume Horel <guillaume.horel@gmail.com> | 2018-11-16 16:41:59 -0500 |
| commit | 69d12d4599b693e6a96d0c46c9dc0845f0a6ac32 (patch) | |
| tree | e037c78a8f8a60d2454f47adcc484d2ea6e81f44 /src | |
| parent | 62545be926e3ed2223620ad78dcb1e7b532a12d0 (diff) | |
| download | lossdistrib-69d12d4599b693e6a96d0c46c9dc0845f0a6ac32.tar.gz | |
one more special case
Diffstat (limited to 'src')
| -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)); } |