I sent the link to your email. Thank you very much! On Wednesday, October 5, 2016 at 9:21:04 PM UTC-4, David Shteynberg wrote: > > I cannot tell from the search parameters if there is a problem. The only > issue I saw was that you didn't enable isotopic offsets, which are > sometimes helpful. If you send me a sample mzML file and possibly the > search results in pep.xml along with the database I can troubleshoot > further. One question, you said you are selecting 2+ 3+ and 4+, why are > only 2+ and 3+ in this dataset? Does any other dataset in your analysis > contain 4+? > > -David > > On Wed, Oct 5, 2016 at 6:16 PM, Ali <sma.ban...@gmail.com <javascript:>> > wrote: > >> Hi David >> >> We are only collecting 2+, 3+ and 4+ data. I am working on the cluster, >> so I am using the GPM version and we were unable to install the k-score >> plugin therefore I am using the X! Tandem native scoring for the searches. >> >> If you mean combining all of the search results with xinteract and then >> using PeptideProphet (as you instructed), I have not yet done so because I >> wanted to make sure that we are getting results from the individual >> searches. Is there an issue in my search parameters which is responsible >> for the problem? >> >> Many Thanks, >> Ali >> >> On Wednesday, October 5, 2016 at 7:06:19 PM UTC-4, David Shteynberg wrote: >>> >>> Hi Ali, >>> >>> >>> Are you setting your instrument to collect only 2+ and 3+ data? Are you >>> using the tandem bundled with TPP or the GPM version of tandem? Have you >>> tried to process all your runs in one file? >>> >>> Thanks, >>> -David >>> >>> >>> On Wed, Oct 5, 2016 at 11:39 AM, Ali <sma.ban...@gmail.com> wrote: >>> >>>> Thank you very much! I really appreciate it! >>>> >>>> On Tuesday, October 4, 2016 at 8:03:27 PM UTC-4, David Shteynberg wrote: >>>>> >>>>> PeptideProphet is doing it's job here by telling you that there is no >>>>> distinct population of correct results in this search. I suspect your >>>>> search parameters are the problem. This is orbitrap data so the charges >>>>> should be known, but it seems Tandem only searched 2+ and 3+ for you, >>>>> this >>>>> is somewhat troubling. Tomorrow, I will check your parameters more >>>>> carefully to try to identify the actual problem with your search. >>>>> >>>>> -David >>>>> >>>>> >>>>> On Tue, Oct 4, 2016 at 2:12 PM, Ali <sma.ban...@gmail.com> wrote: >>>>> >>>>>> It appears that PeptideProphet fails, I am getting this message: >>>>>> >>>>>> using Accurate Mass Bins >>>>>> using PPM mass difference >>>>>> (X! Tandem) (minprob 0) >>>>>> adding Accurate Mass mixture distr >>>>>> init with X! Tandem trypsin >>>>>> MS Instrument info: Manufacturer: Thermo Scientific, Model: UNKNOWN, >>>>>> Ionization: nanoelectrospray, Analyzer: orbitrap, Detector: inductive >>>>>> detector >>>>>> >>>>>> PeptideProphet (TPP v4.8.0 PHILAE, Build 201608051257-exported >>>>>> (RedHatEnterpriseServer-x86_64)) AKeller@ISB >>>>>> read in 0 1+, 14418 2+, 11080 3+, 0 4+, 0 5+, 0 6+, and 0 7+ spectra. >>>>>> Initialising statistical models ... >>>>>> negmean = 0.0533258 >>>>>> negmean = 1 >>>>>> negmean = 0.0376944 >>>>>> negmean = -0.2 >>>>>> negmean = -0.2 >>>>>> negmean = -0.2 >>>>>> negmean = -0.2 >>>>>> Iterations: .........10.........20. >>>>>> WARNING: Mixture model quality test failed for charge (2+). >>>>>> WARNING: Mixture model quality test failed for charge (3+). >>>>>> >>>>>> On Monday, October 3, 2016 at 5:57:39 PM UTC-4, David Shteynberg >>>>>> wrote: >>>>>>> >>>>>>> I highly recommend you start employing decoys in your search >>>>>>> databases. I noticed you are not combining together the runs in the >>>>>>> analysis. If the data was generated and searched in the same way it >>>>>>> should >>>>>>> be analysed together for better statistical power. I suggest you >>>>>>> combine >>>>>>> the results from all your runs and filter by the spectrum names later. >>>>>>> Try >>>>>>> this: >>>>>>> >>>>>>> InteractParser ./ALL.interact.xml ./*.tandem.pep.xml -S -Etrypsin >>>>>>> PeptideProphetParser ./ALL.interact.xml ACCMASS PPM NOICAT MINPROB=0 >>>>>>> >>>>>>> >>>>>>> -David >>>>>>> >>>>>>> On Mon, Oct 3, 2016 at 2:39 PM, Ali <sma.ban...@gmail.com> wrote: >>>>>>> >>>>>>>> Hi David >>>>>>>> >>>>>>>> Thanks for your reply. I don't think so, it is not just one set, as >>>>>>>> I said almost a third of my search results are returned with zero >>>>>>>> probabilities. All of the searches are done with the same setting and >>>>>>>> using >>>>>>>> the same database. The fact that for some of the searched files I am >>>>>>>> getting zero probabilities and some of them are just fine is a bit >>>>>>>> strange. >>>>>>>> I don't have decoys in my data base. If it helps, I can send you one >>>>>>>> of the >>>>>>>> search results with or without PeptideProphet probabilities. >>>>>>>> >>>>>>>> Sincerely, >>>>>>>> Ali >>>>>>>> >>>>>>>> >>>>>>>> >>>>>>>> On Monday, October 3, 2016 at 4:19:05 PM UTC-4, David Shteynberg >>>>>>>> wrote: >>>>>>>>> >>>>>>>>> Is it possible there are no correct results in that set? Common >>>>>>>>> pitfalls are incorrect search parameters (e.g. wrong mods enabled or >>>>>>>>> disabled) and wrong database. If you have decoys in your database >>>>>>>>> these >>>>>>>>> can be utilized for a more accurate mixture model generation using >>>>>>>>> options >>>>>>>>> "DECOY=<your decoy prefix> DECOYPROBS NONPARAM" in additions to the >>>>>>>>> options you already have enabled. >>>>>>>>> >>>>>>>>> -David >>>>>>>>> >>>>>>>>> On Mon, Oct 3, 2016 at 10:39 AM, Ali <sma.ban...@gmail.com> wrote: >>>>>>>>> >>>>>>>>>> Hi >>>>>>>>>> >>>>>>>>>> I am using PeptideProphet on my X!Tandem search results. There >>>>>>>>>> are ~150 experiments and PeptideProphet works fine on most of them >>>>>>>>>> but for >>>>>>>>>> ~50 experiments, PeptideProphet returns zero probability for all of >>>>>>>>>> the >>>>>>>>>> peptides. Can some one please tell me what might be the issue? >>>>>>>>>> >>>>>>>>>> *More Info:* >>>>>>>>>> I am using TPP v4.8 on Linux. PeptideProphet does not show an >>>>>>>>>> error when it produces zero probability for all of the search >>>>>>>>>> results. >>>>>>>>>> This is what I'm doing: >>>>>>>>>> -search with X!Tandem >>>>>>>>>> -Tandem2XML ./100.tandem.xml ./100.tandem.pep.xml >>>>>>>>>> -InteractParser ./100.interact.xml ./100.tandem.pep.xml -S >>>>>>>>>> -Etrypsin >>>>>>>>>> -PeptideProphetParser ./100.interact.xml ACCMASS PPM NOICAT >>>>>>>>>> MINPROB=0 >>>>>>>>>> >>>>>>>>>> >>>>>>>>>> Thank you very much, >>>>>>>>>> *Ali* >>>>>>>>>> >>>>>>>>>> -- >>>>>>>>>> You received this message because you are subscribed to the >>>>>>>>>> Google Groups "spctools-discuss" group. >>>>>>>>>> To unsubscribe from this group and stop receiving emails from it, >>>>>>>>>> send an email to spctools-discu...@googlegroups.com. >>>>>>>>>> To post to this group, send email to spctools...@googlegroups.com >>>>>>>>>> . >>>>>>>>>> Visit this group at >>>>>>>>>> https://groups.google.com/group/spctools-discuss. >>>>>>>>>> For more options, visit https://groups.google.com/d/optout. >>>>>>>>>> >>>>>>>>> >>>>>>>>> -- >>>>>>>> You received this message because you are subscribed to the Google >>>>>>>> Groups "spctools-discuss" group. >>>>>>>> To unsubscribe from this group and stop receiving emails from it, >>>>>>>> send an email to spctools-discu...@googlegroups.com. >>>>>>>> To post to this group, send email to spctools...@googlegroups.com. >>>>>>>> Visit this group at >>>>>>>> https://groups.google.com/group/spctools-discuss. >>>>>>>> For more options, visit https://groups.google.com/d/optout. >>>>>>>> >>>>>>> >>>>>>> -- >>>>>> You received this message because you are subscribed to the Google >>>>>> Groups "spctools-discuss" group. >>>>>> To unsubscribe from this group and stop receiving emails from it, >>>>>> send an email to spctools-discu...@googlegroups.com. >>>>>> To post to this group, send email to spctools...@googlegroups.com. >>>>>> Visit this group at https://groups.google.com/group/spctools-discuss. >>>>>> For more options, visit https://groups.google.com/d/optout. >>>>>> >>>>> >>>>> -- >>>> You received this message because you are subscribed to the Google >>>> Groups "spctools-discuss" group. >>>> To unsubscribe from this group and stop receiving emails from it, send >>>> an email to spctools-discu...@googlegroups.com. >>>> To post to this group, send email to spctools...@googlegroups.com. >>>> Visit this group at https://groups.google.com/group/spctools-discuss. >>>> For more options, visit https://groups.google.com/d/optout. >>>> >>> >>> -- >> You received this message because you are subscribed to the Google Groups >> "spctools-discuss" group. >> To unsubscribe from this group and stop receiving emails from it, send an >> email to spctools-discu...@googlegroups.com <javascript:>. >> To post to this group, send email to spctools...@googlegroups.com >> <javascript:>. >> Visit this group at https://groups.google.com/group/spctools-discuss. >> For more options, visit https://groups.google.com/d/optout. >> > >
-- You received this message because you are subscribed to the Google Groups "spctools-discuss" group. To unsubscribe from this group and stop receiving emails from it, send an email to spctools-discuss+unsubscr...@googlegroups.com. To post to this group, send email to spctools-discuss@googlegroups.com. Visit this group at https://groups.google.com/group/spctools-discuss. For more options, visit https://groups.google.com/d/optout.