df: l.smooth=TRUE

pievos101 · Apr 11, 2020 · b5e5a74 · b5e5a74
1 parent 9a09395
commit b5e5a74
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/R/calc_df.R b/R/calc_df.R
@@ -186,7 +186,7 @@ L_baba      <- lbeta(alpha_BABA,beta_BBAA) #(lgamma(alpha_BABA)*lgamma(beta_BBAA
 D_bayes             <- D
 D_bayes[D==0]       <- 1
 D_bayes[D >0]       <- 1 + exp((L_abba/L_baba)) - exp(1)
-D_bayes[D <0]       <- 1 + exp((L_baba/L_abba)) - exp(1)
+D_bayes[D <0]       <- -(1 + exp((L_baba/L_abba)) - exp(1))
 
 #D_bayes            <- sum_ABBA/sum_BABA
 

diff --git a/R/introgression.stats.R b/R/introgression.stats.R
@@ -1,4 +1,4 @@
-setGeneric("introgression.stats", function(object, subsites=FALSE, do.D=TRUE, do.df=TRUE, keep.site.info=TRUE, block.size=FALSE, do.RNDmin=FALSE, l.smooth=FALSE) standardGeneric("introgression.stats"))
+setGeneric("introgression.stats", function(object, subsites=FALSE, do.D=TRUE, do.df=TRUE, keep.site.info=TRUE, block.size=FALSE, do.RNDmin=FALSE, l.smooth=TRUE) standardGeneric("introgression.stats"))
 
 setMethod("introgression.stats","GENOME",function(object, subsites, do.D, do.df, keep.site.info, block.size, do.RNDmin, l.smooth){