From 32f0a614798f42d0ff51a325ba183aaa85030b24 Mon Sep 17 00:00:00 2001 From: Justin Bedo Date: Wed, 29 Mar 2017 09:14:33 +1100 Subject: Correct p-values --- R/pmt.R | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/R/pmt.R b/R/pmt.R index 5921638..927fd9e 100644 --- a/R/pmt.R +++ b/R/pmt.R @@ -54,7 +54,9 @@ pmt <- function(a, b, two.sided=T, signed=F) { t1 <- mapply(function(xa, mu, rho) { i <- 0:xa; logsumexp(i * log(2 * mu) - lfactorial(i) - i * log(1 + rho))}, a, mu, rho) t2 <- mapply(function(xb, mu, rho) { j <- xb:(xb*100); logsumexp(j * log(2 * mu * rho) - lfactorial(j) - j * log(1 + rho))}, b, mu, rho) - ps <- (2*mu-t1-t2) / log(10) + ps <- (2*mu-t1-t2) + ps <- ps - log1p(ps) + ps <- ps / log(10) if(!two.sided) ps[msk] <- NA else if(signed) ps <- ps * (1-2*msk) ps -- cgit v1.2.3