public static void setBaselineComparison(ProteinProphetResult ppr, int numMS2, int ddaNum)
{
if (!baselineComparisonSet.ContainsKey(ddaNum))
{
BaselineComparison bc = new BaselineComparison(ppr, numMS2, ddaNum);
baselineComparisonSet.Add(ddaNum, bc);
}
}
private void postProcessingCalculations(int ddaNum, ProteinProphetResult ppr, ExclusionProfile exclusionProfile)
{
BaselineComparison bc = baselineComparisonSet[ddaNum];
List <String> proteinsIdentifiedByNoExclusion = bc.getProteinsIdentifiedByNoExclusion();
int totalResourcesNaiveExperiment = bc.getTotalResourcesNaiveExperiment();
int numProteinsIdentifiedNaiveExperiment = bc.getNumProteinsIdentifiedNaiveExperiment();
// set proteins identified first
setProteinsIdentified(ppr, proteinsIdentifiedByNoExclusion);
int correctlyExcluded = (int)data[Header.EvaluateExclusion_FoundOnCurrentExclusionList]
+ (int)data[Header.EvaluateExclusion_FoundOnCurrentObservedExclusionList];
int incorrectlyExcluded = (int)data[Header.EvaluateExclusion_NotFoundOnExclusionList];
/*
* found on past observed, found on past exclusion list, and found on future
* exclusion list are not incorrect exclusions, they are retention time being
* predicted incorrectly...
*/
double ratioIncorrectlyExcludedOverCorrectlyExcluded = takeRatio(incorrectlyExcluded, correctlyExcluded);
ChangeValue(Header.CorrectlyExcluded, correctlyExcluded);
ChangeValue(Header.IncorrectlyExcluded, incorrectlyExcluded);
ChangeValue(Header.RatioIncorrectlyExcludedOverCorrectlyExcluded, ratioIncorrectlyExcludedOverCorrectlyExcluded);
// Resources saved in total # available MS2 - # ms2 used foreach analysis
int resourcesSaved = totalResourcesNaiveExperiment - (int)data[Header.NumMS2Analyzed];
double percentResourcesSaved = takeRatio(resourcesSaved, totalResourcesNaiveExperiment);
double percentResourcesUsed = 1 - percentResourcesSaved;
ChangeValue(Header.PercentResourcesSaved, percentResourcesSaved);
ChangeValue(Header.PercentResourcesUsed, percentResourcesUsed);
/*-
* Protein Identification Sensitivity = # proteins identified / # proteins identified in whole experiment
* Protein Identification Fold Change = # proteins identified / # proteins identified by naive approach
* Protein Identification Sensitivity Limited DDA = # proteins identified also identified in naive approach / proteins identified by naive approach
*/
ChangeValue(Header.ProteinIdentificationSensitivity,
takeRatio((int)data[Header.NumProteinsIdentified], numProteinsIdentifiedOriginalExperiment));
ChangeValue(Header.ProteinIdentificationFoldChange,
takeRatio((int)data[Header.NumProteinsIdentified], numProteinsIdentifiedNaiveExperiment));
ChangeValue(Header.ProteinIdentificationSensitivityLimitedDDA,
takeRatio((int)data[Header.ProteinsIdentifiedInLimitedDDA], numProteinsIdentifiedNaiveExperiment));
List <String> inProgramExcludedProteins = exclusionProfile.getDatabase().getExcludedProteins();
int proteinOverlap_inProgramExcluded_vs_NoExclusion = compareProteins(inProgramExcludedProteins, proteinsIdentifiedByNoExclusion);
ChangeValue(Header.NumProteinOverlap_ExcludedProteinsAgainstNoExclusionProteins, proteinOverlap_inProgramExcluded_vs_NoExclusion);
ChangeValue(Header.ProteinGroupsIdentified, ppr.getFilteredProteinGroups().Count);
}
