double clusterCutoff = cutoff;
while (matrix->getSmallDist() <= clusterCutoff && matrix->getNNodes() > 0){
cluster->update(clusterCutoff);
- float dist = matrix->getSmallDist();
}
list->setLabel(toString(cutoff));
for(int j=0;j<nSeqsPerOTU[i];j++){
int index = cumNumSeqs[i] + j;
int nI = seqIndex[index];
- int nIU = mapSeqToUnique[nI];
double tauValue = singleTau[seqNumber[index]];
for(int j=0;j<nSeqsPerOTU[i];j++){
int index = cumNumSeqs[i] + j;
- int nI = seqIndex[index];
double tauValue = singleTau[seqNumber[index]];
weight[i] += tauValue;
}
vector<vector<int> > qualities(numOTUs);
vector<double> pr(HOMOPS, 0);
- int index = 0;
for(int i=0;i<numOTUs;i++){
int index = 0;