From 2bf4fc94f09c350b239fe4b67ddfe6ef34715d2b Mon Sep 17 00:00:00 2001 From: fromer Date: Fri, 10 Dec 2010 02:58:46 +0000 Subject: [PATCH] Try to use more sampling to get a "correct" estimate of multivariate integral git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@4815 348d0f76-0448-11de-a6fe-93d51630548a --- R/phasing/RBP_theoretical.R | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/R/phasing/RBP_theoretical.R b/R/phasing/RBP_theoretical.R index 334c7f153..cf9d8c382 100644 --- a/R/phasing/RBP_theoretical.R +++ b/R/phasing/RBP_theoretical.R @@ -230,7 +230,7 @@ pDirectlyPhaseHetPairUsingWindow <- function(meanDepth, nReadsToPhase, L, theta, lower = as.vector(matrix(data=MIN_DISTANCE, nrow=1, ncol=ndim)) upper = as.vector(matrix(data=MAX_DISTANCE, nrow=1, ncol=ndim)) - N = 10^4 * ndim + N = 10^4 * ndim^2 high_dimensional_integrate(ndim, lower, upper, integrandFunction, N, DEBUG = TRUE, PRINT_EVERY = 10^2) }