These common commands can be used on the default data set obtained when downloading microPITA, simply cut and paste them into a commandline in the downloaded microPITA directory.
I. PCL file or BIOM file
BIOM file definition: For BIOM file definition please see http://biom-format.org/
PCL file definition: Although some defaults can be changed, microPITA expects a PCL file as an input file. Several PCL files are supplied by default in the input directory. A PCL file is a TEXT delimited file similar to an excel spread sheet with the following characteristics.
- Rows represent metadata and features (bugs), columns represent samples.
- The first row by default should be the sample ids.
- Metadata rows should be next.
- Lastly, rows containing features (bugs) measurements (like abundance) should be after metadata rows.
- The first column should contain the ID describing the column. For metadata this may be, for example, "Age" for a row containing the age of the patients donating the samples. For measurements, this should be the feature name (bug name).
- By default the file is expected to be TAB delimited.
- If a consensus lineage or hierarchy of taxonomy is contained in the feature name, the default delimiter between clades is the pipe ("|").
II. Targeted feature file If using the targeted feature methodology, you will need to provide a txt file listing the feature(s) of interest. Each feature should be on it's own line and should be written as found in the input PCL file.
Please note, all calls to microPITA should work interchangeably with PCL or BIOM files. BIOM files do not require the --lastmeta or --id arguments.
There are four unsupervised methods which can be performed: diverse (maximum diversity), extreme (most dissimilar), representative (representative dissimilarity) and features (targeted feature).
The first three methods are performed as follows (selecting a default 10 samples):
python MicroPITA.py --lastmeta Label -m representative input/Test.pcl output.txt python MicroPITA.py -m representative input/Test.biom output.txt
python MicroPITA.py --lastmeta Label -m diverse input/Test.pcl output.txt python MicroPITA.py -m diverse input/Test.biom output.txt
python MicroPITA.py --lastmeta Label -m extreme input/Test.pcl output.txt python MicroPITA.py -m extreme input/Test.biom output.txt
Each of the previous methods are made up of the following pieces:
- python MicroPITA.py to call the microPITA script.
- --lastmeta which indicates the keyword (first column value) of the last row that contains metadata (PCL input only).
- -m which indicates the method to use in selection.
- input/Test.pcl or input/Test.biom which is the first positional argument indicating an input file
- output.txt which is the second positional argument indicating the location to write to the output file.
Selecting specific features has additional arguments to consider --targets (required) and --feature_method (optional).
python MicroPITA.py --lastmeta Label -m features --targets input/TestFeatures.taxa input/Test.pcl output.txt python MicroPITA.py -m features --targets input/TestFeatures.taxa input/Test.biom output.txt
python MicroPITA.py --lastmeta Label -m features --feature_method abundance --targets input/TestFeatures.taxa input/Test.pcl output.txt python MicroPITA.py -m features --feature_method abundance --targets input/TestFeatures.taxa input/Test.biom output.txt
These additional arguments are described as:
- --targets The path to the file that has the features (bugs or clades) of interest. Make sure they are written as they appear in your input file!
- --feature_method is the method of selection used and can be based on ranked abundance ("rank") or abundance ("abundance"). The default value is rank. To differentiate the methods, rank tends to select samples in which the feature dominates the samples regardless of it's abundance. Abundance tends to select samples in which the feature is most abundant without a guarantee that the feature is the most abundant feature in the sample.
Two supervised methods are also available: distinct and discriminant
These methods require an additional argument --label which is the first column keyword of the row used to classify samples for the supervised methods. These methods can be performed as follows:
python MicroPITA.py --lastmeta Label --label Label -m distinct input/Test.pcl output.txt python MicroPITA.py --label Label -m distinct input/Test.biom output.txt
python MicroPITA.py --lastmeta Label --label Label -m discriminant input/Test.pcl output.txt python MicroPITA.py --label Label -m discriminant input/Test.biom output.txt
The default alpha diversity for the maximum diversity sampling method is inverse simpson; the default beta-diversity for representative and most dissimilar selection is bray-curtis dissimilarity. There are several mechanisms that allow one to change this. You may:
- Choose from a selection of alpha-diveristy metrics. Note when supplying an alpha diversity. This will affect the maximum diveristy sampling method only. Please make sure to use a diversity metric where the larger number indicates a higher diversity. If this is not the case make sure to use the -f or --invertDiversity flag to invert the metric. The inversion is multiplicative (1/alpha-metric).
python MicroPITA.py --lastmeta Label -m diverse -a simpson input/Test.pcl output.txt python MicroPITA.py -m diverse -a simpson input/Test.biom output.txt
A case where inverting the metric is needed.
python MicroPITA.py --lastmeta Label -m diverse -a dominance -f input/Test.pcl output.txt python MicroPITA.py -m diverse -a dominance -f input/Test.biom output.txt
- Choose from a selection of beta-diversity metrics. Note when supplying a beta-diversity. This will effect both the representative and most dissimilar sampling methods. The metric as given will be used for the representative method while 1-beta-metric is used for the most dissimilar.
python MicroPITA.py --lastmeta Label -m representative -b euclidean input/Test.pcl output.txt python MicroPITA.py -m representative -b euclidean input/Test.biom output.txt
python MicroPITA.py --lastmeta Label -m extreme -b euclidean input/Test.pcl output.txt python MicroPITA.py -m extreme -b euclidean input/Test.biom output.txt
Note for using Unifrac. Both Weighted and Unweighted unifrac are available for use. Make sure to supply the associated tree (-o, --tree) and environment files (-i,--envr) as well as indicate using Unifrac with (-b,--beta)
python MicroPITA.py --lastmeta Label -m extreme -b unifrac_weighted -o input/Test.tree -i input/Test-env.txt input/Test.pcl output.txt python MicroPITA.py -m extreme -b unifrac_weighted -o input/Test.tree -i input/Test-env.txt input/Test.biom output.txt python MicroPITA.py --lastmeta Label -m extreme -b unifrac_unweighted -o input/Test.tree -i input/Test-env.txt input/Test.pcl output.txt python MicroPITA.py -m extreme -b unifrac_unweighted -o input/Test.tree -i input/Test-env.txt input/Test.biom output.txt python MicroPITA.py --lastmeta Label -m representative -b unifrac_weighted -o input/Test.tree -i input/Test-env.txt input/Test.pcl output.txt python MicroPITA.py -m representative -b unifrac_weighted -o input/Test.tree -i input/Test-env.txt input/Test.biom output.txt python MicroPITA.py --lastmeta Label -m representative -b unifrac_unweighted -o input/Test.tree -i input/Test-env.txt input/Test.pcl output.txt python MicroPITA.py -m representative -b unifrac_unweighted -o input/Test.tree -i input/Test-env.txt input/Test.biom output.txt
- Supply your own custom alpha-diversity per sample as a metadata (row) in your pcl file.
python MicroPITA.py --lastmeta Label -m diverse -q alpha_custom input/Test.pcl output.txt python MicroPITA.py -m diverse -q alpha_custom input/Test2.biom output.txt
- Supply your own custom beta diversity as a matrix (provided in a seperate file).
python MicroPITA.py --lastmeta Label -m representative -x input/Test_Matrix.txt input/Test.pcl output.txt python MicroPITA.py -m representative -x input/Test_Matrix.txt input/Test.biom output.txt python MicroPITA.py --lastmeta Label -m extreme -x input/Test_Matrix.txt input/Test.pcl output.txt python MicroPITA.py -m extreme -x input/Test_Matrix.txt input/Test.biom output.txt
Sample Selection: To change the number of selected samples for any method use the -n argument. This example selects 6 representative samples instead of the default 10.
python MicroPITA.py --lastmeta Label -m representative -n 6 input/Test.pcl output.txt python MicroPITA.py -m representative -n 6 input/Test.biom output.txt
When using a supervised method this indicates how many samples will be selected per class of sample. For example if you are performing supervised selection of 6 samples (-n 6) on a dataset with 2 classes (values) in it's label row, you will get 6 x 2 = 12 samples. If a class does not have 6 samples in it, you will get the max possible for that class. In a scenario where you are selecting 6 samples (-n 6) and have two classes but one class has only 3 samples then you will get 6 + 3 = 9 selected samples.
Stratification: To stratify any method use the --stratify argument which is the first column keyword of the metadata row used to stratify samples before selection occurs. (Selection will occur independently within each strata). This example stratifies diverse selection by the "Label".
python MicroPITA.py --lastmeta Label --stratify Label -m representative input/Test.pcl output.txt python MicroPITA.py --stratify Label -m representative input/Test.biom output.txt
python MicroPITA.py --lastmeta Label --label Label --stratify StratifyLabel -m distinct input/Test.pcl output.txt python MicroPITA.py --label Label --stratify StratifyLabel -m distinct input/Test2.biom output.txt
Changing PCL file defaults: Some PCL files have feature metadata. These are columns of data that comment on bug features (rows) in the file. An example of this could be a certain taxonomy clade for different bug features. If this type of data exists please use -w or --lastFeatureMetadata to indicate the last column of feature metadata before the first column which is a sample. For an example please look in docs for PCL-Description.txt. This only applys to PCL files.
python MicroPITA.py --lastmeta Label -m representative -w taxonomy_5 input/FeatureMetadata.pcl output.txt
MicroPITA assumes the first row of the input file is the sample IDs, if it is not you may use --id to indicate the row. --id expects the entry in the first column of your input file that matches the row used as Sample Ids. See the input file and the following command as an example. This only applys to PCL files.
python MicroPITA.py --id Sample --lastmeta Label -m representative input/Test.pcl output.txt
MicroPITA assumes the input file is TAB delimited, we strongly recommend you use this convention. If not, you can use --delim to change the delimiter used to read in the file. Here is an example of reading the comma delimited file micropita/input/CommaDelim.pcl This only applys to PCL files.
python MicroPITA.py --delim , --lastmeta Label -m representative input/CommaDelim.pcl output.txt
MicroPITA assumes the input file has feature names in which, if the name contains the consensus lineage or full taxonomic hierarchy, it is delimited with a pipe "|". We strongly recommend you use this default. The delimiter of the feature name can be changed using --featdelim. Here is an example of reading in a file with periods as the delimiter. This only applys to PCL files.
python MicroPITA.py --featdelim . --lastmeta Label -m representative input/PeriodDelim.pcl output.txt
Please note the following dependencies need to be installed for micropita to run.
- Python 2.x http://www.python.org/download/
- blist http://pypi.python.org/pypi/blist/
- NumPy http://numpy.scipy.org/
- SciPy http://www.scipy.org/
- PyCogent http://pycogent.sourceforge.net/install.html
- mlpy http://mlpy.sourceforge.net/
- mpi4py http://mpi4py.scipy.org/
- biome support http://biom-format.org/
This covers how to use microPITA. Thank you for using this software and good luck with all your endeavors!
Thanks go to these wonderful people: