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,45 +0,0 @@
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
2
|
-
% Please edit documentation in R/recplot2.R
|
3
|
-
\name{enve.recplot2.extractWindows}
|
4
|
-
\alias{enve.recplot2.extractWindows}
|
5
|
-
\title{Enveomics: Recruitment Plot (2) Extract Windows}
|
6
|
-
\usage{
|
7
|
-
enve.recplot2.extractWindows(
|
8
|
-
rp,
|
9
|
-
peak,
|
10
|
-
lower.tail = TRUE,
|
11
|
-
significance = 0.05,
|
12
|
-
seq.names = FALSE
|
13
|
-
)
|
14
|
-
}
|
15
|
-
\arguments{
|
16
|
-
\item{rp}{Recruitment plot, a \code{\link{enve.RecPlot2}} object.}
|
17
|
-
|
18
|
-
\item{peak}{Peak, an \code{\link{enve.RecPlot2.Peak}} object. If list, it is assumed to be a
|
19
|
-
list of \code{\link{enve.RecPlot2.Peak}} objects, in which case the core peak is
|
20
|
-
used (see \code{\link{enve.recplot2.corePeak}}).}
|
21
|
-
|
22
|
-
\item{lower.tail}{If \code{FALSE}, it returns windows significantly above the peak in
|
23
|
-
sequencing depth.}
|
24
|
-
|
25
|
-
\item{significance}{Significance threshold (alpha) to select windows.}
|
26
|
-
|
27
|
-
\item{seq.names}{Returns subject sequence names instead of a vector of Booleans. If
|
28
|
-
the recruitment plot was generated with named position bins (e.g, using
|
29
|
-
\code{pos.breaks=0} or a two-column \code{pos.breaks.tsv}), it returns a
|
30
|
-
vector of characters (the sequence identifiers), otherwise it returns a
|
31
|
-
data.frame with a name column and two columns of coordinates.}
|
32
|
-
}
|
33
|
-
\value{
|
34
|
-
Returns a vector of logicals if \code{seq.names = FALSE}.
|
35
|
-
If \code{seq.names = TRUE}, it returns a data.frame with five columns:
|
36
|
-
\code{name.from}, \code{name.to}, \code{pos.from}, \code{pos.to}, and
|
37
|
-
\code{seq.name} (see \code{\link{enve.recplot2.coordinates}}).
|
38
|
-
}
|
39
|
-
\description{
|
40
|
-
Extract windows significantly below (or above) the peak in sequencing
|
41
|
-
depth.
|
42
|
-
}
|
43
|
-
\author{
|
44
|
-
Luis M. Rodriguez-R [aut, cre]
|
45
|
-
}
|
@@ -1,36 +0,0 @@
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
2
|
-
% Please edit documentation in R/recplot2.R
|
3
|
-
\name{enve.recplot2.findPeaks}
|
4
|
-
\alias{enve.recplot2.findPeaks}
|
5
|
-
\title{Enveomics: Recruitment Plot (2) Peak Finder}
|
6
|
-
\usage{
|
7
|
-
enve.recplot2.findPeaks(x, method = "emauto", ...)
|
8
|
-
}
|
9
|
-
\arguments{
|
10
|
-
\item{x}{An \code{\link{enve.RecPlot2}} object.}
|
11
|
-
|
12
|
-
\item{method}{Peak-finder method. This should be one of:
|
13
|
-
\itemize{
|
14
|
-
\item \strong{emauto}
|
15
|
-
(Expectation-Maximization with auto-selection of components)
|
16
|
-
\item \strong{em}
|
17
|
-
(Expectation-Maximization)
|
18
|
-
\item \strong{mower}
|
19
|
-
(Custom distribution-mowing method)
|
20
|
-
}}
|
21
|
-
|
22
|
-
\item{...}{Any additional parameters supported by
|
23
|
-
\code{\link{enve.recplot2.findPeaks}}.}
|
24
|
-
}
|
25
|
-
\value{
|
26
|
-
Returns a list of \code{\link{enve.RecPlot2.Peak}} objects.
|
27
|
-
}
|
28
|
-
\description{
|
29
|
-
Identifies peaks in the population histogram potentially indicating
|
30
|
-
sub-population mixtures.
|
31
|
-
}
|
32
|
-
\author{
|
33
|
-
Luis M. Rodriguez-R [aut, cre]
|
34
|
-
|
35
|
-
export
|
36
|
-
}
|
@@ -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.findPeaks.__em_e}
|
4
|
-
\alias{enve.recplot2.findPeaks.__em_e}
|
5
|
-
\title{Enveomics: Recruitment Plot (2) EM Peak Finder - Internal Ancillary Function Expectation}
|
6
|
-
\usage{
|
7
|
-
enve.recplot2.findPeaks.__em_e(x, theta)
|
8
|
-
}
|
9
|
-
\arguments{
|
10
|
-
\item{x}{Vector of log-transformed sequencing depths}
|
11
|
-
|
12
|
-
\item{theta}{Parameters list}
|
13
|
-
}
|
14
|
-
\description{
|
15
|
-
Internal ancillary function (see \code{\link{enve.recplot2.findPeaks.em}}).
|
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.findPeaks.__em_m}
|
4
|
-
\alias{enve.recplot2.findPeaks.__em_m}
|
5
|
-
\title{Enveomics: Recruitment Plot (2) Em Peak Finder - Internal Ancillary Function Maximization}
|
6
|
-
\usage{
|
7
|
-
enve.recplot2.findPeaks.__em_m(x, posterior)
|
8
|
-
}
|
9
|
-
\arguments{
|
10
|
-
\item{x}{Vector of log-transformed sequencing depths}
|
11
|
-
|
12
|
-
\item{posterior}{Posterior probability}
|
13
|
-
}
|
14
|
-
\description{
|
15
|
-
Internal ancillary function (see \code{\link{enve.recplot2.findPeaks.em}}).
|
16
|
-
}
|
17
|
-
\author{
|
18
|
-
Luis M. Rodriguez-R [aut, cre]
|
19
|
-
}
|
@@ -1,27 +0,0 @@
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
2
|
-
% Please edit documentation in R/recplot2.R
|
3
|
-
\name{enve.recplot2.findPeaks.__emauto_one}
|
4
|
-
\alias{enve.recplot2.findPeaks.__emauto_one}
|
5
|
-
\title{Enveomics: Recruitment Plot (2) EMauto Peak Finder - Internal Ancillary Function}
|
6
|
-
\usage{
|
7
|
-
enve.recplot2.findPeaks.__emauto_one(x, comp, do_crit, best, verbose, ...)
|
8
|
-
}
|
9
|
-
\arguments{
|
10
|
-
\item{x}{\code{\link{enve.RecPlot2}} object}
|
11
|
-
|
12
|
-
\item{comp}{Components}
|
13
|
-
|
14
|
-
\item{do_crit}{Function estimating the criterion}
|
15
|
-
|
16
|
-
\item{best}{Best solution thus far}
|
17
|
-
|
18
|
-
\item{verbose}{If verbose}
|
19
|
-
|
20
|
-
\item{...}{Additional parameters for \code{\link{enve.recplot2.findPeaks.em}}}
|
21
|
-
}
|
22
|
-
\description{
|
23
|
-
Internal ancillary function (see \code{\link{enve.recplot2.findPeaks.emauto}}).
|
24
|
-
}
|
25
|
-
\author{
|
26
|
-
Luis M. Rodriguez-R [aut, cre]
|
27
|
-
}
|
@@ -1,52 +0,0 @@
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
2
|
-
% Please edit documentation in R/recplot2.R
|
3
|
-
\name{enve.recplot2.findPeaks.__mow_one}
|
4
|
-
\alias{enve.recplot2.findPeaks.__mow_one}
|
5
|
-
\title{Enveomics: Recruitment Plot (2) Mowing Peak Finder - Internal Ancillary Function 1}
|
6
|
-
\usage{
|
7
|
-
enve.recplot2.findPeaks.__mow_one(
|
8
|
-
lsd1,
|
9
|
-
min.points,
|
10
|
-
quant.est,
|
11
|
-
mlv.opts,
|
12
|
-
fitdist.opts,
|
13
|
-
with.skewness,
|
14
|
-
optim.rounds,
|
15
|
-
optim.epsilon,
|
16
|
-
n.total,
|
17
|
-
merge.logdist,
|
18
|
-
verbose,
|
19
|
-
log
|
20
|
-
)
|
21
|
-
}
|
22
|
-
\arguments{
|
23
|
-
\item{lsd1}{Vector of log-transformed sequencing depths}
|
24
|
-
|
25
|
-
\item{min.points}{Minimum number of points}
|
26
|
-
|
27
|
-
\item{quant.est}{Quantile estimate}
|
28
|
-
|
29
|
-
\item{mlv.opts}{List of options for \code{mlv}}
|
30
|
-
|
31
|
-
\item{fitdist.opts}{List of options for \code{fitdist}}
|
32
|
-
|
33
|
-
\item{with.skewness}{If skewed-normal should be used}
|
34
|
-
|
35
|
-
\item{optim.rounds}{Maximum number of optimization rounds}
|
36
|
-
|
37
|
-
\item{optim.epsilon}{Minimum difference considered negligible}
|
38
|
-
|
39
|
-
\item{n.total}{Global number of windows}
|
40
|
-
|
41
|
-
\item{merge.logdist}{Attempted \code{merge.logdist} parameter}
|
42
|
-
|
43
|
-
\item{verbose}{If verbose}
|
44
|
-
|
45
|
-
\item{log}{If log-transformed depths}
|
46
|
-
}
|
47
|
-
\description{
|
48
|
-
Internall ancillary function (see \code{\link{enve.recplot2.findPeaks.mower}}).
|
49
|
-
}
|
50
|
-
\author{
|
51
|
-
Luis M. Rodriguez-R [aut, cre]
|
52
|
-
}
|
@@ -1,17 +0,0 @@
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
2
|
-
% Please edit documentation in R/recplot2.R
|
3
|
-
\name{enve.recplot2.findPeaks.__mower}
|
4
|
-
\alias{enve.recplot2.findPeaks.__mower}
|
5
|
-
\title{Enveomics: Recruitment Plot (2) Mowing Peak Finder - Internal Ancillary Function 2}
|
6
|
-
\usage{
|
7
|
-
enve.recplot2.findPeaks.__mower(peaks.opts)
|
8
|
-
}
|
9
|
-
\arguments{
|
10
|
-
\item{peaks.opts}{List of options for \code{\link{enve.recplot2.findPeaks.__mow_one}}}
|
11
|
-
}
|
12
|
-
\description{
|
13
|
-
Internal ancillary function (see \code{\link{enve.recplot2.findPeaks.mower}}).
|
14
|
-
}
|
15
|
-
\author{
|
16
|
-
Luis M. Rodriguez-R [aut, cre]
|
17
|
-
}
|
@@ -1,51 +0,0 @@
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
2
|
-
% Please edit documentation in R/recplot2.R
|
3
|
-
\name{enve.recplot2.findPeaks.em}
|
4
|
-
\alias{enve.recplot2.findPeaks.em}
|
5
|
-
\title{Enveomics: Recruitment Plot (2) Em Peak Finder}
|
6
|
-
\usage{
|
7
|
-
enve.recplot2.findPeaks.em(
|
8
|
-
x,
|
9
|
-
max.iter = 1000,
|
10
|
-
ll.diff.res = 1e-08,
|
11
|
-
components = 2,
|
12
|
-
rm.top = 0.05,
|
13
|
-
verbose = FALSE,
|
14
|
-
init,
|
15
|
-
log = TRUE
|
16
|
-
)
|
17
|
-
}
|
18
|
-
\arguments{
|
19
|
-
\item{x}{An \code{\link{enve.RecPlot2}} object.}
|
20
|
-
|
21
|
-
\item{max.iter}{Maximum number of EM iterations.}
|
22
|
-
|
23
|
-
\item{ll.diff.res}{Maximum Log-Likelihood difference to be considered as convergent.}
|
24
|
-
|
25
|
-
\item{components}{Number of distributions assumed in the mixture.}
|
26
|
-
|
27
|
-
\item{rm.top}{Top-values to remove before finding peaks, as a quantile probability.
|
28
|
-
This step is useful to remove highly conserved regions, but can be
|
29
|
-
turned off by setting \code{rm.top=0}. The quantile is determined
|
30
|
-
\strong{after} removing zero-coverage windows.}
|
31
|
-
|
32
|
-
\item{verbose}{Display (mostly debugging) information.}
|
33
|
-
|
34
|
-
\item{init}{Initialization parameters. By default, these are derived from k-means
|
35
|
-
clustering. A named list with vectors for \code{mu}, \code{sd}, and
|
36
|
-
\code{alpha}, each of length \code{components}.}
|
37
|
-
|
38
|
-
\item{log}{Logical value indicating if the estimations should be performed in
|
39
|
-
natural logarithm units. Do not change unless you know what you're
|
40
|
-
doing.}
|
41
|
-
}
|
42
|
-
\value{
|
43
|
-
Returns a list of \code{\link{enve.RecPlot2.Peak}} objects.
|
44
|
-
}
|
45
|
-
\description{
|
46
|
-
Identifies peaks in the population histogram using a Gaussian Mixture
|
47
|
-
Model Expectation Maximization (GMM-EM) method.
|
48
|
-
}
|
49
|
-
\author{
|
50
|
-
Luis M. Rodriguez-R [aut, cre]
|
51
|
-
}
|
@@ -1,43 +0,0 @@
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
2
|
-
% Please edit documentation in R/recplot2.R
|
3
|
-
\name{enve.recplot2.findPeaks.emauto}
|
4
|
-
\alias{enve.recplot2.findPeaks.emauto}
|
5
|
-
\title{Enveomics: Recruitment Plot (2) Emauto Peak Finder}
|
6
|
-
\usage{
|
7
|
-
enve.recplot2.findPeaks.emauto(
|
8
|
-
x,
|
9
|
-
components = seq(1, 5),
|
10
|
-
criterion = "aic",
|
11
|
-
merge.tol = 2L,
|
12
|
-
verbose = FALSE,
|
13
|
-
...
|
14
|
-
)
|
15
|
-
}
|
16
|
-
\arguments{
|
17
|
-
\item{x}{An \code{\link{enve.RecPlot2}} object.}
|
18
|
-
|
19
|
-
\item{components}{A vector of number of components to evaluate.}
|
20
|
-
|
21
|
-
\item{criterion}{Criterion to use for components selection. Must be one of:
|
22
|
-
\code{aic} (Akaike Information Criterion), \code{bic} or \code{sbc}
|
23
|
-
(Bayesian Information Criterion or Schwarz Criterion).}
|
24
|
-
|
25
|
-
\item{merge.tol}{When attempting to merge peaks with very similar sequencing depth, use
|
26
|
-
this number of significant digits (in log-scale).}
|
27
|
-
|
28
|
-
\item{verbose}{Display (mostly debugging) information.}
|
29
|
-
|
30
|
-
\item{...}{Any additional parameters supported by
|
31
|
-
\code{\link{enve.recplot2.findPeaks.em}}.}
|
32
|
-
}
|
33
|
-
\value{
|
34
|
-
Returns a list of \code{\link{enve.RecPlot2.Peak}} objects.
|
35
|
-
}
|
36
|
-
\description{
|
37
|
-
Identifies peaks in the population histogram using a Gaussian Mixture
|
38
|
-
Model Expectation Maximization (GMM-EM) method with number of components
|
39
|
-
automatically detected.
|
40
|
-
}
|
41
|
-
\author{
|
42
|
-
Luis M. Rodriguez-R [aut, cre]
|
43
|
-
}
|
@@ -1,82 +0,0 @@
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
2
|
-
% Please edit documentation in R/recplot2.R
|
3
|
-
\name{enve.recplot2.findPeaks.mower}
|
4
|
-
\alias{enve.recplot2.findPeaks.mower}
|
5
|
-
\title{Enveomics: Recruitment Plot (2) Mowing Peak Finder}
|
6
|
-
\usage{
|
7
|
-
enve.recplot2.findPeaks.mower(
|
8
|
-
x,
|
9
|
-
min.points = 10,
|
10
|
-
quant.est = c(0.002, 0.998),
|
11
|
-
mlv.opts = list(method = "parzen"),
|
12
|
-
fitdist.opts.sn = list(distr = "sn", method = "qme", probs = c(0.1, 0.5, 0.8), start =
|
13
|
-
list(omega = 1, alpha = -1), lower = c(0, -Inf, -Inf)),
|
14
|
-
fitdist.opts.norm = list(distr = "norm", method = "qme", probs = c(0.4, 0.6), start =
|
15
|
-
list(sd = 1), lower = c(0, -Inf)),
|
16
|
-
rm.top = 0.05,
|
17
|
-
with.skewness = TRUE,
|
18
|
-
optim.rounds = 200,
|
19
|
-
optim.epsilon = 1e-04,
|
20
|
-
merge.logdist = log(1.75),
|
21
|
-
verbose = FALSE,
|
22
|
-
log = TRUE
|
23
|
-
)
|
24
|
-
}
|
25
|
-
\arguments{
|
26
|
-
\item{x}{An \code{\link{enve.RecPlot2}} object.}
|
27
|
-
|
28
|
-
\item{min.points}{Minimum number of points in the quantile-estimation-range
|
29
|
-
\code{(quant.est)} to estimate a peak.}
|
30
|
-
|
31
|
-
\item{quant.est}{Range of quantiles to be used in the estimation of a peak's
|
32
|
-
parameters.}
|
33
|
-
|
34
|
-
\item{mlv.opts}{Ignored. For backwards compatibility.}
|
35
|
-
|
36
|
-
\item{fitdist.opts.sn}{Options passed to \code{fitdist} to estimate the standard deviation if
|
37
|
-
\code{with.skewness=TRUE}. Note that the \code{start} parameter will be
|
38
|
-
ammended with \code{xi=estimated} mode for each peak.}
|
39
|
-
|
40
|
-
\item{fitdist.opts.norm}{Options passed to \code{fitdist} to estimate the standard deviation if
|
41
|
-
\code{with.skewness=FALSE}. Note that the \code{start} parameter will be
|
42
|
-
ammended with \code{mean=estimated} mode for each peak.}
|
43
|
-
|
44
|
-
\item{rm.top}{Top-values to remove before finding peaks, as a quantile probability.
|
45
|
-
This step is useful to remove highly conserved regions, but can be
|
46
|
-
turned off by setting \code{rm.top=0}. The quantile is determined
|
47
|
-
\strong{after} removing zero-coverage windows.}
|
48
|
-
|
49
|
-
\item{with.skewness}{Allow skewness correction of the peaks. Typically, the
|
50
|
-
sequencing-depth distribution for a single peak is left-skewed, due
|
51
|
-
partly (but not exclusively) to fragmentation and mapping sensitivity.
|
52
|
-
See \emph{Lindner et al 2013, Bioinformatics 29(10):1260-7} for an
|
53
|
-
alternative solution for the first problem (fragmentation) called
|
54
|
-
"tail distribution".}
|
55
|
-
|
56
|
-
\item{optim.rounds}{Maximum rounds of peak optimization.}
|
57
|
-
|
58
|
-
\item{optim.epsilon}{Trace change at which optimization stops (unless \code{optim.rounds} is
|
59
|
-
reached first). The trace change is estimated as the sum of square
|
60
|
-
differences between parameters in one round and those from two rounds
|
61
|
-
earlier (to avoid infinite loops from approximation).}
|
62
|
-
|
63
|
-
\item{merge.logdist}{Maximum value of \code{|log-ratio|} between centrality parameters in peaks
|
64
|
-
to attempt merging. The default of ~0.22 corresponds to a maximum
|
65
|
-
difference of 25\%.}
|
66
|
-
|
67
|
-
\item{verbose}{Display (mostly debugging) information.}
|
68
|
-
|
69
|
-
\item{log}{Logical value indicating if the estimations should be performed in
|
70
|
-
natural logarithm units. Do not change unless you know what you're
|
71
|
-
doing.}
|
72
|
-
}
|
73
|
-
\value{
|
74
|
-
Returns a list of \code{\link{enve.RecPlot2.Peak}} objects.
|
75
|
-
}
|
76
|
-
\description{
|
77
|
-
Identifies peaks in the population histogram potentially indicating
|
78
|
-
sub-population mixtures, using a custom distribution-mowing method.
|
79
|
-
}
|
80
|
-
\author{
|
81
|
-
Luis M. Rodriguez-R [aut, cre]
|
82
|
-
}
|
@@ -1,59 +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.Peak-class}
|
5
|
-
\alias{enve.RecPlot2.Peak-class}
|
6
|
-
\alias{enve.RecPlot2.Peak}
|
7
|
-
\title{Enveomics: Recruitment Plot (2) Peak - S4 Class}
|
8
|
-
\description{
|
9
|
-
Enve-omics representation of a peak in the sequencing depth histogram
|
10
|
-
of a Recruitment plot (see \code{\link{enve.recplot2.findPeaks}}).
|
11
|
-
}
|
12
|
-
\section{Slots}{
|
13
|
-
|
14
|
-
\describe{
|
15
|
-
\item{\code{dist}}{\code{(character)}
|
16
|
-
Distribution of the peak. Currently supported: \code{norm} (normal) and \code{sn}
|
17
|
-
(skew-normal).}
|
18
|
-
|
19
|
-
\item{\code{values}}{\code{(numeric)}
|
20
|
-
Sequencing depth values predicted to conform the peak.}
|
21
|
-
|
22
|
-
\item{\code{values.res}}{\code{(numeric)}
|
23
|
-
Sequencing depth values not explained by this or previously identified
|
24
|
-
peaks.}
|
25
|
-
|
26
|
-
\item{\code{mode}}{\code{(numeric)}
|
27
|
-
Seed-value of mode anchoring the peak.}
|
28
|
-
|
29
|
-
\item{\code{param.hat}}{\code{(list)}
|
30
|
-
Parameters of the distribution. A list of two values if dist=\code{norm} (sd
|
31
|
-
and mean), or three values if dist=\code{sn}(omega=scale, alpha=shape, and
|
32
|
-
xi=location). Note that the "dispersion" parameter is always first and
|
33
|
-
the "location" parameter is always last.}
|
34
|
-
|
35
|
-
\item{\code{n.hat}}{\code{(numeric)}
|
36
|
-
Number of bins estimated to be explained by this peak. This should
|
37
|
-
ideally be equal to the length of \code{values}, but it's not an integer.}
|
38
|
-
|
39
|
-
\item{\code{n.total}}{\code{(numeric)}
|
40
|
-
Total number of bins from which the peak was extracted. I.e., total
|
41
|
-
number of position bins with non-zero sequencing depth in the recruitment
|
42
|
-
plot (regardless of peak count).}
|
43
|
-
|
44
|
-
\item{\code{err.res}}{\code{(numeric)}
|
45
|
-
Error left after adding the peak (mower) or log-likelihood (em or emauto).}
|
46
|
-
|
47
|
-
\item{\code{merge.logdist}}{\code{(numeric)}
|
48
|
-
Attempted \code{merge.logdist} parameter.}
|
49
|
-
|
50
|
-
\item{\code{seq.depth}}{\code{(numeric)}
|
51
|
-
Best estimate available for the sequencing depth of the peak (centrality).}
|
52
|
-
|
53
|
-
\item{\code{log}}{\code{(logical)}
|
54
|
-
Indicates if the estimation was performed in natural logarithm space.}
|
55
|
-
}}
|
56
|
-
|
57
|
-
\author{
|
58
|
-
Luis M. Rodriguez-R [aut, cre]
|
59
|
-
}
|
@@ -1,27 +0,0 @@
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
2
|
-
% Please edit documentation in R/recplot2.R
|
3
|
-
\name{enve.recplot2.seqdepth}
|
4
|
-
\alias{enve.recplot2.seqdepth}
|
5
|
-
\title{Enveomics: Recruitment Plot (2) Sequencing Depth}
|
6
|
-
\usage{
|
7
|
-
enve.recplot2.seqdepth(x, sel, low.identity = FALSE)
|
8
|
-
}
|
9
|
-
\arguments{
|
10
|
-
\item{x}{\code{\link{enve.RecPlot2}} object.}
|
11
|
-
|
12
|
-
\item{sel}{Window(s) for which the sequencing depth is to be calculated. If not
|
13
|
-
passed, it returns the sequencing depth of all windows.}
|
14
|
-
|
15
|
-
\item{low.identity}{A logical indicating if the sequencing depth is to be estimated only
|
16
|
-
with low-identity matches. By default, only high-identity matches are
|
17
|
-
used.}
|
18
|
-
}
|
19
|
-
\value{
|
20
|
-
Returns a numeric vector of sequencing depths (in bp/bp).
|
21
|
-
}
|
22
|
-
\description{
|
23
|
-
Calculate the sequencing depth of the given window(s).
|
24
|
-
}
|
25
|
-
\author{
|
26
|
-
Luis M. Rodriguez-R [aut, cre]
|
27
|
-
}
|
@@ -1,36 +0,0 @@
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
2
|
-
% Please edit documentation in R/recplot2.R
|
3
|
-
\name{enve.recplot2.windowDepthThreshold}
|
4
|
-
\alias{enve.recplot2.windowDepthThreshold}
|
5
|
-
\title{Enveomics: Recruitment Plot (2) Window Depth Threshold}
|
6
|
-
\usage{
|
7
|
-
enve.recplot2.windowDepthThreshold(
|
8
|
-
rp,
|
9
|
-
peak,
|
10
|
-
lower.tail = TRUE,
|
11
|
-
significance = 0.05
|
12
|
-
)
|
13
|
-
}
|
14
|
-
\arguments{
|
15
|
-
\item{rp}{Recruitment plot, an \code{\link{enve.RecPlot2}} object.}
|
16
|
-
|
17
|
-
\item{peak}{Peak, an \code{\link{enve.RecPlot2.Peak}} object. If list, it is assumed to be a
|
18
|
-
list of \code{\link{enve.RecPlot2.Peak}} objects, in which case the core peak is
|
19
|
-
used (see \code{\link{enve.recplot2.corePeak}}).}
|
20
|
-
|
21
|
-
\item{lower.tail}{If \code{FALSE}, it returns windows significantly above the peak in
|
22
|
-
sequencing depth.}
|
23
|
-
|
24
|
-
\item{significance}{Significance threshold (alpha) to select windows.}
|
25
|
-
}
|
26
|
-
\value{
|
27
|
-
Returns a float. The units are depth if the peaks were estimated in
|
28
|
-
linear scale, or log-depth otherwise (\code{peak$log}).
|
29
|
-
}
|
30
|
-
\description{
|
31
|
-
Identifies the threshold below which windows should be identified as
|
32
|
-
variable or absent.
|
33
|
-
}
|
34
|
-
\author{
|
35
|
-
Luis M. Rodriguez-R [aut, cre]
|
36
|
-
}
|
@@ -1,23 +0,0 @@
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
2
|
-
% Please edit documentation in R/utils.R
|
3
|
-
\name{enve.selvector}
|
4
|
-
\alias{enve.selvector}
|
5
|
-
\title{Enveomics: Selection vector}
|
6
|
-
\usage{
|
7
|
-
enve.selvector(sel, dim.names)
|
8
|
-
}
|
9
|
-
\arguments{
|
10
|
-
\item{sel}{A vector of numbers, characters, or booleans.}
|
11
|
-
|
12
|
-
\item{dim.names}{A vector of names from which to select.}
|
13
|
-
}
|
14
|
-
\value{
|
15
|
-
Returns a logical vector with the same length as \code{dim.name}.
|
16
|
-
}
|
17
|
-
\description{
|
18
|
-
Normalizes a selection vector \code{sel} to a logical vector with indexes
|
19
|
-
from \code{dim.names}.
|
20
|
-
}
|
21
|
-
\author{
|
22
|
-
Luis M. Rodriguez-R [aut, cre]
|
23
|
-
}
|
@@ -1,68 +0,0 @@
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
2
|
-
% Please edit documentation in R/tribs.R
|
3
|
-
\name{enve.tribs}
|
4
|
-
\alias{enve.tribs}
|
5
|
-
\title{Enveomics: TRIBS}
|
6
|
-
\usage{
|
7
|
-
enve.tribs(
|
8
|
-
dist,
|
9
|
-
selection = labels(dist),
|
10
|
-
replicates = 1000,
|
11
|
-
summary.fx = median,
|
12
|
-
dist.method = "euclidean",
|
13
|
-
subsamples = seq(0, 1, by = 0.01),
|
14
|
-
dimensions = ceiling(length(selection) * 0.05),
|
15
|
-
metaMDS.opts = list(),
|
16
|
-
threads = 2,
|
17
|
-
verbosity = 1,
|
18
|
-
points,
|
19
|
-
pre.tribs
|
20
|
-
)
|
21
|
-
}
|
22
|
-
\arguments{
|
23
|
-
\item{dist}{Distances as a \code{dist} object.}
|
24
|
-
|
25
|
-
\item{selection}{Objects to include in the subsample. By default, all objects are
|
26
|
-
selected.}
|
27
|
-
|
28
|
-
\item{replicates}{Number of replications per point.}
|
29
|
-
|
30
|
-
\item{summary.fx}{Function to summarize the distance distributions in a given replicate. By
|
31
|
-
default, the median distance is estimated.}
|
32
|
-
|
33
|
-
\item{dist.method}{Distance method between random points and samples in the transformed
|
34
|
-
space. See \code{dist}.}
|
35
|
-
|
36
|
-
\item{subsamples}{Subsampling fractions.}
|
37
|
-
|
38
|
-
\item{dimensions}{Dimensions to use in the NMDS. By default, 5\% of the selection length.}
|
39
|
-
|
40
|
-
\item{metaMDS.opts}{Any additional options to pass to metaMDS, as \code{list}.}
|
41
|
-
|
42
|
-
\item{threads}{Number of threads to use.}
|
43
|
-
|
44
|
-
\item{verbosity}{Verbosity. Use 0 to run quietly, increase for additional information.}
|
45
|
-
|
46
|
-
\item{points}{Optional. If passed, the MDS step is skipped and this object is used
|
47
|
-
instead. It can be the \code{$points} slot of class \code{metaMDS}
|
48
|
-
(from \code{vegan}).
|
49
|
-
It must be a matrix or matrix-coercible object, with samples as rows and
|
50
|
-
dimensions as columns.}
|
51
|
-
|
52
|
-
\item{pre.tribs}{Optional. If passed, the points are recovered from this object (except if
|
53
|
-
\code{points} is also passed. This should be an \code{\link{enve.TRIBS}} object
|
54
|
-
estimated on the same objects (the selection is unimportant).}
|
55
|
-
}
|
56
|
-
\value{
|
57
|
-
Returns an \code{\link{enve.TRIBS}} object.
|
58
|
-
}
|
59
|
-
\description{
|
60
|
-
Subsample any objects in "distance space" to reduce the effect of
|
61
|
-
sample-clustering. This function was originally designed to subsample
|
62
|
-
genomes in "phylogenetic distance space", a clear case of strong
|
63
|
-
clustering bias in sampling, by Luis M. Rodriguez-R and Michael R
|
64
|
-
Weigand.
|
65
|
-
}
|
66
|
-
\author{
|
67
|
-
Luis M. Rodriguez-R [aut, cre]
|
68
|
-
}
|
@@ -1,28 +0,0 @@
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
2
|
-
% Please edit documentation in R/tribs.R
|
3
|
-
\name{enve.tribs.test}
|
4
|
-
\alias{enve.tribs.test}
|
5
|
-
\title{Enveomics: TRIBS Test}
|
6
|
-
\usage{
|
7
|
-
enve.tribs.test(dist, selection, bins = 50, ...)
|
8
|
-
}
|
9
|
-
\arguments{
|
10
|
-
\item{dist}{Distances as \code{dist} object.}
|
11
|
-
|
12
|
-
\item{selection}{Selection defining the subset.}
|
13
|
-
|
14
|
-
\item{bins}{Number of bins to evaluate in the range of distances.}
|
15
|
-
|
16
|
-
\item{...}{Any other parameters supported by \code{\link{enve.tribs}},
|
17
|
-
except \code{subsamples}.}
|
18
|
-
}
|
19
|
-
\value{
|
20
|
-
Returns an \code{\link{enve.TRIBStest}} object.
|
21
|
-
}
|
22
|
-
\description{
|
23
|
-
Estimates the empirical difference between all the distances in a set of
|
24
|
-
objects and a subset, together with its statistical significance.
|
25
|
-
}
|
26
|
-
\author{
|
27
|
-
Luis M. Rodriguez-R [aut, cre]
|
28
|
-
}
|
@@ -1,27 +0,0 @@
|
|
1
|
-
% Generated by roxygen2: do not edit by hand
|
2
|
-
% Please edit documentation in R/utils.R
|
3
|
-
\name{enve.truncate}
|
4
|
-
\alias{enve.truncate}
|
5
|
-
\title{Enveomics: Truncate}
|
6
|
-
\usage{
|
7
|
-
enve.truncate(x, f = 0.95, FUN = mean)
|
8
|
-
}
|
9
|
-
\arguments{
|
10
|
-
\item{x}{A vector of numbers.}
|
11
|
-
|
12
|
-
\item{f}{The fraction of values to retain.}
|
13
|
-
|
14
|
-
\item{FUN}{Summary function to apply to the vectors. To obtain the
|
15
|
-
truncated vector itself, use \code{c}.}
|
16
|
-
}
|
17
|
-
\value{
|
18
|
-
Returns the summary \code{(FUN)} of the truncated vector.
|
19
|
-
}
|
20
|
-
\description{
|
21
|
-
Removes the \code{n} highest and lowest values from a vector, and applies
|
22
|
-
summary function. The value of \code{n} is determined such that the central
|
23
|
-
range is used, corresponding to the \code{f} fraction of values.
|
24
|
-
}
|
25
|
-
\author{
|
26
|
-
Luis M. Rodriguez-R [aut, cre]
|
27
|
-
}
|
@@ -1,14 +0,0 @@
|
|
1
|
-
|
2
|
-
\name{growth.curves}
|
3
|
-
\docType{data}
|
4
|
-
\alias{growth.curves}
|
5
|
-
\title{Bacterial growth curves for three Escherichia coli mutants}
|
6
|
-
\description{
|
7
|
-
This data set provides time (first column) and three triplicated growth
|
8
|
-
curves as optical density at 600nm (OD_600nm) for different mutants of E.
|
9
|
-
coli.
|
10
|
-
}
|
11
|
-
\usage{growth.curves}
|
12
|
-
\format{A data frame with 16 rows (times) and 10 rows (times and OD_600nm).}
|
13
|
-
\keyword{datasets}
|
14
|
-
|