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: