public PeakBoundsMatch(TransitionGroupChromInfo chromInfoTrue, TransitionGroupChromInfo chromInfoPicked, MatchKey key, ChromatogramSet chromSet, MsDataFileUri filePath, TransitionGroupDocNode nodeGroup, PeptideDocNode nodePep, bool hasNoQValues, bool hasNoScores, bool apexPresent) { _chromatogramSet = chromSet; ChromInfoTrue = chromInfoTrue; ChromInfoPicked = chromInfoPicked; Key = key; FilePath = filePath; NodeGroup = nodeGroup; NodePep = nodePep; // Read apex if (apexPresent) { PickedApex = ParsePickedTimeAnnotation(ComparePeakBoundaries.APEX_ANNOTATION); PickedStartBoundary = ParsePickedTimeAnnotation(ComparePeakBoundaries.START_TIME_ANNOTATION); PickedEndBoundary = ParsePickedTimeAnnotation(ComparePeakBoundaries.END_TIME_ANNOTATION); } else { PickedApex = ChromInfoPicked.RetentionTime; PickedStartBoundary = ChromInfoPicked.StartRetentionTime; PickedEndBoundary = ChromInfoPicked.EndRetentionTime; } QValue = GetScoreValue(ChromInfoPicked, MProphetResultsHandler.AnnotationName, ci => ci.QValue); if (QValue == null && !IsMissingPickedPeak && !hasNoQValues) { throw new IOException(string.Format( Resources.PeakBoundsMatch_QValue_Unable_to_read_q_value_annotation_for_peptide__0__of_file__1_, ModifiedSequence, FileName)); } Score = GetScoreValue(ChromInfoPicked, MProphetResultsHandler.MAnnotationName, ci => ci.ZScore); if (Score == null && !IsMissingPickedPeak && !hasNoScores) { throw new IOException(string.Format( Resources.PeakBoundsMatch_PeakBoundsMatch_Unable_to_read_a_score_annotation_for_peptide__0__of_file__1_, ModifiedSequence, FileName)); } }
protected bool Equals(MatchKey other) { return(FileId == other.FileId && TransitionGroupId == other.TransitionGroupId); }
public PeakBoundsMatch(TransitionGroupChromInfo chromInfoTrue, TransitionGroupChromInfo chromInfoPicked, MatchKey key, MsDataFileUri filePath, TransitionGroupDocNode nodeGroup, bool hasNoQValues, bool hasNoScores, bool apexPresent) { ChromInfoTrue = chromInfoTrue; ChromInfoPicked = chromInfoPicked; Key = key; FilePath = filePath; NodeGroup = nodeGroup; // Read apex if (apexPresent) { string apexString = ChromInfoPicked.Annotations.GetAnnotation(ComparePeakBoundaries.APEX_ANNOTATION); double apexDouble; if (!double.TryParse(apexString, out apexDouble)) { if (apexString == null || apexString.Equals(TextUtil.EXCEL_NA)) { PickedApex = null; } else { throw new IOException(string.Format(Resources.PeakBoundsMatch_PeakBoundsMatch_Unable_to_read_apex_retention_time_value_for_peptide__0__of_file__1__, Sequence, FileName)); } } else { PickedApex = apexDouble; } } else { PickedApex = ChromInfoPicked.RetentionTime; } string qValueString = ChromInfoPicked.Annotations.GetAnnotation(MProphetResultsHandler.AnnotationName); double qValueDouble; // qValue MUST be present unless the peak is null or ALL q Values are absent if (!double.TryParse(qValueString, out qValueDouble)) { if (!IsMissingPickedPeak && !hasNoQValues) { throw new IOException(string.Format(Resources.PeakBoundsMatch_QValue_Unable_to_read_q_value_annotation_for_peptide__0__of_file__1_, Sequence, FileName)); } QValue = null; } else { QValue = qValueDouble; } // Same for score string scoreString = ChromInfoPicked.Annotations.GetAnnotation(MProphetResultsHandler.MAnnotationName); double scoreDouble; if (!double.TryParse(scoreString, out scoreDouble)) { if (!IsMissingPickedPeak && !hasNoScores) { throw new IOException(string.Format(Resources.PeakBoundsMatch_QValue_Unable_to_read_q_value_annotation_for_peptide__0__of_file__1_, Sequence, FileName)); } Score = null; } else { Score = scoreDouble; } }
public void GenerateComparison(SrmDocument docOriginal, IProgressMonitor progressMonitor) { DocOriginal = docOriginal; _docCompare = docOriginal; if (IsModel) { var handler = new MProphetResultsHandler(DocOriginal, PeakScoringModel) { AddAnnotation = true, AddMAnnotation = true }; handler.ScoreFeatures(progressMonitor); _docCompare = handler.ChangePeaks(progressMonitor); } else { // Add in the necessary annotation definitions so that ImportPeakBoundaries can read them _docCompare = AddAnnotationIfMissing(_docCompare, MProphetResultsHandler.AnnotationName); _docCompare = AddAnnotationIfMissing(_docCompare, MProphetResultsHandler.MAnnotationName); _docCompare = AddAnnotationIfMissing(_docCompare, APEX_ANNOTATION); // But strip the values of these annotations if there were any var annotationNamesToStrip = new List <string> { MProphetResultsHandler.AnnotationName, MProphetResultsHandler.MAnnotationName, APEX_ANNOTATION }; var annotationNamesToKeep = _docCompare.Settings.DataSettings.AnnotationDefs .Where(def => !annotationNamesToStrip.Contains(def.Name)) .Select(def => def.Name).ToList(); _docCompare = (SrmDocument)_docCompare.StripAnnotationValues(annotationNamesToKeep); Importer = new PeakBoundaryImporter(_docCompare); long lineCount = Helpers.CountLinesInFile(FilePath); // Peek at the file to see if it has a column called "Apex" using (var reader = new StreamReader(FilePath)) { var line = reader.ReadLine(); var fieldNames = PeakBoundaryImporter.FIELD_NAMES.ToList(); fieldNames.Add(new [] { APEX_ANNOTATION }); int fieldsTotal; int[] fieldIndices; var separator = PeakBoundaryImporter.DetermineCorrectSeparator(line, fieldNames, PeakBoundaryImporter.REQUIRED_NO_CHROM, out fieldIndices, out fieldsTotal); ApexPresent = separator != null && fieldIndices.Last() != -1; } _docCompare = Importer.Import(FilePath, progressMonitor, lineCount, true, !ApexPresent); } var matches = new List <PeakBoundsMatch>(); var truePeptides = DocOriginal.Molecules.ToList(); var pickedPeptides = _docCompare.Molecules.ToList(); int nTruePeptides = truePeptides.Count; var chromatogramSets = DocOriginal.Settings.MeasuredResults.Chromatograms; for (int i = 0; i < chromatogramSets.Count; i++) { var set = chromatogramSets[i]; for (int j = 0; j < nTruePeptides; j++) { var truePeptide = truePeptides[j]; var pickedPeptide = pickedPeptides[j]; // We don't care about peak picking performance for decoys if (truePeptide.IsDecoy) { continue; } // We don't care about peak picking performance for standard peptides if (truePeptide.GlobalStandardType != null) { continue; } var trueGroups = truePeptide.TransitionGroups.ToList(); var pickedGroups = pickedPeptide.TransitionGroups.ToList(); for (int k = 0; k < trueGroups.Count; k++) { var trueGroup = trueGroups[k]; var pickedGroup = pickedGroups[k]; var trueResults = trueGroup.Results[i]; var pickedResults = pickedGroup.Results[i]; if (trueResults == null) { continue; } int nChromInfos = trueResults.Count; for (int m = 0; m < nChromInfos; m++) { var trueInfo = trueResults[m]; var pickedInfo = pickedResults[m]; var fileInfo = set.GetFileInfo(trueInfo.FileId); var filePath = fileInfo.FilePath; var key = new MatchKey(trueGroup.Id.GlobalIndex, trueInfo.FileId.GlobalIndex); matches.Add(new PeakBoundsMatch(trueInfo, pickedInfo, key, filePath, trueGroup, HasNoQValues, HasNoScores, ApexPresent)); } } } } Matches = matches; // Detect if the apex time is in seconds, and if so adjust it to minutes if (ApexPresent) { double maxTime = matches.Where(match => match.PickedApex.HasValue) .Select(match => match.PickedApex.Value).Max(); if (maxTime > PeakBoundaryImporter.GetMaxRt(DocOriginal)) { foreach (var match in matches.Where(match => match.PickedApex.HasValue)) { match.PickedApex = match.PickedApex / 60; } } } }
public void GenerateComparison(SrmDocument docOriginal, IProgressMonitor progressMonitor) { DocOriginal = docOriginal; _docCompare = docOriginal; if (IsModel) { var handler = new MProphetResultsHandler(DocOriginal, PeakScoringModel); handler.ScoreFeatures(progressMonitor, true); _docCompare = handler.ChangePeaks(progressMonitor); } else { // Add in the necessary annotation definitions so that ImportPeakBoundaries can read them var annotationNamesToStrip = new List <string> { MProphetResultsHandler.AnnotationName, MProphetResultsHandler.MAnnotationName, APEX_ANNOTATION, START_TIME_ANNOTATION, END_TIME_ANNOTATION }; foreach (var annotationName in annotationNamesToStrip) { _docCompare = AddAnnotationIfMissing(_docCompare, annotationName); } // But strip the values of these annotations if there were any var annotationNamesToKeep = _docCompare.Settings.DataSettings.AnnotationDefs .Where(def => !annotationNamesToStrip.Contains(def.Name)) .Select(def => def.Name).ToList(); _docCompare = (SrmDocument)_docCompare.StripAnnotationValues(annotationNamesToKeep); Importer = new PeakBoundaryImporter(_docCompare); long lineCount = Helpers.CountLinesInFile(FilePath); // Peek at the file to see if it has a column called "Apex" using (var reader = new StreamReader(FilePath)) { var line = reader.ReadLine(); var fieldNames = PeakBoundaryImporter.FIELD_NAMES.ToList(); int fieldsTotal; int[] fieldIndices; var separator = PeakBoundaryImporter.DetermineCorrectSeparator(line, fieldNames, PeakBoundaryImporter.REQUIRED_NO_CHROM, out fieldIndices, out fieldsTotal); ApexPresent = separator != null && fieldIndices[(int)PeakBoundaryImporter.Field.apex_time] != -1; } _docCompare = Importer.Import(FilePath, progressMonitor, lineCount, true, !ApexPresent); if (progressMonitor.IsCanceled) { return; } } var matches = new List <PeakBoundsMatch>(); var truePeptides = DocOriginal.Molecules.ToList(); var pickedPeptides = _docCompare.Molecules.ToList(); int nTruePeptides = truePeptides.Count; var chromatogramSets = DocOriginal.Settings.MeasuredResults.Chromatograms; for (int i = 0; i < chromatogramSets.Count; i++) { var chromSet = chromatogramSets[i]; var set = chromatogramSets[i]; for (int j = 0; j < nTruePeptides; j++) { var truePeptide = truePeptides[j]; var pickedPeptide = pickedPeptides[j]; // We don't care about peak picking performance for decoys if (truePeptide.IsDecoy) { continue; } // We don't care about peak picking performance for standard peptides if (truePeptide.GlobalStandardType != null) { continue; } var trueGroups = truePeptide.TransitionGroups.ToList(); var pickedGroups = pickedPeptide.TransitionGroups.ToList(); for (int k = 0; k < trueGroups.Count; k++) { var trueGroup = trueGroups[k]; var pickedGroup = pickedGroups[k]; var trueResults = trueGroup.Results[i]; var pickedResults = pickedGroup.Results[i]; if (trueResults.IsEmpty) { continue; } int nChromInfos = trueResults.Count; for (int m = 0; m < nChromInfos; m++) { var trueInfo = trueResults[m]; var pickedInfo = pickedResults[m]; // For TRIC testing at the moment // double? qvaluePicked = PeakBoundsMatch.GetScoreValue(pickedInfo, // MProphetResultsHandler.AnnotationName, ci => pickedInfo.QValue); // if (qvaluePicked.HasValue && qvaluePicked.Value > 0.01) // continue; var fileInfo = set.GetFileInfo(trueInfo.FileId); var filePath = fileInfo.FilePath; var key = new MatchKey(trueGroup.Id.GlobalIndex, trueInfo.FileId.GlobalIndex); matches.Add(new PeakBoundsMatch(trueInfo, pickedInfo, key, chromSet, filePath, trueGroup, truePeptide, HasNoQValues, HasNoScores, ApexPresent)); } } } } Matches = matches; // Detect if the apex time is in seconds, and if so adjust it to minutes if (ApexPresent) { // ReSharper disable PossibleMultipleEnumeration var matchesWithApex = matches.Where(match => match.PickedApex.HasValue).ToArray(); double maxTime = 0; if (matchesWithApex.Length > 0) { maxTime = matchesWithApex.Select(match => match.PickedApex.Value).Max(); } if (maxTime > PeakBoundaryImporter.GetMaxRt(DocOriginal)) { foreach (var match in matchesWithApex) { // If the end is greater than the Apex, then it must also be in seconds (otherwise minutes) if (match.PickedEndBoundary.HasValue && match.PickedEndBoundary > match.PickedApex) { match.PickedEndBoundary = match.PickedEndBoundary / 60; } match.PickedApex = match.PickedApex / 60; // If the start is now greater than the apex, then it must also be in seconds (otherwise minutes) if (match.PickedStartBoundary.HasValue && match.PickedStartBoundary > match.PickedApex) { match.PickedStartBoundary = match.PickedStartBoundary / 60; } } } // ReSharper restore PossibleMultipleEnumeration } }
public PeakBoundsMatch(TransitionGroupChromInfo chromInfoTrue, TransitionGroupChromInfo chromInfoPicked, MatchKey key, MsDataFileUri filePath, TransitionGroupDocNode nodeGroup, bool hasNoQValues, bool hasNoScores, bool apexPresent) { ChromInfoTrue = chromInfoTrue; ChromInfoPicked = chromInfoPicked; Key = key; FilePath = filePath; NodeGroup = nodeGroup; // Read apex if (apexPresent) { string apexString = ChromInfoPicked.Annotations.GetAnnotation(ComparePeakBoundaries.APEX_ANNOTATION); double apexDouble; if (!double.TryParse(apexString, out apexDouble)) { if (apexString == null || apexString.Equals(TextUtil.EXCEL_NA)) { PickedApex = null; } else { throw new IOException(string.Format(Resources.PeakBoundsMatch_PeakBoundsMatch_Unable_to_read_apex_retention_time_value_for_peptide__0__of_file__1__, Sequence, FileName)); } } else { PickedApex = apexDouble; } } else { PickedApex = ChromInfoPicked.RetentionTime; } string qValueString = ChromInfoPicked.Annotations.GetAnnotation(MProphetResultsHandler.AnnotationName); double qValueDouble; // qValue MUST be present unless the peak is null or ALL q Values are absent if (!double.TryParse(qValueString, out qValueDouble)) { if (!IsMissingPickedPeak && !hasNoQValues) { throw new IOException(string.Format(Resources.PeakBoundsMatch_QValue_Unable_to_read_q_value_annotation_for_peptide__0__of_file__1_, Sequence, FileName)); } QValue = null; } else { QValue = qValueDouble; } // Same for score string scoreString = ChromInfoPicked.Annotations.GetAnnotation(MProphetResultsHandler.MAnnotationName); double scoreDouble; if (!double.TryParse(scoreString, out scoreDouble)) { if (!IsMissingPickedPeak && !hasNoScores) { throw new IOException(string.Format(Resources.PeakBoundsMatch_QValue_Unable_to_read_q_value_annotation_for_peptide__0__of_file__1_, Sequence, FileName)); } Score = null; } else { Score = scoreDouble; } }
protected bool Equals(MatchKey other) { return FileId == other.FileId && TransitionGroupId == other.TransitionGroupId; }
public void GenerateComparison(SrmDocument docOriginal, IProgressMonitor progressMonitor) { DocOriginal = docOriginal; _docCompare = docOriginal; if (IsModel) { var handler = new MProphetResultsHandler(DocOriginal, PeakScoringModel) { AddAnnotation = true, AddMAnnotation = true }; handler.ScoreFeatures(progressMonitor); _docCompare = handler.ChangePeaks(progressMonitor); } else { // Add in the necessary annotation definitions so that ImportPeakBoundaries can read them _docCompare = AddAnnotationIfMissing(_docCompare, MProphetResultsHandler.AnnotationName); _docCompare = AddAnnotationIfMissing(_docCompare, MProphetResultsHandler.MAnnotationName); _docCompare = AddAnnotationIfMissing(_docCompare, APEX_ANNOTATION); // But strip the values of these annotations if there were any var annotationNamesToStrip = new List<string> { MProphetResultsHandler.AnnotationName, MProphetResultsHandler.MAnnotationName, APEX_ANNOTATION }; var annotationNamesToKeep = _docCompare.Settings.DataSettings.AnnotationDefs .Where(def => !annotationNamesToStrip.Contains(def.Name)) .Select(def => def.Name).ToList(); _docCompare = (SrmDocument)_docCompare.StripAnnotationValues(annotationNamesToKeep); Importer = new PeakBoundaryImporter(_docCompare); long lineCount = Helpers.CountLinesInFile(FilePath); // Peek at the file to see if it has a column called "Apex" using (var reader = new StreamReader(FilePath)) { var line = reader.ReadLine(); var fieldNames = PeakBoundaryImporter.FIELD_NAMES.ToList(); fieldNames.Add(new [] {APEX_ANNOTATION}); int fieldsTotal; int[] fieldIndices; var separator = PeakBoundaryImporter.DetermineCorrectSeparator(line, fieldNames, PeakBoundaryImporter.REQUIRED_NO_CHROM, out fieldIndices, out fieldsTotal); ApexPresent = separator != null && fieldIndices.Last() != -1; } _docCompare = Importer.Import(FilePath, progressMonitor, lineCount, true, !ApexPresent); } var matches = new List<PeakBoundsMatch>(); var truePeptides = DocOriginal.Molecules.ToList(); var pickedPeptides = _docCompare.Molecules.ToList(); int nTruePeptides = truePeptides.Count; var chromatogramSets = DocOriginal.Settings.MeasuredResults.Chromatograms; for (int i = 0; i < chromatogramSets.Count; i++) { var set = chromatogramSets[i]; for (int j = 0; j < nTruePeptides; j++) { var truePeptide = truePeptides[j]; var pickedPeptide = pickedPeptides[j]; // We don't care about peak picking performance for decoys if (truePeptide.IsDecoy) continue; // We don't care about peak picking performance for standard peptides if (truePeptide.GlobalStandardType != null) continue; var trueGroups = truePeptide.TransitionGroups.ToList(); var pickedGroups = pickedPeptide.TransitionGroups.ToList(); for (int k = 0; k < trueGroups.Count; k++) { var trueGroup = trueGroups[k]; var pickedGroup = pickedGroups[k]; var trueResults = trueGroup.Results[i]; var pickedResults = pickedGroup.Results[i]; if (trueResults == null) continue; int nChromInfos = trueResults.Count; for (int m = 0; m < nChromInfos; m++) { var trueInfo = trueResults[m]; var pickedInfo = pickedResults[m]; var fileInfo = set.GetFileInfo(trueInfo.FileId); var filePath = fileInfo.FilePath; var key = new MatchKey(trueGroup.Id.GlobalIndex, trueInfo.FileId.GlobalIndex); matches.Add(new PeakBoundsMatch(trueInfo, pickedInfo, key, filePath, trueGroup, HasNoQValues, HasNoScores, ApexPresent)); } } } } Matches = matches; // Detect if the apex time is in seconds, and if so adjust it to minutes if (ApexPresent) { double maxTime = matches.Where(match => match.PickedApex.HasValue) .Select(match => match.PickedApex.Value).Max(); if (maxTime > PeakBoundaryImporter.GetMaxRt(DocOriginal)) { foreach (var match in matches.Where(match => match.PickedApex.HasValue)) { match.PickedApex = match.PickedApex / 60; } } } }