From 60ee644d6f1add6722eb8444385f203c6110975a Mon Sep 17 00:00:00 2001 From: Calvin Morrison Date: Wed, 12 Jun 2013 09:10:25 -0400 Subject: Matlab: move multifasta2otu into the main quikr folder --- src/matlab/README_for_multifasta2otu | 27 +++++++++++++++++++++++++++ 1 file changed, 27 insertions(+) create mode 100644 src/matlab/README_for_multifasta2otu (limited to 'src/matlab/README_for_multifasta2otu') diff --git a/src/matlab/README_for_multifasta2otu b/src/matlab/README_for_multifasta2otu new file mode 100644 index 0000000..05a26c6 --- /dev/null +++ b/src/matlab/README_for_multifasta2otu @@ -0,0 +1,27 @@ +* Quikr multifasta->otu_table_(for_qiime_use) wrapper code written by Gail Rosen -- 2/1/2013 + +Usage tips: +* Please name fasta files of sample reads with .fa<*> and place them into one directory without any other file in that directory (for example, no hidden files that the operating system may generate, are allowed in that directory) +* Note: When making your QIIME Metadata file, the sample id's must match the fasta file prefix names +* Fasta files of reads must have a suffix that starts with .fa (e.g.: .fasta and .fa are valid while .fna is NOT) +* Modify the top of the Matlab/Octave scripts for , , , and + +To use with QIIME, one must run the QIIME conversion tool on our OTU table output: +convert_biom.py -i -o .biom --biom_table_type="otu table" + +--------------------------- + +4-step QIIME procedure after using Quikr to obtain 3D PCoA graphs: +(Note: Our code works much better with WEIGHTED Unifrac as opposed to +Unweighted.) + +Pre-requisites: +1. +2. the tree of the database sequences that were used (e.g. rdp7_mafft.fasttree, gg_94_otus_4feb2011.tre. they are in the data directory) +3. your-defined + +1. convert_biom.py -i -o .biom --biom_table_type="otu table" +2. beta_diversity.py -i .biom -m weighted_unifrac -o beta_div -t +3. principal_coordinates.py -i beta_div/weighted_unifrac_.txt -o _weighted.txt +4. make_3d_plots.py -i _weighted.txt -o <3d_pcoa_plotdirectory> -m + -- cgit v1.2.3