Commit 55049d55 authored by Victor's avatar Victor
Browse files

metrics bug fixed

parent cb0a09c7
Pipeline #74698 passed with stage
in 23 minutes and 36 seconds
......@@ -124,11 +124,12 @@ function get_beta_div(world::World,trait=1)
if trait == 0
# need to convert to Float64, otherwise infinite variance
sbar_i = [mean(Float64.(get_geo(World(subw,space(world),parameters(world))))) for subw in values(g)]
return var(sbar_i,corrected=false)
m = [mean(World(subw,space(world),parameters(world)),trait=trait) for subw in values(g)]
return mean(var(h,dims=1,corrected=false))
return mean(var(h,dims=1,corrected=false))
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment