Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

FlashLFQ now enables selection of peptides to be quantified #779

Merged
merged 7 commits into from
May 20, 2024

Conversation

Alexander-Sol
Copy link
Contributor

Previously, flashLFQ would report peptides that had a peptide level FDR greater than 1%. This is because FlashLFQ usually only considers PSMs for input (and therefore, PSM level FDR).

I added an optional parameter to FlashLFQ engine that allows users to pass in a list of strings representing that peptides that should be quantified/used to calculate protein intensities.

Copy link

codecov bot commented May 6, 2024

Codecov Report

All modified and coverable lines are covered by tests ✅

Project coverage is 75.25%. Comparing base (a93dfa2) to head (b1f0d0a).

Additional details and impacted files

Impacted file tree graph

@@            Coverage Diff             @@
##           master     #779      +/-   ##
==========================================
+ Coverage   75.23%   75.25%   +0.02%     
==========================================
  Files         195      195              
  Lines       30455    30484      +29     
  Branches     3053     3058       +5     
==========================================
+ Hits        22913    22942      +29     
  Misses       7027     7027              
  Partials      515      515              
Files Coverage Δ
mzLib/FlashLFQ/FlashLFQResults.cs 91.84% <100.00%> (+0.08%) ⬆️
mzLib/FlashLFQ/FlashLfqEngine.cs 88.21% <100.00%> (+0.28%) ⬆️


public FlashLfqResults(List<SpectraFileInfo> spectraFiles, List<Identification> identifications)
public FlashLfqResults(List<SpectraFileInfo> spectraFiles, List<Identification> identifications, HashSet<string> peptides = null)
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

should peptides be peptideFullSequencesToQuantify?


foreach (SpectraFileInfo file in spectraFiles)
{
Peaks.Add(file, new List<ChromatographicPeak>());
}

foreach (Identification id in identifications)

// Only quantify peptides within the set of valid peptide sequences. This is done to enable pepitde-level FDR control of reported results
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

"full" sequences

/// Other peptides may appear in the QuantifiedPeaks output, but this list is used to enable
/// peptide-level FDR filtering
/// </summary>
public HashSet<string> PeptidesToQuantify { get; init; }
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

PeptideFullSequences

@@ -479,19 +480,29 @@ public static void TestFlashLfqMergeResults()
public static void TestFlashLfqMatchBetweenRuns()
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

can you add a comment for this test

nbollis
nbollis previously approved these changes May 20, 2024
Copy link
Contributor

@nbollis nbollis left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I agree with Shortreeds requested changes to the naming convention

@nbollis nbollis merged commit 84b2e9d into smith-chem-wisc:master May 20, 2024
3 checks passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

4 participants