/** * Creates the recalibration report. Report can then be written to a stream via * GATKReport.print(PrintStream). * * @return newly created recalibration report */ public GATKReport createGATKReport() { return RecalUtils.createRecalibrationGATKReport( argumentTable, quantizationInfo, recalibrationTables, requestedCovariates, RAC.SORT_BY_ALL_COLUMNS); }
/** * Combines two recalibration reports by adding all observations and errors * * <p>Note: This method DOES NOT recalculate the empirical qualities and quantized qualities. You * have to recalculate them after combining. The reason for not calculating it is because this * function is intended for combining a series of recalibration reports, and it only makes sense * to calculate the empirical qualities and quantized qualities after all the recalibration * reports have been combined. Having the user recalculate when appropriate, makes this method * faster * * <p>Note2: The empirical quality reported, however, is recalculated given its simplicity. * * @param other the recalibration report to combine with this one */ public void combine(final RecalibrationReport other) { for (int tableIndex = 0; tableIndex < recalibrationTables.numTables(); tableIndex++) { final NestedIntegerArray<RecalDatum> myTable = recalibrationTables.getTable(tableIndex); final NestedIntegerArray<RecalDatum> otherTable = other.recalibrationTables.getTable(tableIndex); RecalUtils.combineTables(myTable, otherTable); } }
public RecalibrationReport(final File recalFile, final SortedSet<String> allReadGroups) { final GATKReport report = new GATKReport(recalFile); argumentTable = report.getTable(RecalUtils.ARGUMENT_REPORT_TABLE_TITLE); RAC = initializeArgumentCollectionTable(argumentTable); GATKReportTable quantizedTable = report.getTable(RecalUtils.QUANTIZED_REPORT_TABLE_TITLE); quantizationInfo = initializeQuantizationTable(quantizedTable); Pair<ArrayList<Covariate>, ArrayList<Covariate>> covariates = RecalUtils.initializeCovariates(RAC); // initialize the required and optional covariates ArrayList<Covariate> requiredCovariates = covariates.getFirst(); ArrayList<Covariate> optionalCovariates = covariates.getSecond(); requestedCovariates = new Covariate[requiredCovariates.size() + optionalCovariates.size()]; optionalCovariateIndexes = new HashMap<String, Integer>(optionalCovariates.size()); int covariateIndex = 0; for (final Covariate covariate : requiredCovariates) requestedCovariates[covariateIndex++] = covariate; for (final Covariate covariate : optionalCovariates) { requestedCovariates[covariateIndex] = covariate; final String covariateName = covariate .getClass() .getSimpleName() .split("Covariate")[ 0]; // get the name of the covariate (without the "covariate" part of it) so we can // match with the GATKReport optionalCovariateIndexes.put(covariateName, covariateIndex - 2); covariateIndex++; } for (Covariate cov : requestedCovariates) cov.initialize( RAC); // initialize any covariate member variables using the shared argument collection recalibrationTables = new RecalibrationTables(requestedCovariates, allReadGroups.size()); initializeReadGroupCovariates(allReadGroups); parseReadGroupTable( report.getTable(RecalUtils.READGROUP_REPORT_TABLE_TITLE), recalibrationTables.getReadGroupTable()); parseQualityScoreTable( report.getTable(RecalUtils.QUALITY_SCORE_REPORT_TABLE_TITLE), recalibrationTables.getQualityScoreTable()); parseAllCovariatesTable( report.getTable(RecalUtils.ALL_COVARIATES_REPORT_TABLE_TITLE), recalibrationTables); }