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

Trivial bug fix 

git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@4848 348d0f76-0448-11de-a6fe-93d51630548a
This commit is contained in:
fromer 2010-12-15 18:29:15 +00:00
parent dabdeb729e
commit aea481ae01
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

@ -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. # 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) { pHetPairAtDistance <- function(k, theta) {
pOneSiteIsHetTheta = pOneSiteIsHet(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: # Since the geometric/exponential distribution is "memory-free", can simply multiply the (independent) probabilities for the distances: