From aea481ae01610e2a21283a8df8b3953484b4e4dc Mon Sep 17 00:00:00 2001 From: fromer Date: Wed, 15 Dec 2010 18:29:15 +0000 Subject: [PATCH] Trivial bug fix git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@4848 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 e93550b33..371b2d1e0 100644 --- a/R/phasing/RBP_theoretical.R +++ b/R/phasing/RBP_theoretical.R @@ -26,7 +26,7 @@ meanIntraHetDistanceToTheta <- function(d) { # That is, het site x already "exists", and we want to know what the probability that the NEXT het site (y) is k bases away. pHetPairAtDistance <- function(k, theta) { pOneSiteIsHetTheta = pOneSiteIsHet(theta) - dexp(k, pOneSiteIsHet) + dexp(k, pOneSiteIsHetTheta) } # Since the geometric/exponential distribution is "memory-free", can simply multiply the (independent) probabilities for the distances: