my $file= "a.mgf"; my $format="mgf"; my $sp=InSilicoSpectro::Spectra::MSSpectra->new(source=>$file, format=>$format); $sp->open();
my $sf = new InSilicoSpectro::Spectra::Filter::MSFilter(); $sf->readXml('a.xml'); $sf->filterSpectra($sp);
This class allows for quality assasement and filtering of spectra or peaks using direct (or very simple) accesible informations.
create a new object.
opens the provided XML-file containing the information which filter to use and its parameters.
reads the TwigEl passed by $el and executes twig_addSpectrumFilter() with this values.
apply the XML-file previously loaded on the Spectra (can be MSRun, MS, MSMS, or MSCmpd-Spectra).
applies the action on the currentSpectra under the threshold conditions.
only prepared for further classes inheriting from this one.
checks if the different values provided by the xml are valid.
stores the at the moment just considers the size of fragments
sums up the fragment intensities or moz-values and stores the value for each compound in filterValues. The method summ() corresponds to the filter sum. It is renamed, because List::Utils allready uses this function-name.
the natural log of the filterValue. This method is renamed to log_() because this name is already used.
gets the values (moz, intensity) from the fragments of the currentSpectra (has to be an MSCmpd). The results are stored into filterValue. Which values should be extracted can be selected by $val, otherwise the next name in the filterName is used.
gets the values (moz, intensity) from the precursors of the currentSpectra (has to be an MSMSSpectra). The The results are stored into filterValue. Which values should be extracted can be selected by $val, otherwise the next name in the filterName is used.
calls fragment("intensity"). Allows you just to write intensity as a filter-name instead.
calls fragment("moz"). Allows you just to write moz as a filter-name instead.
calculates the rank from the intensity of the fragments from the smallest first and biggest last (1, 3, 3, 4, 5, ..) and divides it by the total number of peaks.
makes the normalisation of the int-peaks as described in the paper from M. Bern et al. Because we can make a selection of peaks to be tested before, it doesnt make a lot of sense in this context.
set or get a current spectra.
set or get the threshold value.
set or get the filter value.
see MSFilterDirecValue.t for an example.
Copyright (C) 2004-2005 Geneva Bioinformatics www.genebio.com
This library is free software; you can redistribute it and/or modify it under the terms of the GNU Lesser General Public License as published by the Free Software Foundation; either version 2.1 of the License, or (at your option) any later version.
This library is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for more details.
You should have received a copy of the GNU Lesser General Public License along with this library; if not, write to the Free Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
Roman Mylonas, www.genebio.com