Skip to content
Browse files

bug in D_Jost gave wrong harmonic mean

  • Loading branch information...
1 parent 477c9d8 commit c4362ab0dcd4f30df3c866b938bddd6ae4da95d6 @dwinter committed Feb 1, 2012
Showing with 1 addition and 1 deletion.
  1. +1 −1 R/D_Jost.R
  2. BIN mmod_0.2.tar.gz
View
2 R/D_Jost.R
@@ -38,7 +38,7 @@ D_Jost <- function(x){
global_Hs <- mean(loci[,1], na.rm=T)
global_Ht <- mean(loci[,2], na.rm=T)
global_D <- (global_Ht - global_Hs)/(1 - global_Hs ) * (n/(n-1))
- harm_D <- harmonic_mean(loci)
+ harm_D <- harmonic_mean(loci[,3])
return(list("per.locus"=loci[,3],
"global.het"=global_D,
"global.harm_mean" = harm_D
View
BIN mmod_0.2.tar.gz
Binary file not shown.

0 comments on commit c4362ab

Please sign in to comment.
Something went wrong with that request. Please try again.