X-Git-Url: https://git.donarmstrong.com/?p=mothur.git;a=blobdiff_plain;f=simpson.cpp;h=edc35bd6c04f42c3d687dce377df6b6442e36a18;hp=bf3365ddeadd942a53733e85cbb1eeb30cdc0bb7;hb=cf9987b67aa49777a4c91c2d21f96e58bf17aa82;hpb=afed323beddfc099c365104b9071600aea9731bd diff --git a/simpson.cpp b/simpson.cpp index bf3365d..edc35bd 100644 --- a/simpson.cpp +++ b/simpson.cpp @@ -28,13 +28,13 @@ EstOutput Simpson::getValues(SAbundVector* rank){ if(sobs != 0){ double simnum=0.0000; - for(unsigned long int i=1;i<=maxRank;i++){ + for(unsigned long long i=1;i<=maxRank;i++){ simnum += (double)(rank->get(i)*i*(i-1)); } simpson = simnum / (sampled*(sampled-1)); - for(unsigned long int i=1;i<=maxRank;i++){ + for(unsigned long long i=1;i<=maxRank;i++){ double piI = (double) i / (double)sampled; firstTerm += rank->get(i) * pow(piI, 3); secondTerm += rank->get(i) * pow(piI, 2);