zzh
/
gatk-3.8
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

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
This commit is contained in:
fromer 2010-12-10 02:58:46 +00:00
parent c64bf80b57
commit 2bf4fc94f0
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
R/phasing/RBP_theoretical.R
View File

@ -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)
}