for(int i=0;i<NUMBINS;i++){
for(int j=0;j<NUMBINS;j++){
- float minSum = 10000000000;
+ double minSum = 100000000;
for(int k=0;k<HOMOPS;k++){
- float sum = singleLookUp[k * NUMBINS + i] + singleLookUp[k * NUMBINS + j];
+ double sum = singleLookUp[k * NUMBINS + i] + singleLookUp[k * NUMBINS + j];
if(sum < minSum) { minSum = sum; }
}
/**************************************************************************************************/
-float ShhherCommand::getProbIntensity(int intIntensity){
+double ShhherCommand::getProbIntensity(int intIntensity){
try{
- float minNegLogProb = 10000000000;
+
+ double minNegLogProb = 100000000;
+
for(int i=0;i<HOMOPS;i++){//loop signal strength
float negLogProb = singleLookUp[i * NUMBINS + intIntensity];
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;