weighted jackknife addons

pievos101 · Apr 6, 2018 · 133acd7 · 133acd7
1 parent 200093b
commit 133acd7
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 1 deletion.
diff --git a/R/GENOME.R b/R/GENOME.R
@@ -129,7 +129,8 @@ RNDmin 		       =   "matrix",
 D.z		       =   "matrix",              # introgression slots
 D.pval		       =   "matrix",
 BDF.z		       =   "matrix",              # introgression slots
-BDF.pval	       =   "matrix",                            
+BDF.pval	       =   "matrix",
+BDF.SE		       =   "matrix",                            
 jack.knife             =   "logical",
 missing.freqs          =   "matrix",
 

diff --git a/R/weighted.jackknife.R b/R/weighted.jackknife.R
@@ -20,6 +20,7 @@ if(per.region){ #drop-one jackknife
   D.pval   <- matrix(NaN,n.region.names,1) # jacknife
   BDF.z    <- matrix(NaN,n.region.names,1) # jacknife
   BDF.pval <- matrix(NaN,n.region.names,1) # jacknife
+  BDF.SE   <- matrix(NaN,n.region.names,1) # jacknife
 
   # get site specific values 
   site.D    <- object@region.stats@D
@@ -149,6 +150,7 @@ object@BDF.z    <- as.matrix(BDF/SE_BDF)
 object@BDF.pval <- as.matrix(2*pnorm(-abs(object@BDF.z)))
 object@D.z      <- as.matrix(D/SE_D)
 object@D.pval   <- as.matrix(2*pnorm(-abs(object@D.z)))
+object@BDF.SE   <- as.matrix(SE_BDF)
 
 return(object)