INNER CODE UNIT · C

avg

aquaskyline/SOAPdenovo2 · fusion/attachPEinfo.c:484

              long long int avg = sum / weight;
              sum = 0;

              for(ii = 0; ii < weight; ii++)
                {
                  sum += ((tmp->PE[ii] - avg) * (tmp->PE[ii] - avg));
                }

              double SD = (sqrt((double)sum / (weight - 1))) * 3; //just for fast
              sum = 0;
              int num = 0;

              for(ii = 0; ii < weight; ii++)
                {
                  if(abs(tmp->PE[ii] - avg) <= SD)
                    {
                      sum += tmp->PE[ii];
                      num++;

View source record →