miga-base 1.2.15.2 → 1.2.15.4
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- checksums.yaml +4 -4
- data/lib/miga/cli/action/download/gtdb.rb +4 -1
- data/lib/miga/cli/action/gtdb_get.rb +4 -0
- data/lib/miga/daemon.rb +4 -1
- data/lib/miga/lair.rb +6 -4
- data/lib/miga/remote_dataset/download.rb +3 -2
- data/lib/miga/remote_dataset.rb +25 -7
- data/lib/miga/taxonomy.rb +6 -0
- data/lib/miga/version.rb +2 -2
- metadata +6 -302
- data/utils/FastAAI/00.Libraries/01.SCG_HMMs/Archaea_SCG.hmm +0 -41964
- data/utils/FastAAI/00.Libraries/01.SCG_HMMs/Bacteria_SCG.hmm +0 -32439
- data/utils/FastAAI/00.Libraries/01.SCG_HMMs/Complete_SCG_DB.hmm +0 -62056
- data/utils/FastAAI/FastAAI +0 -3659
- data/utils/FastAAI/FastAAI-legacy/FastAAI +0 -1336
- data/utils/FastAAI/FastAAI-legacy/kAAI_v1.0_virus.py +0 -1296
- data/utils/FastAAI/README.md +0 -84
- data/utils/enveomics/Docs/recplot2.md +0 -244
- data/utils/enveomics/Examples/aai-matrix.bash +0 -66
- data/utils/enveomics/Examples/ani-matrix.bash +0 -66
- data/utils/enveomics/Examples/essential-phylogeny.bash +0 -105
- data/utils/enveomics/Examples/unus-genome-phylogeny.bash +0 -100
- data/utils/enveomics/LICENSE.txt +0 -73
- data/utils/enveomics/Makefile +0 -52
- data/utils/enveomics/Manifest/Tasks/aasubs.json +0 -103
- data/utils/enveomics/Manifest/Tasks/blasttab.json +0 -790
- data/utils/enveomics/Manifest/Tasks/distances.json +0 -161
- data/utils/enveomics/Manifest/Tasks/fasta.json +0 -802
- data/utils/enveomics/Manifest/Tasks/fastq.json +0 -291
- data/utils/enveomics/Manifest/Tasks/graphics.json +0 -126
- data/utils/enveomics/Manifest/Tasks/mapping.json +0 -137
- data/utils/enveomics/Manifest/Tasks/ogs.json +0 -382
- data/utils/enveomics/Manifest/Tasks/other.json +0 -906
- data/utils/enveomics/Manifest/Tasks/remote.json +0 -355
- data/utils/enveomics/Manifest/Tasks/sequence-identity.json +0 -650
- data/utils/enveomics/Manifest/Tasks/tables.json +0 -308
- data/utils/enveomics/Manifest/Tasks/trees.json +0 -68
- data/utils/enveomics/Manifest/Tasks/variants.json +0 -111
- data/utils/enveomics/Manifest/categories.json +0 -165
- data/utils/enveomics/Manifest/examples.json +0 -162
- data/utils/enveomics/Manifest/tasks.json +0 -4
- data/utils/enveomics/Pipelines/assembly.pbs/CONFIG.mock.bash +0 -69
- data/utils/enveomics/Pipelines/assembly.pbs/FastA.N50.pl +0 -1
- data/utils/enveomics/Pipelines/assembly.pbs/FastA.filterN.pl +0 -1
- data/utils/enveomics/Pipelines/assembly.pbs/FastA.length.pl +0 -1
- data/utils/enveomics/Pipelines/assembly.pbs/README.md +0 -189
- data/utils/enveomics/Pipelines/assembly.pbs/RUNME-2.bash +0 -112
- data/utils/enveomics/Pipelines/assembly.pbs/RUNME-3.bash +0 -23
- data/utils/enveomics/Pipelines/assembly.pbs/RUNME-4.bash +0 -44
- data/utils/enveomics/Pipelines/assembly.pbs/RUNME.bash +0 -50
- data/utils/enveomics/Pipelines/assembly.pbs/kSelector.R +0 -37
- data/utils/enveomics/Pipelines/assembly.pbs/newbler.pbs +0 -68
- data/utils/enveomics/Pipelines/assembly.pbs/newbler_preparator.pl +0 -49
- data/utils/enveomics/Pipelines/assembly.pbs/soap.pbs +0 -80
- data/utils/enveomics/Pipelines/assembly.pbs/stats.pbs +0 -57
- data/utils/enveomics/Pipelines/assembly.pbs/velvet.pbs +0 -63
- data/utils/enveomics/Pipelines/blast.pbs/01.pbs.bash +0 -38
- data/utils/enveomics/Pipelines/blast.pbs/02.pbs.bash +0 -73
- data/utils/enveomics/Pipelines/blast.pbs/03.pbs.bash +0 -21
- data/utils/enveomics/Pipelines/blast.pbs/BlastTab.recover_job.pl +0 -72
- data/utils/enveomics/Pipelines/blast.pbs/CONFIG.mock.bash +0 -98
- data/utils/enveomics/Pipelines/blast.pbs/FastA.split.pl +0 -1
- data/utils/enveomics/Pipelines/blast.pbs/README.md +0 -127
- data/utils/enveomics/Pipelines/blast.pbs/RUNME.bash +0 -109
- data/utils/enveomics/Pipelines/blast.pbs/TASK.check.bash +0 -128
- data/utils/enveomics/Pipelines/blast.pbs/TASK.dry.bash +0 -16
- data/utils/enveomics/Pipelines/blast.pbs/TASK.eo.bash +0 -22
- data/utils/enveomics/Pipelines/blast.pbs/TASK.pause.bash +0 -26
- data/utils/enveomics/Pipelines/blast.pbs/TASK.run.bash +0 -89
- data/utils/enveomics/Pipelines/blast.pbs/sentinel.pbs.bash +0 -29
- data/utils/enveomics/Pipelines/idba.pbs/README.md +0 -49
- data/utils/enveomics/Pipelines/idba.pbs/RUNME.bash +0 -95
- data/utils/enveomics/Pipelines/idba.pbs/run.pbs +0 -56
- data/utils/enveomics/Pipelines/trim.pbs/README.md +0 -54
- data/utils/enveomics/Pipelines/trim.pbs/RUNME.bash +0 -70
- data/utils/enveomics/Pipelines/trim.pbs/run.pbs +0 -130
- data/utils/enveomics/README.md +0 -42
- data/utils/enveomics/Scripts/AAsubs.log2ratio.rb +0 -171
- data/utils/enveomics/Scripts/Aln.cat.rb +0 -221
- data/utils/enveomics/Scripts/Aln.convert.pl +0 -35
- data/utils/enveomics/Scripts/AlphaDiversity.pl +0 -152
- data/utils/enveomics/Scripts/BedGraph.tad.rb +0 -93
- data/utils/enveomics/Scripts/BedGraph.window.rb +0 -71
- data/utils/enveomics/Scripts/BlastPairwise.AAsubs.pl +0 -102
- data/utils/enveomics/Scripts/BlastTab.addlen.rb +0 -63
- data/utils/enveomics/Scripts/BlastTab.advance.bash +0 -48
- data/utils/enveomics/Scripts/BlastTab.best_hit_sorted.pl +0 -55
- data/utils/enveomics/Scripts/BlastTab.catsbj.pl +0 -104
- data/utils/enveomics/Scripts/BlastTab.cogCat.rb +0 -76
- data/utils/enveomics/Scripts/BlastTab.filter.pl +0 -47
- data/utils/enveomics/Scripts/BlastTab.kegg_pep2path_rest.pl +0 -194
- data/utils/enveomics/Scripts/BlastTab.metaxaPrep.pl +0 -104
- data/utils/enveomics/Scripts/BlastTab.pairedHits.rb +0 -157
- data/utils/enveomics/Scripts/BlastTab.recplot2.R +0 -48
- data/utils/enveomics/Scripts/BlastTab.seqdepth.pl +0 -86
- data/utils/enveomics/Scripts/BlastTab.seqdepth_ZIP.pl +0 -119
- data/utils/enveomics/Scripts/BlastTab.seqdepth_nomedian.pl +0 -86
- data/utils/enveomics/Scripts/BlastTab.subsample.pl +0 -47
- data/utils/enveomics/Scripts/BlastTab.sumPerHit.pl +0 -114
- data/utils/enveomics/Scripts/BlastTab.taxid2taxrank.pl +0 -90
- data/utils/enveomics/Scripts/BlastTab.topHits_sorted.rb +0 -123
- data/utils/enveomics/Scripts/Chao1.pl +0 -97
- data/utils/enveomics/Scripts/CharTable.classify.rb +0 -234
- data/utils/enveomics/Scripts/EBIseq2tax.rb +0 -83
- data/utils/enveomics/Scripts/FastA.N50.pl +0 -60
- data/utils/enveomics/Scripts/FastA.extract.rb +0 -152
- data/utils/enveomics/Scripts/FastA.filter.pl +0 -52
- data/utils/enveomics/Scripts/FastA.filterLen.pl +0 -28
- data/utils/enveomics/Scripts/FastA.filterN.pl +0 -60
- data/utils/enveomics/Scripts/FastA.fragment.rb +0 -100
- data/utils/enveomics/Scripts/FastA.gc.pl +0 -42
- data/utils/enveomics/Scripts/FastA.interpose.pl +0 -93
- data/utils/enveomics/Scripts/FastA.length.pl +0 -38
- data/utils/enveomics/Scripts/FastA.mask.rb +0 -89
- data/utils/enveomics/Scripts/FastA.per_file.pl +0 -36
- data/utils/enveomics/Scripts/FastA.qlen.pl +0 -57
- data/utils/enveomics/Scripts/FastA.rename.pl +0 -65
- data/utils/enveomics/Scripts/FastA.revcom.pl +0 -23
- data/utils/enveomics/Scripts/FastA.sample.rb +0 -98
- data/utils/enveomics/Scripts/FastA.slider.pl +0 -85
- data/utils/enveomics/Scripts/FastA.split.pl +0 -55
- data/utils/enveomics/Scripts/FastA.split.rb +0 -79
- data/utils/enveomics/Scripts/FastA.subsample.pl +0 -131
- data/utils/enveomics/Scripts/FastA.tag.rb +0 -65
- data/utils/enveomics/Scripts/FastA.toFastQ.rb +0 -69
- data/utils/enveomics/Scripts/FastA.wrap.rb +0 -48
- data/utils/enveomics/Scripts/FastQ.filter.pl +0 -54
- data/utils/enveomics/Scripts/FastQ.interpose.pl +0 -90
- data/utils/enveomics/Scripts/FastQ.maskQual.rb +0 -89
- data/utils/enveomics/Scripts/FastQ.offset.pl +0 -90
- data/utils/enveomics/Scripts/FastQ.split.pl +0 -53
- data/utils/enveomics/Scripts/FastQ.tag.rb +0 -70
- data/utils/enveomics/Scripts/FastQ.test-error.rb +0 -81
- data/utils/enveomics/Scripts/FastQ.toFastA.awk +0 -24
- data/utils/enveomics/Scripts/GFF.catsbj.pl +0 -127
- data/utils/enveomics/Scripts/GenBank.add_fields.rb +0 -84
- data/utils/enveomics/Scripts/HMM.essential.rb +0 -351
- data/utils/enveomics/Scripts/HMM.haai.rb +0 -168
- data/utils/enveomics/Scripts/HMMsearch.extractIds.rb +0 -83
- data/utils/enveomics/Scripts/JPlace.distances.rb +0 -88
- data/utils/enveomics/Scripts/JPlace.to_iToL.rb +0 -320
- data/utils/enveomics/Scripts/M5nr.getSequences.rb +0 -81
- data/utils/enveomics/Scripts/MeTaxa.distribution.pl +0 -198
- data/utils/enveomics/Scripts/MyTaxa.fragsByTax.pl +0 -35
- data/utils/enveomics/Scripts/MyTaxa.seq-taxrank.rb +0 -49
- data/utils/enveomics/Scripts/NCBIacc2tax.rb +0 -92
- data/utils/enveomics/Scripts/Newick.autoprune.R +0 -27
- data/utils/enveomics/Scripts/RAxML-EPA.to_iToL.pl +0 -228
- data/utils/enveomics/Scripts/RecPlot2.compareIdentities.R +0 -32
- data/utils/enveomics/Scripts/RefSeq.download.bash +0 -48
- data/utils/enveomics/Scripts/SRA.download.bash +0 -55
- data/utils/enveomics/Scripts/TRIBS.plot-test.R +0 -36
- data/utils/enveomics/Scripts/TRIBS.test.R +0 -39
- data/utils/enveomics/Scripts/Table.barplot.R +0 -31
- data/utils/enveomics/Scripts/Table.df2dist.R +0 -30
- data/utils/enveomics/Scripts/Table.filter.pl +0 -61
- data/utils/enveomics/Scripts/Table.merge.pl +0 -77
- data/utils/enveomics/Scripts/Table.prefScore.R +0 -60
- data/utils/enveomics/Scripts/Table.replace.rb +0 -69
- data/utils/enveomics/Scripts/Table.round.rb +0 -63
- data/utils/enveomics/Scripts/Table.split.pl +0 -57
- data/utils/enveomics/Scripts/Taxonomy.silva2ncbi.rb +0 -227
- data/utils/enveomics/Scripts/VCF.KaKs.rb +0 -147
- data/utils/enveomics/Scripts/VCF.SNPs.rb +0 -88
- data/utils/enveomics/Scripts/aai.rb +0 -421
- data/utils/enveomics/Scripts/ani.rb +0 -362
- data/utils/enveomics/Scripts/anir.rb +0 -137
- data/utils/enveomics/Scripts/clust.rand.rb +0 -102
- data/utils/enveomics/Scripts/gi2tax.rb +0 -103
- data/utils/enveomics/Scripts/in_silico_GA_GI.pl +0 -96
- data/utils/enveomics/Scripts/lib/data/dupont_2012_essential.hmm.gz +0 -0
- data/utils/enveomics/Scripts/lib/data/lee_2019_essential.hmm.gz +0 -0
- data/utils/enveomics/Scripts/lib/enveomics.R +0 -1
- data/utils/enveomics/Scripts/lib/enveomics_rb/anir.rb +0 -293
- data/utils/enveomics/Scripts/lib/enveomics_rb/bm_set.rb +0 -175
- data/utils/enveomics/Scripts/lib/enveomics_rb/enveomics.rb +0 -24
- data/utils/enveomics/Scripts/lib/enveomics_rb/errors.rb +0 -17
- data/utils/enveomics/Scripts/lib/enveomics_rb/gmm_em.rb +0 -30
- data/utils/enveomics/Scripts/lib/enveomics_rb/jplace.rb +0 -253
- data/utils/enveomics/Scripts/lib/enveomics_rb/match.rb +0 -88
- data/utils/enveomics/Scripts/lib/enveomics_rb/og.rb +0 -182
- data/utils/enveomics/Scripts/lib/enveomics_rb/rbm.rb +0 -49
- data/utils/enveomics/Scripts/lib/enveomics_rb/remote_data.rb +0 -74
- data/utils/enveomics/Scripts/lib/enveomics_rb/seq_range.rb +0 -237
- data/utils/enveomics/Scripts/lib/enveomics_rb/stats/rand.rb +0 -31
- data/utils/enveomics/Scripts/lib/enveomics_rb/stats/sample.rb +0 -152
- data/utils/enveomics/Scripts/lib/enveomics_rb/stats.rb +0 -3
- data/utils/enveomics/Scripts/lib/enveomics_rb/utils.rb +0 -74
- data/utils/enveomics/Scripts/lib/enveomics_rb/vcf.rb +0 -135
- data/utils/enveomics/Scripts/ogs.annotate.rb +0 -88
- data/utils/enveomics/Scripts/ogs.core-pan.rb +0 -160
- data/utils/enveomics/Scripts/ogs.extract.rb +0 -125
- data/utils/enveomics/Scripts/ogs.mcl.rb +0 -186
- data/utils/enveomics/Scripts/ogs.rb +0 -104
- data/utils/enveomics/Scripts/ogs.stats.rb +0 -131
- data/utils/enveomics/Scripts/rbm-legacy.rb +0 -172
- data/utils/enveomics/Scripts/rbm.rb +0 -108
- data/utils/enveomics/Scripts/sam.filter.rb +0 -148
- data/utils/enveomics/Tests/Makefile +0 -10
- data/utils/enveomics/Tests/Mgen_M2288.faa +0 -3189
- data/utils/enveomics/Tests/Mgen_M2288.fna +0 -8282
- data/utils/enveomics/Tests/Mgen_M2321.fna +0 -8288
- data/utils/enveomics/Tests/Nequ_Kin4M.faa +0 -2970
- data/utils/enveomics/Tests/Xanthomonas_oryzae-PilA.tribs.Rdata +0 -0
- data/utils/enveomics/Tests/Xanthomonas_oryzae-PilA.txt +0 -7
- data/utils/enveomics/Tests/Xanthomonas_oryzae.aai-mat.tsv +0 -17
- data/utils/enveomics/Tests/Xanthomonas_oryzae.aai.tsv +0 -137
- data/utils/enveomics/Tests/a_mg.cds-go.blast.tsv +0 -123
- data/utils/enveomics/Tests/a_mg.reads-cds.blast.tsv +0 -200
- data/utils/enveomics/Tests/a_mg.reads-cds.counts.tsv +0 -55
- data/utils/enveomics/Tests/alkB.nwk +0 -1
- data/utils/enveomics/Tests/anthrax-cansnp-data.tsv +0 -13
- data/utils/enveomics/Tests/anthrax-cansnp-key.tsv +0 -17
- data/utils/enveomics/Tests/hiv1.faa +0 -59
- data/utils/enveomics/Tests/hiv1.fna +0 -134
- data/utils/enveomics/Tests/hiv2.faa +0 -70
- data/utils/enveomics/Tests/hiv_mix-hiv1.blast.tsv +0 -233
- data/utils/enveomics/Tests/hiv_mix-hiv1.blast.tsv.lim +0 -1
- data/utils/enveomics/Tests/hiv_mix-hiv1.blast.tsv.rec +0 -233
- data/utils/enveomics/Tests/phyla_counts.tsv +0 -10
- data/utils/enveomics/Tests/primate_lentivirus.ogs +0 -11
- data/utils/enveomics/Tests/primate_lentivirus.rbm/hiv1-hiv1.rbm +0 -9
- data/utils/enveomics/Tests/primate_lentivirus.rbm/hiv1-hiv2.rbm +0 -8
- data/utils/enveomics/Tests/primate_lentivirus.rbm/hiv1-siv.rbm +0 -6
- data/utils/enveomics/Tests/primate_lentivirus.rbm/hiv2-hiv2.rbm +0 -9
- data/utils/enveomics/Tests/primate_lentivirus.rbm/hiv2-siv.rbm +0 -6
- data/utils/enveomics/Tests/primate_lentivirus.rbm/siv-siv.rbm +0 -6
- data/utils/enveomics/build_enveomics_r.bash +0 -45
- data/utils/enveomics/enveomics.R/DESCRIPTION +0 -31
- data/utils/enveomics/enveomics.R/NAMESPACE +0 -39
- data/utils/enveomics/enveomics.R/R/autoprune.R +0 -155
- data/utils/enveomics/enveomics.R/R/barplot.R +0 -184
- data/utils/enveomics/enveomics.R/R/cliopts.R +0 -135
- data/utils/enveomics/enveomics.R/R/df2dist.R +0 -154
- data/utils/enveomics/enveomics.R/R/growthcurve.R +0 -331
- data/utils/enveomics/enveomics.R/R/prefscore.R +0 -79
- data/utils/enveomics/enveomics.R/R/recplot.R +0 -354
- data/utils/enveomics/enveomics.R/R/recplot2.R +0 -1631
- data/utils/enveomics/enveomics.R/R/tribs.R +0 -583
- data/utils/enveomics/enveomics.R/R/utils.R +0 -80
- data/utils/enveomics/enveomics.R/README.md +0 -81
- data/utils/enveomics/enveomics.R/data/growth.curves.rda +0 -0
- data/utils/enveomics/enveomics.R/data/phyla.counts.rda +0 -0
- data/utils/enveomics/enveomics.R/man/cash-enve.GrowthCurve-method.Rd +0 -16
- data/utils/enveomics/enveomics.R/man/cash-enve.RecPlot2-method.Rd +0 -16
- data/utils/enveomics/enveomics.R/man/cash-enve.RecPlot2.Peak-method.Rd +0 -16
- data/utils/enveomics/enveomics.R/man/enve.GrowthCurve-class.Rd +0 -25
- data/utils/enveomics/enveomics.R/man/enve.TRIBS-class.Rd +0 -46
- data/utils/enveomics/enveomics.R/man/enve.TRIBS.merge.Rd +0 -23
- data/utils/enveomics/enveomics.R/man/enve.TRIBStest-class.Rd +0 -47
- data/utils/enveomics/enveomics.R/man/enve.__prune.iter.Rd +0 -23
- data/utils/enveomics/enveomics.R/man/enve.__prune.reduce.Rd +0 -23
- data/utils/enveomics/enveomics.R/man/enve.__tribs.Rd +0 -40
- data/utils/enveomics/enveomics.R/man/enve.barplot.Rd +0 -103
- data/utils/enveomics/enveomics.R/man/enve.cliopts.Rd +0 -67
- data/utils/enveomics/enveomics.R/man/enve.col.alpha.Rd +0 -24
- data/utils/enveomics/enveomics.R/man/enve.col2alpha.Rd +0 -19
- data/utils/enveomics/enveomics.R/man/enve.df2dist.Rd +0 -45
- data/utils/enveomics/enveomics.R/man/enve.df2dist.group.Rd +0 -44
- data/utils/enveomics/enveomics.R/man/enve.df2dist.list.Rd +0 -47
- data/utils/enveomics/enveomics.R/man/enve.growthcurve.Rd +0 -75
- data/utils/enveomics/enveomics.R/man/enve.prefscore.Rd +0 -50
- data/utils/enveomics/enveomics.R/man/enve.prune.dist.Rd +0 -44
- data/utils/enveomics/enveomics.R/man/enve.recplot.Rd +0 -139
- data/utils/enveomics/enveomics.R/man/enve.recplot2-class.Rd +0 -45
- data/utils/enveomics/enveomics.R/man/enve.recplot2.ANIr.Rd +0 -24
- data/utils/enveomics/enveomics.R/man/enve.recplot2.Rd +0 -77
- data/utils/enveomics/enveomics.R/man/enve.recplot2.__counts.Rd +0 -25
- data/utils/enveomics/enveomics.R/man/enve.recplot2.__peakHist.Rd +0 -21
- data/utils/enveomics/enveomics.R/man/enve.recplot2.__whichClosestPeak.Rd +0 -19
- data/utils/enveomics/enveomics.R/man/enve.recplot2.changeCutoff.Rd +0 -19
- data/utils/enveomics/enveomics.R/man/enve.recplot2.compareIdentities.Rd +0 -47
- data/utils/enveomics/enveomics.R/man/enve.recplot2.coordinates.Rd +0 -29
- data/utils/enveomics/enveomics.R/man/enve.recplot2.corePeak.Rd +0 -18
- data/utils/enveomics/enveomics.R/man/enve.recplot2.extractWindows.Rd +0 -45
- data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.Rd +0 -36
- data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.__em_e.Rd +0 -19
- data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.__em_m.Rd +0 -19
- data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.__emauto_one.Rd +0 -27
- data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.__mow_one.Rd +0 -52
- data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.__mower.Rd +0 -17
- data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.em.Rd +0 -51
- data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.emauto.Rd +0 -43
- data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.mower.Rd +0 -82
- data/utils/enveomics/enveomics.R/man/enve.recplot2.peak-class.Rd +0 -59
- data/utils/enveomics/enveomics.R/man/enve.recplot2.seqdepth.Rd +0 -27
- data/utils/enveomics/enveomics.R/man/enve.recplot2.windowDepthThreshold.Rd +0 -36
- data/utils/enveomics/enveomics.R/man/enve.selvector.Rd +0 -23
- data/utils/enveomics/enveomics.R/man/enve.tribs.Rd +0 -68
- data/utils/enveomics/enveomics.R/man/enve.tribs.test.Rd +0 -28
- data/utils/enveomics/enveomics.R/man/enve.truncate.Rd +0 -27
- data/utils/enveomics/enveomics.R/man/growth.curves.Rd +0 -14
- data/utils/enveomics/enveomics.R/man/phyla.counts.Rd +0 -13
- data/utils/enveomics/enveomics.R/man/plot.enve.GrowthCurve.Rd +0 -78
- data/utils/enveomics/enveomics.R/man/plot.enve.TRIBS.Rd +0 -46
- data/utils/enveomics/enveomics.R/man/plot.enve.TRIBStest.Rd +0 -45
- data/utils/enveomics/enveomics.R/man/plot.enve.recplot2.Rd +0 -125
- data/utils/enveomics/enveomics.R/man/summary.enve.GrowthCurve.Rd +0 -19
- data/utils/enveomics/enveomics.R/man/summary.enve.TRIBS.Rd +0 -19
- data/utils/enveomics/enveomics.R/man/summary.enve.TRIBStest.Rd +0 -19
- data/utils/enveomics/globals.mk +0 -8
- data/utils/enveomics/manifest.json +0 -9
- data/utils/multitrim/Multitrim How-To.pdf +0 -0
- data/utils/multitrim/README.md +0 -67
- data/utils/multitrim/multitrim.py +0 -1555
- data/utils/multitrim/multitrim.yml +0 -13
|
@@ -1,47 +0,0 @@
|
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
|
2
|
-
% Please edit documentation in R/df2dist.R
|
|
3
|
-
\name{enve.df2dist.list}
|
|
4
|
-
\alias{enve.df2dist.list}
|
|
5
|
-
\title{Enveomics: Data Frame to Dist (List)}
|
|
6
|
-
\usage{
|
|
7
|
-
enve.df2dist.list(
|
|
8
|
-
x,
|
|
9
|
-
groups,
|
|
10
|
-
obj1.index = 1,
|
|
11
|
-
obj2.index = 2,
|
|
12
|
-
dist.index = 3,
|
|
13
|
-
empty.rm = TRUE,
|
|
14
|
-
...
|
|
15
|
-
)
|
|
16
|
-
}
|
|
17
|
-
\arguments{
|
|
18
|
-
\item{x}{A dataframe (or coercible object) with at least three columns:
|
|
19
|
-
\enumerate{
|
|
20
|
-
\item ID of the object 1,
|
|
21
|
-
\item ID of the object 2, and
|
|
22
|
-
\item distance between the two objects.}}
|
|
23
|
-
|
|
24
|
-
\item{groups}{Named array where the IDs correspond to the object IDs,
|
|
25
|
-
and the values correspond to the group.}
|
|
26
|
-
|
|
27
|
-
\item{obj1.index}{Index of the column containing the ID of the object 1.}
|
|
28
|
-
|
|
29
|
-
\item{obj2.index}{Index of the column containing the ID of the object 2.}
|
|
30
|
-
|
|
31
|
-
\item{dist.index}{Index of the column containing the distance.}
|
|
32
|
-
|
|
33
|
-
\item{empty.rm}{Remove incomplete matrices.}
|
|
34
|
-
|
|
35
|
-
\item{...}{Any other parameters supported by
|
|
36
|
-
\code{\link{enve.df2dist.group}}.}
|
|
37
|
-
}
|
|
38
|
-
\value{
|
|
39
|
-
Returns a \strong{list} of \strong{dist} objects.
|
|
40
|
-
}
|
|
41
|
-
\description{
|
|
42
|
-
Transform a dataframe (or coercible object, like a table)
|
|
43
|
-
into a \strong{dist} object.
|
|
44
|
-
}
|
|
45
|
-
\author{
|
|
46
|
-
Luis M. Rodriguez-R [aut, cre]
|
|
47
|
-
}
|
|
@@ -1,75 +0,0 @@
|
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
|
2
|
-
% Please edit documentation in R/growthcurve.R
|
|
3
|
-
\name{enve.growthcurve}
|
|
4
|
-
\alias{enve.growthcurve}
|
|
5
|
-
\title{Enveomics: Growth Curve}
|
|
6
|
-
\usage{
|
|
7
|
-
enve.growthcurve(
|
|
8
|
-
x,
|
|
9
|
-
times = 1:nrow(x),
|
|
10
|
-
triplicates = FALSE,
|
|
11
|
-
design,
|
|
12
|
-
new.times = seq(min(times), max(times), length.out = length(times) * 10),
|
|
13
|
-
level = 0.95,
|
|
14
|
-
interval = c("confidence", "prediction"),
|
|
15
|
-
plot = TRUE,
|
|
16
|
-
FUN = function(t, K, r, P0) K * P0 * exp(r * t)/(K + P0 * (exp(r * t) - 1)),
|
|
17
|
-
nls.opt = list(),
|
|
18
|
-
...
|
|
19
|
-
)
|
|
20
|
-
}
|
|
21
|
-
\arguments{
|
|
22
|
-
\item{x}{Data frame (or coercible) containing the observed growth data
|
|
23
|
-
(e.g., O.D. values). Each column is an independent growth curve and each
|
|
24
|
-
row is a time point. \code{NA}'s are allowed.}
|
|
25
|
-
|
|
26
|
-
\item{times}{Vector with the times at which each row was taken. By default,
|
|
27
|
-
all rows are assumed to be part of constantly periodic measurements.}
|
|
28
|
-
|
|
29
|
-
\item{triplicates}{If \code{TRUE}, the columns are assumed to be sorted by
|
|
30
|
-
sample with three replicates by sample. It requires a number of columns
|
|
31
|
-
multiple of 3.}
|
|
32
|
-
|
|
33
|
-
\item{design}{Experimental design of the data. An \strong{array} of mode list
|
|
34
|
-
with sample names as index and the list of column names in each sample as
|
|
35
|
-
the values. By default, each column is assumed to be an independent sample
|
|
36
|
-
if \code{triplicates} is \code{FALSE}, or every three columns are assumed
|
|
37
|
-
to be a sample if \code{triplicates} is \code{TRUE}. In the latter case,
|
|
38
|
-
samples are simply numbered.}
|
|
39
|
-
|
|
40
|
-
\item{new.times}{Values of time for the fitted curve.}
|
|
41
|
-
|
|
42
|
-
\item{level}{Confidence (or prediction) interval in the fitted curve.}
|
|
43
|
-
|
|
44
|
-
\item{interval}{Type of interval to be calculated for the fitted curve.}
|
|
45
|
-
|
|
46
|
-
\item{plot}{Should the growth curve be plotted?}
|
|
47
|
-
|
|
48
|
-
\item{FUN}{Function to fit. By default: logistic growth with paramenters
|
|
49
|
-
\code{K}: carrying capacity,
|
|
50
|
-
\code{r}: intrinsic growth rate, and
|
|
51
|
-
\code{P0}: Initial population.}
|
|
52
|
-
|
|
53
|
-
\item{nls.opt}{Any additional options passed to \code{nls}.}
|
|
54
|
-
|
|
55
|
-
\item{...}{Any additional parameters to be passed to
|
|
56
|
-
\code{plot.enve.GrowthCurve}.}
|
|
57
|
-
}
|
|
58
|
-
\value{
|
|
59
|
-
Returns an \code{\link{enve.GrowthCurve}} object.
|
|
60
|
-
}
|
|
61
|
-
\description{
|
|
62
|
-
Calculates growth curves using the logistic growth function.
|
|
63
|
-
}
|
|
64
|
-
\examples{
|
|
65
|
-
# Load data
|
|
66
|
-
data("growth.curves", package="enveomics.R", envir=environment())
|
|
67
|
-
# Generate growth curves with different colors
|
|
68
|
-
g <- enve.growthcurve(growth.curves[,-1], growth.curves[,1], triplicates=TRUE)
|
|
69
|
-
# Generate black-and-white growth curves with different symbols
|
|
70
|
-
plot(g, pch=15:17, col="black", band.density=45, band.angle=c(-45,45,0))
|
|
71
|
-
|
|
72
|
-
}
|
|
73
|
-
\author{
|
|
74
|
-
Luis M. Rodriguez-R [aut, cre]
|
|
75
|
-
}
|
|
@@ -1,50 +0,0 @@
|
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
|
2
|
-
% Please edit documentation in R/prefscore.R
|
|
3
|
-
\name{enve.prefscore}
|
|
4
|
-
\alias{enve.prefscore}
|
|
5
|
-
\title{Enveomics: Pref Score}
|
|
6
|
-
\usage{
|
|
7
|
-
enve.prefscore(
|
|
8
|
-
x,
|
|
9
|
-
set,
|
|
10
|
-
ignore = NULL,
|
|
11
|
-
signif.thr,
|
|
12
|
-
plot = TRUE,
|
|
13
|
-
col.above = rgb(148, 17, 0, maxColorValue = 255),
|
|
14
|
-
col.equal = rgb(189, 189, 189, maxColorValue = 255),
|
|
15
|
-
col.below = rgb(47, 84, 150, maxColorValue = 255),
|
|
16
|
-
...
|
|
17
|
-
)
|
|
18
|
-
}
|
|
19
|
-
\arguments{
|
|
20
|
-
\item{x}{Occupancy matrix (logical or numeric binary) with species as rows and samples
|
|
21
|
-
as columns}
|
|
22
|
-
|
|
23
|
-
\item{set}{Vector indicating samples in the test set. It can be any selection vector:
|
|
24
|
-
boolean (same length as the number of columns in \code{x}), or numeric or
|
|
25
|
-
character vector with indexes of the \code{x} columns.}
|
|
26
|
-
|
|
27
|
-
\item{ignore}{Vector indicating species to ignore. It can be any selection vector with
|
|
28
|
-
respect to the rows in \code{x} (see \code{set}).}
|
|
29
|
-
|
|
30
|
-
\item{signif.thr}{Absolute value of the significance threshold}
|
|
31
|
-
|
|
32
|
-
\item{plot}{Indicates if a plot should be generated}
|
|
33
|
-
|
|
34
|
-
\item{col.above}{Color for points significantly above zero}
|
|
35
|
-
|
|
36
|
-
\item{col.equal}{Color for points not significantly different from zero}
|
|
37
|
-
|
|
38
|
-
\item{col.below}{Color for points significantly below zero}
|
|
39
|
-
|
|
40
|
-
\item{...}{Any additional parameters supported by \code{plot}}
|
|
41
|
-
}
|
|
42
|
-
\value{
|
|
43
|
-
Returns a named vector of preference scores.
|
|
44
|
-
}
|
|
45
|
-
\description{
|
|
46
|
-
Estimate preference score of species based on occupancy in biased sample sets
|
|
47
|
-
}
|
|
48
|
-
\author{
|
|
49
|
-
Luis M. Rodriguez-R [aut, cre]
|
|
50
|
-
}
|
|
@@ -1,44 +0,0 @@
|
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
|
2
|
-
% Please edit documentation in R/autoprune.R
|
|
3
|
-
\name{enve.prune.dist}
|
|
4
|
-
\alias{enve.prune.dist}
|
|
5
|
-
\title{Enveomics: Prune Dist}
|
|
6
|
-
\usage{
|
|
7
|
-
enve.prune.dist(
|
|
8
|
-
t,
|
|
9
|
-
dist.quantile = 0.25,
|
|
10
|
-
min_dist,
|
|
11
|
-
quiet = FALSE,
|
|
12
|
-
max_iters = 100,
|
|
13
|
-
min_nodes_random = 40000,
|
|
14
|
-
random_nodes_frx = 1
|
|
15
|
-
)
|
|
16
|
-
}
|
|
17
|
-
\arguments{
|
|
18
|
-
\item{t}{A \strong{phylo} object or a path to the Newick file.}
|
|
19
|
-
|
|
20
|
-
\item{dist.quantile}{The quantile of edge lengths.}
|
|
21
|
-
|
|
22
|
-
\item{min_dist}{The minimum distance to allow between two tips.
|
|
23
|
-
If not set, \code{dist.quantile} is used instead to calculate it.}
|
|
24
|
-
|
|
25
|
-
\item{quiet}{Boolean indicating if the function must run without output.}
|
|
26
|
-
|
|
27
|
-
\item{max_iters}{Maximum number of iterations.}
|
|
28
|
-
|
|
29
|
-
\item{min_nodes_random}{Minimum number of nodes to trigger \emph{tip-pairs} nodes sampling.
|
|
30
|
-
This sampling is less reproducible and more computationally expensive,
|
|
31
|
-
but it's the only solution if the cophenetic matrix exceeds \code{2^31-1}
|
|
32
|
-
entries; above that, it cannot be represented in R.}
|
|
33
|
-
|
|
34
|
-
\item{random_nodes_frx}{Fraction of the nodes to be sampled if more than \code{min_nodes_random}.}
|
|
35
|
-
}
|
|
36
|
-
\value{
|
|
37
|
-
Returns a pruned \strong{phylo} object.
|
|
38
|
-
}
|
|
39
|
-
\description{
|
|
40
|
-
Automatically prunes a tree, to keep representatives of each clade.
|
|
41
|
-
}
|
|
42
|
-
\author{
|
|
43
|
-
Luis M. Rodriguez-R [aut, cre]
|
|
44
|
-
}
|
|
@@ -1,139 +0,0 @@
|
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
|
2
|
-
% Please edit documentation in R/recplot.R
|
|
3
|
-
\name{enve.recplot}
|
|
4
|
-
\alias{enve.recplot}
|
|
5
|
-
\title{Enveomics: Recruitment Plots}
|
|
6
|
-
\usage{
|
|
7
|
-
enve.recplot(
|
|
8
|
-
prefix,
|
|
9
|
-
id.min = NULL,
|
|
10
|
-
id.max = NULL,
|
|
11
|
-
id.binsize = NULL,
|
|
12
|
-
id.splines = 0,
|
|
13
|
-
id.metric = "id",
|
|
14
|
-
id.summary = "sum",
|
|
15
|
-
pos.min = 1,
|
|
16
|
-
pos.max = NULL,
|
|
17
|
-
pos.binsize = 1000,
|
|
18
|
-
pos.splines = 0,
|
|
19
|
-
rec.col1 = "white",
|
|
20
|
-
rec.col2 = "black",
|
|
21
|
-
main = NULL,
|
|
22
|
-
contig.col = grey(0.85),
|
|
23
|
-
ret.recplot = FALSE,
|
|
24
|
-
ret.hist = FALSE,
|
|
25
|
-
ret.mode = FALSE,
|
|
26
|
-
id.cutoff = NULL,
|
|
27
|
-
verbose = TRUE,
|
|
28
|
-
...
|
|
29
|
-
)
|
|
30
|
-
}
|
|
31
|
-
\arguments{
|
|
32
|
-
\item{prefix}{Path to the prefix of the BlastTab.catsbj.pl output files. At
|
|
33
|
-
least the files \strong{.rec} and \strong{.lim} must exist with this prefix.}
|
|
34
|
-
|
|
35
|
-
\item{id.min}{Minimum identity to be considered. By default, the minimum detected
|
|
36
|
-
identity. This value is a percentage.}
|
|
37
|
-
|
|
38
|
-
\item{id.max}{Maximum identity to be considered. By default, 100\%.}
|
|
39
|
-
|
|
40
|
-
\item{id.binsize}{Size of the identity bins (vertical histograms). By default, 0.1 for
|
|
41
|
-
identity metrics and 5 for bit score.}
|
|
42
|
-
|
|
43
|
-
\item{id.splines}{Smoothing parameter for the splines in the identity histogram. Zero (0) for no
|
|
44
|
-
splines. A generally good value is 1/2. If non-zero, requires the \pkg{stats} package.}
|
|
45
|
-
|
|
46
|
-
\item{id.metric}{Metric of identity to be used (Y-axis).
|
|
47
|
-
It can be any unambiguous prefix of:
|
|
48
|
-
\itemize{
|
|
49
|
-
\item "identity"
|
|
50
|
-
\item "corrected identity"
|
|
51
|
-
\item "bit score"}}
|
|
52
|
-
|
|
53
|
-
\item{id.summary}{Method used to build the identity histogram (Horizontal axis of the right panel).
|
|
54
|
-
It can be any unambiguous prefix of:
|
|
55
|
-
\itemize{
|
|
56
|
-
\item "sum"
|
|
57
|
-
\item "average"
|
|
58
|
-
\item "median"
|
|
59
|
-
\item "90\% lower bound"
|
|
60
|
-
\item "90\% upper bound"
|
|
61
|
-
\item "95\% lower bound"
|
|
62
|
-
\item "95\% upper bound" }
|
|
63
|
-
The last four options
|
|
64
|
-
correspond to the upper and lower boundaries of the 90\% and 95\% empirical confidence
|
|
65
|
-
intervals.}
|
|
66
|
-
|
|
67
|
-
\item{pos.min}{Minimum (leftmost) position in the reference (concatenated) genome (in bp).}
|
|
68
|
-
|
|
69
|
-
\item{pos.max}{Maximum (rightmost) position in the reference (concatenated) genome (in bp).
|
|
70
|
-
By default: Length of the genome.}
|
|
71
|
-
|
|
72
|
-
\item{pos.binsize}{Size of the position bins (horizontal histograms) in bp.}
|
|
73
|
-
|
|
74
|
-
\item{pos.splines}{Smoothing parameter for the splines in the position histogram. Zero (0) for no splines.
|
|
75
|
-
If non-zero, requires the stats package.}
|
|
76
|
-
|
|
77
|
-
\item{rec.col1}{Lightest color in the recruitment plot.}
|
|
78
|
-
|
|
79
|
-
\item{rec.col2}{Darkest color in the recruitment plot.}
|
|
80
|
-
|
|
81
|
-
\item{main}{Title of the plot.}
|
|
82
|
-
|
|
83
|
-
\item{contig.col}{Color of the Contig boundaries. Set to \code{NA} to ignore Contig boundaries.}
|
|
84
|
-
|
|
85
|
-
\item{ret.recplot}{Indicates if the matrix of the recruitment plot is to be returned.}
|
|
86
|
-
|
|
87
|
-
\item{ret.hist}{Ignored, for backwards compatibility.}
|
|
88
|
-
|
|
89
|
-
\item{ret.mode}{Indicates if the mode of the identity is to be computed. It requires the
|
|
90
|
-
\pkg{modeest} package.}
|
|
91
|
-
|
|
92
|
-
\item{id.cutoff}{Minimum identity to consider an alignment as "top". By default, it is 0.95 for the
|
|
93
|
-
identity metrics and 95\% of the best scoring alignment for bit score.}
|
|
94
|
-
|
|
95
|
-
\item{verbose}{Indicates if the function should report the advance.}
|
|
96
|
-
|
|
97
|
-
\item{...}{Any additional graphic parameters to be passed to plot for all panels except the
|
|
98
|
-
recruitment plot (lower-left).}
|
|
99
|
-
}
|
|
100
|
-
\value{
|
|
101
|
-
Returns a list with the following elements:
|
|
102
|
-
|
|
103
|
-
\describe{
|
|
104
|
-
\item{\code{pos.marks}}{Midpoints of the position histogram.}
|
|
105
|
-
\item{\code{id.matrix}}{Midpoints of the identity histogram.}
|
|
106
|
-
\item{\code{recplot}}{Matrix containing the recruitment plot values
|
|
107
|
-
(if \code{ret.recplot=TRUE}).}
|
|
108
|
-
\item{\code{id.mean}}{Mean identity.}
|
|
109
|
-
\item{\code{id.median}}{Median identity.}
|
|
110
|
-
\item{\code{id.mode}}{Mode of the identity (if \code{ret.mode=TRUE}). Deprecated.}
|
|
111
|
-
\item{\code{id.hist}}{Values of the identity histogram (if \code{ret.hist=TRUE}).}
|
|
112
|
-
\item{\code{pos.hist.low}}{Values of the position histogram (depth) with "low"
|
|
113
|
-
identity (i.e., below id.cutoff) (if \code{ret.hist=TRUE}).}
|
|
114
|
-
\item{\code{pos.hist.top}}{Values of the position histogram (depth) with "top"
|
|
115
|
-
identity (i.e., above id.cutoff) (if \code{ret.hist=TRUE}).}
|
|
116
|
-
\item{\code{id.max}}{Value of \code{id.max}. This is returned because
|
|
117
|
-
\code{id.max=NULL} may vary.}
|
|
118
|
-
\item{\code{id.cutoff}}{Value of \code{id.cutoff}.
|
|
119
|
-
This is returned because \code{id.cutoff=NULL} may vary.}
|
|
120
|
-
\item{\code{seqdepth.mean.top}}{Average sequencing depth with identity above
|
|
121
|
-
\code{id.cutoff}.}
|
|
122
|
-
\item{\code{seqdepth.mean.low}}{Average sequencing depth with identity below
|
|
123
|
-
\code{id.cutoff}.}
|
|
124
|
-
\item{\code{seqdepth.mean.all}}{Average sequencing depth without identity filtering.}
|
|
125
|
-
\item{\code{seqdepth.median.top}}{Median sequencing depth with identity above
|
|
126
|
-
\code{id.cutoff}.}
|
|
127
|
-
\item{\code{seqdepth.median.low}}{Median sequencing depth with identity below
|
|
128
|
-
\code{id.cutoff}.}
|
|
129
|
-
\item{\code{seqdepth.median.all}}{Median sequencing depth without identity filtering.}
|
|
130
|
-
\item{\code{id.metric}}{Full name of the used identity metric.}
|
|
131
|
-
\item{\code{id.summary}}{Full name of the summary method used to build the identity plot.}}
|
|
132
|
-
}
|
|
133
|
-
\description{
|
|
134
|
-
Produces recruitment plots provided that BlastTab.catsbj.pl has
|
|
135
|
-
been previously executed. Requires the \pkg{gplots} library.
|
|
136
|
-
}
|
|
137
|
-
\author{
|
|
138
|
-
Luis M. Rodriguez-R [aut, cre]
|
|
139
|
-
}
|
|
@@ -1,45 +0,0 @@
|
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
|
2
|
-
% Please edit documentation in R/recplot2.R
|
|
3
|
-
\docType{class}
|
|
4
|
-
\name{enve.RecPlot2-class}
|
|
5
|
-
\alias{enve.RecPlot2-class}
|
|
6
|
-
\alias{enve.RecPlot2}
|
|
7
|
-
\title{Enveomics: Recruitment Plot (2) - S4 Class}
|
|
8
|
-
\description{
|
|
9
|
-
Enve-omics representation of Recruitment plots. This object can
|
|
10
|
-
be produced by \code{\link{enve.recplot2}} and supports S4 method plot.
|
|
11
|
-
}
|
|
12
|
-
\section{Slots}{
|
|
13
|
-
|
|
14
|
-
\describe{
|
|
15
|
-
\item{\code{counts}}{\code{(matrix)} Counts as a two-dimensional histogram.}
|
|
16
|
-
|
|
17
|
-
\item{\code{pos.counts.in}}{\code{(numeric)} Counts of in-group hits per position bin.}
|
|
18
|
-
|
|
19
|
-
\item{\code{pos.counts.out}}{\code{(numeric)} Counts of out-group hits per position bin.}
|
|
20
|
-
|
|
21
|
-
\item{\code{id.counts}}{\code{(numeric)} Counts per ID bin.}
|
|
22
|
-
|
|
23
|
-
\item{\code{id.breaks}}{\code{(numeric)} Breaks of identity bins.}
|
|
24
|
-
|
|
25
|
-
\item{\code{pos.breaks}}{\code{(numeric)} Breaks of position bins.}
|
|
26
|
-
|
|
27
|
-
\item{\code{pos.names}}{\code{(character)} Names of the position bins.}
|
|
28
|
-
|
|
29
|
-
\item{\code{seq.breaks}}{\code{(numeric)} Breaks of input sequences.}
|
|
30
|
-
|
|
31
|
-
\item{\code{peaks}}{\code{(list)} Peaks identified in the recplot.
|
|
32
|
-
Limits of the subject sequences after concatenation.}
|
|
33
|
-
|
|
34
|
-
\item{\code{seq.names}}{\code{(character}) Names of the subject sequences.}
|
|
35
|
-
|
|
36
|
-
\item{\code{id.metric}}{\code{(character}) Metric used as 'identity'.}
|
|
37
|
-
|
|
38
|
-
\item{\code{id.ingroup}}{\code{(logical}) Identity bins considered in-group.}
|
|
39
|
-
|
|
40
|
-
\item{\code{call}}{\code{(call)} Call producing this object.}
|
|
41
|
-
}}
|
|
42
|
-
|
|
43
|
-
\author{
|
|
44
|
-
Luis M. Rodriguez-R [aut, cre]
|
|
45
|
-
}
|
|
@@ -1,24 +0,0 @@
|
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
|
2
|
-
% Please edit documentation in R/recplot2.R
|
|
3
|
-
\name{enve.recplot2.ANIr}
|
|
4
|
-
\alias{enve.recplot2.ANIr}
|
|
5
|
-
\title{Enveomics: Recruitment Plot (2) ANI Estimate}
|
|
6
|
-
\usage{
|
|
7
|
-
enve.recplot2.ANIr(x, range = c(0, Inf))
|
|
8
|
-
}
|
|
9
|
-
\arguments{
|
|
10
|
-
\item{x}{\code{\link{enve.RecPlot2}} object.}
|
|
11
|
-
|
|
12
|
-
\item{range}{Range of identities to be considered. By default, the full range
|
|
13
|
-
is used (note that the upper boundary is \code{Inf} and not 100 because
|
|
14
|
-
recruitment plots can also be built with bit-scores). To use only
|
|
15
|
-
intra-population matches (with identities), use c(95,100). To use only
|
|
16
|
-
inter-population values, use c(0,95).}
|
|
17
|
-
}
|
|
18
|
-
\description{
|
|
19
|
-
Estimate the Average Nucleotide Identity from reads (ANIr) from a
|
|
20
|
-
recruitment plot.
|
|
21
|
-
}
|
|
22
|
-
\author{
|
|
23
|
-
Luis M. Rodriguez-R [aut, cre]
|
|
24
|
-
}
|
|
@@ -1,77 +0,0 @@
|
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
|
2
|
-
% Please edit documentation in R/recplot2.R
|
|
3
|
-
\name{enve.recplot2}
|
|
4
|
-
\alias{enve.recplot2}
|
|
5
|
-
\title{Enveomics: Recruitment Plot (2)}
|
|
6
|
-
\usage{
|
|
7
|
-
enve.recplot2(
|
|
8
|
-
prefix,
|
|
9
|
-
plot = TRUE,
|
|
10
|
-
pos.breaks = 1000,
|
|
11
|
-
pos.breaks.tsv = NA,
|
|
12
|
-
id.breaks = 60,
|
|
13
|
-
id.free.range = FALSE,
|
|
14
|
-
id.metric = c("identity", "corrected identity", "bit score"),
|
|
15
|
-
id.summary = sum,
|
|
16
|
-
id.cutoff = 95,
|
|
17
|
-
threads = 2,
|
|
18
|
-
verbose = TRUE,
|
|
19
|
-
...
|
|
20
|
-
)
|
|
21
|
-
}
|
|
22
|
-
\arguments{
|
|
23
|
-
\item{prefix}{Path to the prefix of the \code{BlastTab.catsbj.pl} output files. At
|
|
24
|
-
least the files .rec and .lim must exist with this prefix.}
|
|
25
|
-
|
|
26
|
-
\item{plot}{Should the object be plotted?}
|
|
27
|
-
|
|
28
|
-
\item{pos.breaks}{Breaks in the positions histogram. It can also be a vector of break
|
|
29
|
-
points, and values outside the range are ignored. If zero (0), it
|
|
30
|
-
uses the sequence breaks as defined in the .lim file, which means
|
|
31
|
-
one bin per contig (or gene, if the mapping is agains genes). Ignored
|
|
32
|
-
if `pos.breaks.tsv` is passed.}
|
|
33
|
-
|
|
34
|
-
\item{pos.breaks.tsv}{Path to a list of (absolute) coordinates to use as position breaks.
|
|
35
|
-
This tab-delimited file can be produced by \code{GFF.catsbj.pl}, and it
|
|
36
|
-
must contain at least one column: coordinates of the break positions of
|
|
37
|
-
each position bin. If it has a second column, this is used as the name
|
|
38
|
-
of the position bin that ends at the given coordinate (the first row is
|
|
39
|
-
ignored). Any additional columns are currently ignored. If \code{NA},
|
|
40
|
-
position bins are determined by \code{pos.breaks}.}
|
|
41
|
-
|
|
42
|
-
\item{id.breaks}{Breaks in the identity histogram. It can also be a vector of break
|
|
43
|
-
points, and values outside the range are ignored.}
|
|
44
|
-
|
|
45
|
-
\item{id.free.range}{Indicates that the range should be freely set from the observed
|
|
46
|
-
values. Otherwise, 70-100\% is included in the identity histogram
|
|
47
|
-
(default).}
|
|
48
|
-
|
|
49
|
-
\item{id.metric}{Metric of identity to be used (Y-axis). Corrected identity is only
|
|
50
|
-
supported if the original BLAST file included sequence lengths.}
|
|
51
|
-
|
|
52
|
-
\item{id.summary}{Function summarizing the identity bins. Other recommended options
|
|
53
|
-
include: \code{median} to estimate the median instead of total bins, and
|
|
54
|
-
\code{function(x) mlv(x,method='parzen')$M} to estimate the mode.}
|
|
55
|
-
|
|
56
|
-
\item{id.cutoff}{Cutoff of identity metric above which the hits are considered
|
|
57
|
-
\code{in-group}. The 95\% identity corresponds to the expectation of
|
|
58
|
-
ANI<95\% within species.}
|
|
59
|
-
|
|
60
|
-
\item{threads}{Number of threads to use.}
|
|
61
|
-
|
|
62
|
-
\item{verbose}{Indicates if the function should report the advance.}
|
|
63
|
-
|
|
64
|
-
\item{...}{Any additional parameters supported by \code{\link{plot.enve.RecPlot2}}.}
|
|
65
|
-
}
|
|
66
|
-
\value{
|
|
67
|
-
Returns an object of class \code{\link{enve.RecPlot2}}.
|
|
68
|
-
}
|
|
69
|
-
\description{
|
|
70
|
-
Produces recruitment plots provided that \code{BlastTab.catsbj.pl} has
|
|
71
|
-
been previously executed.
|
|
72
|
-
}
|
|
73
|
-
\author{
|
|
74
|
-
Luis M. Rodriguez-R [aut, cre]
|
|
75
|
-
|
|
76
|
-
Kenji Gerhardt [aut]
|
|
77
|
-
}
|
|
@@ -1,25 +0,0 @@
|
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
|
2
|
-
% Please edit documentation in R/recplot2.R
|
|
3
|
-
\name{enve.recplot2.__counts}
|
|
4
|
-
\alias{enve.recplot2.__counts}
|
|
5
|
-
\title{Enveomics: Recruitment Plot (2) Internal Ancillary Function}
|
|
6
|
-
\usage{
|
|
7
|
-
enve.recplot2.__counts(x, pos.breaks, id.breaks, rec.idcol)
|
|
8
|
-
}
|
|
9
|
-
\arguments{
|
|
10
|
-
\item{x}{\code{\link{enve.RecPlot2}} object}
|
|
11
|
-
|
|
12
|
-
\item{pos.breaks}{Position breaks}
|
|
13
|
-
|
|
14
|
-
\item{id.breaks}{Identity breaks}
|
|
15
|
-
|
|
16
|
-
\item{rec.idcol}{Identity column to use}
|
|
17
|
-
}
|
|
18
|
-
\description{
|
|
19
|
-
Internal ancillary function (see \code{\link{enve.recplot2}}).
|
|
20
|
-
}
|
|
21
|
-
\author{
|
|
22
|
-
Luis M. Rodriguez-R [aut, cre]
|
|
23
|
-
|
|
24
|
-
Kenji Gerhardt [aut]
|
|
25
|
-
}
|
|
@@ -1,21 +0,0 @@
|
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
|
2
|
-
% Please edit documentation in R/recplot2.R
|
|
3
|
-
\name{enve.recplot2.__peakHist}
|
|
4
|
-
\alias{enve.recplot2.__peakHist}
|
|
5
|
-
\title{Enveomics: Recruitment Plot (2) Peak S4 Class - Internal Ancillary Function}
|
|
6
|
-
\usage{
|
|
7
|
-
enve.recplot2.__peakHist(x, mids, counts = TRUE)
|
|
8
|
-
}
|
|
9
|
-
\arguments{
|
|
10
|
-
\item{x}{\code{\link{enve.RecPlot2.Peak}} object}
|
|
11
|
-
|
|
12
|
-
\item{mids}{Midpoints}
|
|
13
|
-
|
|
14
|
-
\item{counts}{Counts}
|
|
15
|
-
}
|
|
16
|
-
\description{
|
|
17
|
-
Internal ancillary function (see \code{\link{enve.RecPlot2.Peak}}).
|
|
18
|
-
}
|
|
19
|
-
\author{
|
|
20
|
-
Luis M. Rodriguez-R [aut, cre]
|
|
21
|
-
}
|
|
@@ -1,19 +0,0 @@
|
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
|
2
|
-
% Please edit documentation in R/recplot2.R
|
|
3
|
-
\name{enve.recplot2.__whichClosestPeak}
|
|
4
|
-
\alias{enve.recplot2.__whichClosestPeak}
|
|
5
|
-
\title{Enveomics: Recruitment Plot (2) Peak Finder - Internal Ancillary Function}
|
|
6
|
-
\usage{
|
|
7
|
-
enve.recplot2.__whichClosestPeak(peak, peaks)
|
|
8
|
-
}
|
|
9
|
-
\arguments{
|
|
10
|
-
\item{peak}{Query \code{\link{enve.RecPlot2.Peak}} object}
|
|
11
|
-
|
|
12
|
-
\item{peaks}{list of \code{\link{enve.RecPlot2.Peak}} objects}
|
|
13
|
-
}
|
|
14
|
-
\description{
|
|
15
|
-
Internal ancillary function (see \code{\link{enve.recplot2.findPeaks}}).
|
|
16
|
-
}
|
|
17
|
-
\author{
|
|
18
|
-
Luis M. Rodriguez-R [aut, cre]
|
|
19
|
-
}
|
|
@@ -1,19 +0,0 @@
|
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
|
2
|
-
% Please edit documentation in R/recplot2.R
|
|
3
|
-
\name{enve.recplot2.changeCutoff}
|
|
4
|
-
\alias{enve.recplot2.changeCutoff}
|
|
5
|
-
\title{Enveomics: Recruitment Plot (2) Change Cutoff}
|
|
6
|
-
\usage{
|
|
7
|
-
enve.recplot2.changeCutoff(rp, new.cutoff = 98)
|
|
8
|
-
}
|
|
9
|
-
\arguments{
|
|
10
|
-
\item{rp}{\code{\link{enve.RecPlot2}} object.}
|
|
11
|
-
|
|
12
|
-
\item{new.cutoff}{New cutoff to use.}
|
|
13
|
-
}
|
|
14
|
-
\description{
|
|
15
|
-
Change the intra-species cutoff of an existing recruitment plot.
|
|
16
|
-
}
|
|
17
|
-
\author{
|
|
18
|
-
Luis M. Rodriguez-R [aut, cre]
|
|
19
|
-
}
|
|
@@ -1,47 +0,0 @@
|
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
|
2
|
-
% Please edit documentation in R/recplot2.R
|
|
3
|
-
\name{enve.recplot2.compareIdentities}
|
|
4
|
-
\alias{enve.recplot2.compareIdentities}
|
|
5
|
-
\title{Enveomics: Recruitment Plot (2) Compare Identities}
|
|
6
|
-
\usage{
|
|
7
|
-
enve.recplot2.compareIdentities(
|
|
8
|
-
x,
|
|
9
|
-
y,
|
|
10
|
-
method = "hellinger",
|
|
11
|
-
smooth.par = NULL,
|
|
12
|
-
pseudocounts = 0,
|
|
13
|
-
max.deviation = 0.75
|
|
14
|
-
)
|
|
15
|
-
}
|
|
16
|
-
\arguments{
|
|
17
|
-
\item{x}{First \code{\link{enve.RecPlot2}} object.}
|
|
18
|
-
|
|
19
|
-
\item{y}{Second \code{\link{enve.RecPlot2}} object.}
|
|
20
|
-
|
|
21
|
-
\item{method}{Distance method to use. This should be (an unambiguous abbreviation of)
|
|
22
|
-
one of:
|
|
23
|
-
\itemize{
|
|
24
|
-
\item{"hellinger" (\emph{Hellinger, 1090, doi:10.1515/crll.1909.136.210}),}
|
|
25
|
-
\item{"bhattacharyya" (\emph{Bhattacharyya, 1943, Bull. Calcutta Math. Soc. 35}),}
|
|
26
|
-
\item{"kl" or "kullback-leibler" (\emph{Kullback & Leibler, 1951,
|
|
27
|
-
doi:10.1214/aoms/1177729694}), or}
|
|
28
|
-
\item{"euclidean"}
|
|
29
|
-
}}
|
|
30
|
-
|
|
31
|
-
\item{smooth.par}{Smoothing parameter for cubic spline smoothing. Use 0 for no smoothing.
|
|
32
|
-
Use \code{NULL} to automatically determine this value using leave-one-out
|
|
33
|
-
cross-validation (see \code{smooth.spline} parameter \code{spar}).}
|
|
34
|
-
|
|
35
|
-
\item{pseudocounts}{Smoothing parameter for Laplace smoothing. Use 0 for no smoothing, or
|
|
36
|
-
1 for add-one smoothing.}
|
|
37
|
-
|
|
38
|
-
\item{max.deviation}{Maximum mean deviation between identity breaks tolerated (as percent
|
|
39
|
-
identity). Difference in number of \code{id.breaks} is never tolerated.}
|
|
40
|
-
}
|
|
41
|
-
\description{
|
|
42
|
-
Compare the distribution of identities between two
|
|
43
|
-
\code{\link{enve.RecPlot2}} objects.
|
|
44
|
-
}
|
|
45
|
-
\author{
|
|
46
|
-
Luis M. Rodriguez-R [aut, cre]
|
|
47
|
-
}
|
|
@@ -1,29 +0,0 @@
|
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
|
2
|
-
% Please edit documentation in R/recplot2.R
|
|
3
|
-
\name{enve.recplot2.coordinates}
|
|
4
|
-
\alias{enve.recplot2.coordinates}
|
|
5
|
-
\title{Enveomics: Recruitment Plot (2) Coordinates}
|
|
6
|
-
\usage{
|
|
7
|
-
enve.recplot2.coordinates(x, bins)
|
|
8
|
-
}
|
|
9
|
-
\arguments{
|
|
10
|
-
\item{x}{\code{\link{enve.RecPlot2}} object.}
|
|
11
|
-
|
|
12
|
-
\item{bins}{Vector of selected bins to return. It can be a vector of logical values
|
|
13
|
-
with the same length as \code{x$pos.breaks-1} or a vector of integers. If
|
|
14
|
-
missing, returns the coordinates of all windows.}
|
|
15
|
-
}
|
|
16
|
-
\value{
|
|
17
|
-
Returns a data.frame with five columns: \code{name.from} (character),
|
|
18
|
-
\code{pos.from} (numeric), \code{name.to} (character), \code{pos.to}
|
|
19
|
-
(numeric), and \code{seq.name} (character).
|
|
20
|
-
The first two correspond to sequence and position of the start point of the
|
|
21
|
-
bin. The next two correspond to the sequence and position of the end point of
|
|
22
|
-
the bin. The last one indicates the name of the sequence (if defined).
|
|
23
|
-
}
|
|
24
|
-
\description{
|
|
25
|
-
Returns the sequence name and coordinates of the requested position bins.
|
|
26
|
-
}
|
|
27
|
-
\author{
|
|
28
|
-
Luis M. Rodriguez-R [aut, cre]
|
|
29
|
-
}
|
|
@@ -1,18 +0,0 @@
|
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
|
2
|
-
% Please edit documentation in R/recplot2.R
|
|
3
|
-
\name{enve.recplot2.corePeak}
|
|
4
|
-
\alias{enve.recplot2.corePeak}
|
|
5
|
-
\title{Enveomics: Recruitment Plot (2) Core Peak Finder}
|
|
6
|
-
\usage{
|
|
7
|
-
enve.recplot2.corePeak(x)
|
|
8
|
-
}
|
|
9
|
-
\arguments{
|
|
10
|
-
\item{x}{\code{list} of \code{\link{enve.RecPlot2.Peak}} objects.}
|
|
11
|
-
}
|
|
12
|
-
\description{
|
|
13
|
-
Finds the peak in a list of peaks that is most likely to represent the
|
|
14
|
-
"core genome" of a population.
|
|
15
|
-
}
|
|
16
|
-
\author{
|
|
17
|
-
Luis M. Rodriguez-R [aut, cre]
|
|
18
|
-
}
|