/**
* Fills binStats array (second pass through data file).
*
* @param in object providing access to the data
* @throws IOException if an IO error occurs
*/
private void fillBinStats(Object in) throws IOException {
// Load array of bin upper bounds -- evenly spaced from min - max
double min = sampleStats.getMin();
double max = sampleStats.getMax();
double delta = (max - min) / (Double.valueOf(binCount)).doubleValue();
double[] binUpperBounds = new double[binCount];
binUpperBounds[0] = min + delta;
for (int i = 1; i < binCount - 1; i++) {
binUpperBounds[i] = binUpperBounds[i - 1] + delta;
}
binUpperBounds[binCount - 1] = max;
// Initialize binStats ArrayList
if (!binStats.isEmpty()) {
binStats.clear();
}
for (int i = 0; i < binCount; i++) {
SummaryStatistics stats = new SummaryStatistics();
binStats.add(i, stats);
}
// Filling data in binStats Array
DataAdapterFactory aFactory = new DataAdapterFactory();
DataAdapter da = aFactory.getAdapter(in);
try {
da.computeBinStats(min, delta);
} catch (IOException ioe) {
// don't wrap exceptions which are already IOException
throw ioe;
} catch (RuntimeException rte) {
// don't wrap RuntimeExceptions
throw rte;
} catch (Exception e) {
throw MathRuntimeException.createIOException(e);
}
// Assign upperBounds based on bin counts
upperBounds = new double[binCount];
upperBounds[0] = ((double) binStats.get(0).getN()) / (double) sampleStats.getN();
for (int i = 1; i < binCount - 1; i++) {
upperBounds[i] =
upperBounds[i - 1] + ((double) binStats.get(i).getN()) / (double) sampleStats.getN();
}
upperBounds[binCount - 1] = 1.0d;
}
private GeoSegment getQQLineSegment() {
SummaryStatistics stats = new SummaryStatistics();
for (int i = 0; i < sortedData.length; i++) {
stats.addValue(sortedData[i]);
}
double sd = stats.getStandardDeviation();
double mean = stats.getMean();
double min = stats.getMin();
double max = stats.getMax();
// qq line: y = (1/sd)x - mean/sd
GeoPoint startPoint = new GeoPoint(cons);
startPoint.setCoords(min, (min / sd) - mean / sd, 1.0);
GeoPoint endPoint = new GeoPoint(cons);
endPoint.setCoords(max, (max / sd) - mean / sd, 1.0);
GeoSegment seg = new GeoSegment(cons, startPoint, endPoint);
seg.calcLength();
return seg;
}