miga-base 0.3.0.0 → 0.3.0.1

Sign up to get free protection for your applications and to get access to all the features.
Files changed (260) hide show
  1. checksums.yaml +4 -4
  2. data/README.md +21 -4
  3. data/actions/init.rb +258 -0
  4. data/actions/run_local.rb +1 -2
  5. data/actions/test_taxonomy.rb +4 -1
  6. data/bin/miga +8 -1
  7. data/lib/miga/dataset.rb +4 -4
  8. data/lib/miga/dataset_result.rb +7 -4
  9. data/lib/miga/version.rb +2 -2
  10. data/scripts/_distances_noref_nomulti.bash +3 -1
  11. data/scripts/clade_finding.bash +1 -1
  12. data/scripts/init.bash +1 -1
  13. data/scripts/miga.bash +1 -1
  14. data/scripts/mytaxa.bash +78 -72
  15. data/scripts/mytaxa_scan.bash +67 -62
  16. data/scripts/ogs.bash +1 -1
  17. data/scripts/trimmed_fasta.bash +4 -3
  18. data/utils/enveomics/Examples/aai-matrix.bash +66 -0
  19. data/utils/enveomics/Examples/ani-matrix.bash +66 -0
  20. data/utils/enveomics/Examples/essential-phylogeny.bash +105 -0
  21. data/utils/enveomics/Examples/unus-genome-phylogeny.bash +100 -0
  22. data/utils/enveomics/LICENSE.txt +73 -0
  23. data/utils/enveomics/Makefile +52 -0
  24. data/utils/enveomics/Manifest/Tasks/aasubs.json +103 -0
  25. data/utils/enveomics/Manifest/Tasks/blasttab.json +703 -0
  26. data/utils/enveomics/Manifest/Tasks/distances.json +161 -0
  27. data/utils/enveomics/Manifest/Tasks/fasta.json +571 -0
  28. data/utils/enveomics/Manifest/Tasks/fastq.json +208 -0
  29. data/utils/enveomics/Manifest/Tasks/graphics.json +126 -0
  30. data/utils/enveomics/Manifest/Tasks/ogs.json +339 -0
  31. data/utils/enveomics/Manifest/Tasks/other.json +746 -0
  32. data/utils/enveomics/Manifest/Tasks/remote.json +355 -0
  33. data/utils/enveomics/Manifest/Tasks/sequence-identity.json +454 -0
  34. data/utils/enveomics/Manifest/Tasks/tables.json +308 -0
  35. data/utils/enveomics/Manifest/Tasks/trees.json +68 -0
  36. data/utils/enveomics/Manifest/Tasks/variants.json +111 -0
  37. data/utils/enveomics/Manifest/categories.json +132 -0
  38. data/utils/enveomics/Manifest/examples.json +154 -0
  39. data/utils/enveomics/Manifest/tasks.json +4 -0
  40. data/utils/enveomics/Pipelines/assembly.pbs/CONFIG.mock.bash +69 -0
  41. data/utils/enveomics/Pipelines/assembly.pbs/FastA.N50.pl +56 -0
  42. data/utils/enveomics/Pipelines/assembly.pbs/FastA.filterN.pl +60 -0
  43. data/utils/enveomics/Pipelines/assembly.pbs/FastA.length.pl +38 -0
  44. data/utils/enveomics/Pipelines/assembly.pbs/README.md +189 -0
  45. data/utils/enveomics/Pipelines/assembly.pbs/RUNME-2.bash +112 -0
  46. data/utils/enveomics/Pipelines/assembly.pbs/RUNME-3.bash +23 -0
  47. data/utils/enveomics/Pipelines/assembly.pbs/RUNME-4.bash +44 -0
  48. data/utils/enveomics/Pipelines/assembly.pbs/RUNME.bash +50 -0
  49. data/utils/enveomics/Pipelines/assembly.pbs/kSelector.R +37 -0
  50. data/utils/enveomics/Pipelines/assembly.pbs/newbler.pbs +68 -0
  51. data/utils/enveomics/Pipelines/assembly.pbs/newbler_preparator.pl +49 -0
  52. data/utils/enveomics/Pipelines/assembly.pbs/soap.pbs +80 -0
  53. data/utils/enveomics/Pipelines/assembly.pbs/stats.pbs +57 -0
  54. data/utils/enveomics/Pipelines/assembly.pbs/velvet.pbs +63 -0
  55. data/utils/enveomics/Pipelines/blast.pbs/01.pbs.bash +38 -0
  56. data/utils/enveomics/Pipelines/blast.pbs/02.pbs.bash +73 -0
  57. data/utils/enveomics/Pipelines/blast.pbs/03.pbs.bash +21 -0
  58. data/utils/enveomics/Pipelines/blast.pbs/BlastTab.recover_job.pl +72 -0
  59. data/utils/enveomics/Pipelines/blast.pbs/CONFIG.mock.bash +98 -0
  60. data/utils/enveomics/Pipelines/blast.pbs/FastA.split.pl +55 -0
  61. data/utils/enveomics/Pipelines/blast.pbs/README.md +127 -0
  62. data/utils/enveomics/Pipelines/blast.pbs/RUNME.bash +109 -0
  63. data/utils/enveomics/Pipelines/blast.pbs/TASK.check.bash +128 -0
  64. data/utils/enveomics/Pipelines/blast.pbs/TASK.dry.bash +16 -0
  65. data/utils/enveomics/Pipelines/blast.pbs/TASK.eo.bash +22 -0
  66. data/utils/enveomics/Pipelines/blast.pbs/TASK.pause.bash +26 -0
  67. data/utils/enveomics/Pipelines/blast.pbs/TASK.run.bash +89 -0
  68. data/utils/enveomics/Pipelines/blast.pbs/sentinel.pbs.bash +29 -0
  69. data/utils/enveomics/Pipelines/idba.pbs/README.md +49 -0
  70. data/utils/enveomics/Pipelines/idba.pbs/RUNME.bash +95 -0
  71. data/utils/enveomics/Pipelines/idba.pbs/run.pbs +56 -0
  72. data/utils/enveomics/Pipelines/trim.pbs/README.md +54 -0
  73. data/utils/enveomics/Pipelines/trim.pbs/RUNME.bash +70 -0
  74. data/utils/enveomics/Pipelines/trim.pbs/run.pbs +130 -0
  75. data/utils/enveomics/README.md +40 -0
  76. data/utils/enveomics/Scripts/AAsubs.log2ratio.rb +171 -0
  77. data/utils/enveomics/Scripts/Aln.cat.rb +162 -0
  78. data/utils/enveomics/Scripts/Aln.convert.pl +35 -0
  79. data/utils/enveomics/Scripts/AlphaDiversity.pl +152 -0
  80. data/utils/enveomics/Scripts/BlastPairwise.AAsubs.pl +102 -0
  81. data/utils/enveomics/Scripts/BlastTab.addlen.rb +61 -0
  82. data/utils/enveomics/Scripts/BlastTab.advance.bash +48 -0
  83. data/utils/enveomics/Scripts/BlastTab.best_hit_sorted.pl +55 -0
  84. data/utils/enveomics/Scripts/BlastTab.catsbj.pl +106 -0
  85. data/utils/enveomics/Scripts/BlastTab.cogCat.rb +76 -0
  86. data/utils/enveomics/Scripts/BlastTab.filter.pl +47 -0
  87. data/utils/enveomics/Scripts/BlastTab.kegg_pep2path_rest.pl +194 -0
  88. data/utils/enveomics/Scripts/BlastTab.metaxaPrep.pl +104 -0
  89. data/utils/enveomics/Scripts/BlastTab.pairedHits.rb +157 -0
  90. data/utils/enveomics/Scripts/BlastTab.recplot2.R +40 -0
  91. data/utils/enveomics/Scripts/BlastTab.seqdepth.pl +86 -0
  92. data/utils/enveomics/Scripts/BlastTab.seqdepth_ZIP.pl +119 -0
  93. data/utils/enveomics/Scripts/BlastTab.seqdepth_nomedian.pl +86 -0
  94. data/utils/enveomics/Scripts/BlastTab.subsample.pl +47 -0
  95. data/utils/enveomics/Scripts/BlastTab.sumPerHit.pl +114 -0
  96. data/utils/enveomics/Scripts/BlastTab.taxid2taxrank.pl +90 -0
  97. data/utils/enveomics/Scripts/BlastTab.topHits_sorted.rb +101 -0
  98. data/utils/enveomics/Scripts/Chao1.pl +97 -0
  99. data/utils/enveomics/Scripts/CharTable.classify.rb +234 -0
  100. data/utils/enveomics/Scripts/EBIseq2tax.rb +83 -0
  101. data/utils/enveomics/Scripts/FastA.N50.pl +56 -0
  102. data/utils/enveomics/Scripts/FastA.filter.pl +52 -0
  103. data/utils/enveomics/Scripts/FastA.filterLen.pl +28 -0
  104. data/utils/enveomics/Scripts/FastA.filterN.pl +60 -0
  105. data/utils/enveomics/Scripts/FastA.fragment.rb +92 -0
  106. data/utils/enveomics/Scripts/FastA.gc.pl +42 -0
  107. data/utils/enveomics/Scripts/FastA.interpose.pl +87 -0
  108. data/utils/enveomics/Scripts/FastA.length.pl +38 -0
  109. data/utils/enveomics/Scripts/FastA.per_file.pl +36 -0
  110. data/utils/enveomics/Scripts/FastA.qlen.pl +57 -0
  111. data/utils/enveomics/Scripts/FastA.rename.pl +65 -0
  112. data/utils/enveomics/Scripts/FastA.revcom.pl +23 -0
  113. data/utils/enveomics/Scripts/FastA.slider.pl +85 -0
  114. data/utils/enveomics/Scripts/FastA.split.pl +55 -0
  115. data/utils/enveomics/Scripts/FastA.subsample.pl +131 -0
  116. data/utils/enveomics/Scripts/FastA.tag.rb +64 -0
  117. data/utils/enveomics/Scripts/FastA.wrap.rb +48 -0
  118. data/utils/enveomics/Scripts/FastQ.filter.pl +54 -0
  119. data/utils/enveomics/Scripts/FastQ.interpose.pl +90 -0
  120. data/utils/enveomics/Scripts/FastQ.offset.pl +90 -0
  121. data/utils/enveomics/Scripts/FastQ.split.pl +53 -0
  122. data/utils/enveomics/Scripts/FastQ.tag.rb +63 -0
  123. data/utils/enveomics/Scripts/FastQ.toFastA.awk +24 -0
  124. data/utils/enveomics/Scripts/GenBank.add_fields.rb +84 -0
  125. data/utils/enveomics/Scripts/HMM.essential.rb +254 -0
  126. data/utils/enveomics/Scripts/HMMsearch.extractIds.rb +83 -0
  127. data/utils/enveomics/Scripts/JPlace.distances.rb +88 -0
  128. data/utils/enveomics/Scripts/JPlace.to_iToL.rb +306 -0
  129. data/utils/enveomics/Scripts/M5nr.getSequences.rb +81 -0
  130. data/utils/enveomics/Scripts/MeTaxa.distribution.pl +198 -0
  131. data/utils/enveomics/Scripts/MyTaxa.fragsByTax.pl +35 -0
  132. data/utils/enveomics/Scripts/MyTaxa.seq-taxrank.rb +49 -0
  133. data/utils/enveomics/Scripts/NCBIacc2tax.rb +92 -0
  134. data/utils/enveomics/Scripts/Newick.autoprune.R +27 -0
  135. data/utils/enveomics/Scripts/RAxML-EPA.to_iToL.pl +228 -0
  136. data/utils/enveomics/Scripts/RefSeq.download.bash +48 -0
  137. data/utils/enveomics/Scripts/SRA.download.bash +50 -0
  138. data/utils/enveomics/Scripts/TRIBS.plot-test.R +36 -0
  139. data/utils/enveomics/Scripts/TRIBS.test.R +39 -0
  140. data/utils/enveomics/Scripts/Table.barplot.R +30 -0
  141. data/utils/enveomics/Scripts/Table.df2dist.R +30 -0
  142. data/utils/enveomics/Scripts/Table.filter.pl +61 -0
  143. data/utils/enveomics/Scripts/Table.merge.pl +77 -0
  144. data/utils/enveomics/Scripts/Table.replace.rb +69 -0
  145. data/utils/enveomics/Scripts/Table.round.rb +63 -0
  146. data/utils/enveomics/Scripts/Table.split.pl +57 -0
  147. data/utils/enveomics/Scripts/Taxonomy.silva2ncbi.rb +227 -0
  148. data/utils/enveomics/Scripts/VCF.KaKs.rb +147 -0
  149. data/utils/enveomics/Scripts/VCF.SNPs.rb +88 -0
  150. data/utils/enveomics/Scripts/aai.rb +373 -0
  151. data/utils/enveomics/Scripts/ani.rb +362 -0
  152. data/utils/enveomics/Scripts/gi2tax.rb +103 -0
  153. data/utils/enveomics/Scripts/in_silico_GA_GI.pl +96 -0
  154. data/utils/enveomics/Scripts/lib/data/essential.hmm.gz +0 -0
  155. data/utils/enveomics/Scripts/lib/enveomics_rb/enveomics.rb +26 -0
  156. data/utils/enveomics/Scripts/lib/enveomics_rb/jplace.rb +253 -0
  157. data/utils/enveomics/Scripts/lib/enveomics_rb/og.rb +182 -0
  158. data/utils/enveomics/Scripts/lib/enveomics_rb/remote_data.rb +74 -0
  159. data/utils/enveomics/Scripts/lib/enveomics_rb/seq_range.rb +237 -0
  160. data/utils/enveomics/Scripts/lib/enveomics_rb/stat.rb +30 -0
  161. data/utils/enveomics/Scripts/lib/enveomics_rb/vcf.rb +135 -0
  162. data/utils/enveomics/Scripts/ogs.annotate.rb +88 -0
  163. data/utils/enveomics/Scripts/ogs.core-pan.rb +160 -0
  164. data/utils/enveomics/Scripts/ogs.extract.rb +125 -0
  165. data/utils/enveomics/Scripts/ogs.mcl.rb +186 -0
  166. data/utils/enveomics/Scripts/ogs.rb +104 -0
  167. data/utils/enveomics/Scripts/ogs.stats.rb +131 -0
  168. data/utils/enveomics/Scripts/rbm.rb +137 -0
  169. data/utils/enveomics/Tests/Makefile +10 -0
  170. data/utils/enveomics/Tests/Mgen_M2288.faa +3189 -0
  171. data/utils/enveomics/Tests/Mgen_M2288.fna +8282 -0
  172. data/utils/enveomics/Tests/Mgen_M2321.fna +8288 -0
  173. data/utils/enveomics/Tests/Nequ_Kin4M.faa +2970 -0
  174. data/utils/enveomics/Tests/Xanthomonas_oryzae-PilA.tribs.Rdata +0 -0
  175. data/utils/enveomics/Tests/Xanthomonas_oryzae-PilA.txt +7 -0
  176. data/utils/enveomics/Tests/Xanthomonas_oryzae.aai-mat.tsv +17 -0
  177. data/utils/enveomics/Tests/Xanthomonas_oryzae.aai.tsv +137 -0
  178. data/utils/enveomics/Tests/a_mg.cds-go.blast.tsv +123 -0
  179. data/utils/enveomics/Tests/a_mg.reads-cds.blast.tsv +200 -0
  180. data/utils/enveomics/Tests/a_mg.reads-cds.counts.tsv +55 -0
  181. data/utils/enveomics/Tests/alkB.nwk +1 -0
  182. data/utils/enveomics/Tests/anthrax-cansnp-data.tsv +13 -0
  183. data/utils/enveomics/Tests/anthrax-cansnp-key.tsv +17 -0
  184. data/utils/enveomics/Tests/hiv1.faa +59 -0
  185. data/utils/enveomics/Tests/hiv1.fna +134 -0
  186. data/utils/enveomics/Tests/hiv2.faa +70 -0
  187. data/utils/enveomics/Tests/hiv_mix-hiv1.blast.tsv +233 -0
  188. data/utils/enveomics/Tests/hiv_mix-hiv1.blast.tsv.lim +1 -0
  189. data/utils/enveomics/Tests/hiv_mix-hiv1.blast.tsv.rec +233 -0
  190. data/utils/enveomics/Tests/phyla_counts.tsv +10 -0
  191. data/utils/enveomics/Tests/primate_lentivirus.ogs +11 -0
  192. data/utils/enveomics/Tests/primate_lentivirus.rbm/hiv1-hiv1.rbm +9 -0
  193. data/utils/enveomics/Tests/primate_lentivirus.rbm/hiv1-hiv2.rbm +8 -0
  194. data/utils/enveomics/Tests/primate_lentivirus.rbm/hiv1-siv.rbm +6 -0
  195. data/utils/enveomics/Tests/primate_lentivirus.rbm/hiv2-hiv2.rbm +9 -0
  196. data/utils/enveomics/Tests/primate_lentivirus.rbm/hiv2-siv.rbm +6 -0
  197. data/utils/enveomics/Tests/primate_lentivirus.rbm/siv-siv.rbm +6 -0
  198. data/utils/enveomics/build_enveomics_r.bash +44 -0
  199. data/utils/enveomics/enveomics.R/DESCRIPTION +31 -0
  200. data/utils/enveomics/enveomics.R/NAMESPACE +35 -0
  201. data/utils/enveomics/enveomics.R/R/autoprune.R +121 -0
  202. data/utils/enveomics/enveomics.R/R/barplot.R +165 -0
  203. data/utils/enveomics/enveomics.R/R/cliopts.R +119 -0
  204. data/utils/enveomics/enveomics.R/R/df2dist.R +117 -0
  205. data/utils/enveomics/enveomics.R/R/growthcurve.R +263 -0
  206. data/utils/enveomics/enveomics.R/R/recplot.R +320 -0
  207. data/utils/enveomics/enveomics.R/R/recplot2.R +745 -0
  208. data/utils/enveomics/enveomics.R/R/tribs.R +423 -0
  209. data/utils/enveomics/enveomics.R/R/utils.R +16 -0
  210. data/utils/enveomics/enveomics.R/README.md +52 -0
  211. data/utils/enveomics/enveomics.R/data/growth.curves.rda +0 -0
  212. data/utils/enveomics/enveomics.R/data/phyla.counts.rda +0 -0
  213. data/utils/enveomics/enveomics.R/man/enve.GrowthCurve-class.Rd +30 -0
  214. data/utils/enveomics/enveomics.R/man/enve.TRIBS-class.Rd +43 -0
  215. data/utils/enveomics/enveomics.R/man/enve.TRIBS.merge.Rd +19 -0
  216. data/utils/enveomics/enveomics.R/man/enve.TRIBStest-class.Rd +37 -0
  217. data/utils/enveomics/enveomics.R/man/enve.__prune.iter.Rd +24 -0
  218. data/utils/enveomics/enveomics.R/man/enve.__prune.reduce.Rd +24 -0
  219. data/utils/enveomics/enveomics.R/man/enve.__tribs.Rd +33 -0
  220. data/utils/enveomics/enveomics.R/man/enve.barplot.Rd +64 -0
  221. data/utils/enveomics/enveomics.R/man/enve.cliopts.Rd +37 -0
  222. data/utils/enveomics/enveomics.R/man/enve.col.alpha.Rd +19 -0
  223. data/utils/enveomics/enveomics.R/man/enve.col2alpha.Rd +18 -0
  224. data/utils/enveomics/enveomics.R/man/enve.df2dist.Rd +26 -0
  225. data/utils/enveomics/enveomics.R/man/enve.df2dist.group.Rd +25 -0
  226. data/utils/enveomics/enveomics.R/man/enve.df2dist.list.Rd +26 -0
  227. data/utils/enveomics/enveomics.R/man/enve.growthcurve.Rd +49 -0
  228. data/utils/enveomics/enveomics.R/man/enve.prune.dist.Rd +28 -0
  229. data/utils/enveomics/enveomics.R/man/enve.recplot.Rd +97 -0
  230. data/utils/enveomics/enveomics.R/man/enve.recplot2-class.Rd +40 -0
  231. data/utils/enveomics/enveomics.R/man/enve.recplot2.Rd +40 -0
  232. data/utils/enveomics/enveomics.R/man/enve.recplot2.__counts.Rd +24 -0
  233. data/utils/enveomics/enveomics.R/man/enve.recplot2.__findPeak.Rd +40 -0
  234. data/utils/enveomics/enveomics.R/man/enve.recplot2.__findPeaks.Rd +18 -0
  235. data/utils/enveomics/enveomics.R/man/enve.recplot2.__peakHist.Rd +22 -0
  236. data/utils/enveomics/enveomics.R/man/enve.recplot2.__whichClosestPeak.Rd +20 -0
  237. data/utils/enveomics/enveomics.R/man/enve.recplot2.changeCutoff.Rd +18 -0
  238. data/utils/enveomics/enveomics.R/man/enve.recplot2.corePeak.Rd +18 -0
  239. data/utils/enveomics/enveomics.R/man/enve.recplot2.extractWindows.Rd +27 -0
  240. data/utils/enveomics/enveomics.R/man/enve.recplot2.findPeaks.Rd +53 -0
  241. data/utils/enveomics/enveomics.R/man/enve.recplot2.peak-class.Rd +44 -0
  242. data/utils/enveomics/enveomics.R/man/enve.tribs.Rd +44 -0
  243. data/utils/enveomics/enveomics.R/man/enve.tribs.test.Rd +21 -0
  244. data/utils/enveomics/enveomics.R/man/enveomics.R-package.Rd +15 -0
  245. data/utils/enveomics/enveomics.R/man/growth.curves.Rd +14 -0
  246. data/utils/enveomics/enveomics.R/man/phyla.counts.Rd +13 -0
  247. data/utils/enveomics/enveomics.R/man/plot.enve.GrowthCurve.Rd +43 -0
  248. data/utils/enveomics/enveomics.R/man/plot.enve.TRIBS.Rd +29 -0
  249. data/utils/enveomics/enveomics.R/man/plot.enve.TRIBStest.Rd +30 -0
  250. data/utils/enveomics/enveomics.R/man/plot.enve.recplot2.Rd +71 -0
  251. data/utils/enveomics/enveomics.R/man/summary.enve.GrowthCurve.Rd +18 -0
  252. data/utils/enveomics/enveomics.R/man/summary.enve.TRIBS.Rd +18 -0
  253. data/utils/enveomics/enveomics.R/man/summary.enve.TRIBStest.Rd +18 -0
  254. data/utils/enveomics/enveomics.R/man/z$-methods.Rd +27 -0
  255. data/utils/enveomics/globals.mk +8 -0
  256. data/utils/enveomics/manifest.json +9 -0
  257. data/utils/index_metadata.rb +0 -0
  258. data/utils/plot-taxdist.R +0 -0
  259. data/utils/requirements.txt +19 -19
  260. metadata +242 -2
@@ -0,0 +1,30 @@
1
+ \name{enve.GrowthCurve-class}
2
+ \Rdversion{1.1}
3
+ \docType{class}
4
+ \alias{enve.GrowthCurve-class}
5
+ \alias{enve.GrowthCurve}
6
+ %% \alias{$}
7
+
8
+ \title{enve.GrowthCurve S4 class}
9
+ \description{Enve-omics representation of fitted growth curves.}
10
+ \section{Objects from the Class}{Objects can be created by calls of the form \code{new(enve.GrowthCurve ...)}}
11
+ \section{Slots}{
12
+ \describe{
13
+ \item{\code{design}:}{(\code{array}) Experimental design of the experiment.}
14
+ \item{\code{models}:}{(\code{list}) Fitted growth curve models.}
15
+ \item{\code{predict}:}{(\code{list}) Fitted growth curve values.}
16
+ \item{\code{call}:}{(\code{call}) Call producing this object.}
17
+ }
18
+ }
19
+ \section{Methods}{
20
+ \describe{
21
+ \item{$}{\code{signature(x = "enve.GrowthCurve")}: ... }
22
+ }
23
+ }
24
+
25
+ \author{Luis M. Rodriguez-R [aut, cre]}
26
+
27
+
28
+
29
+
30
+
@@ -0,0 +1,43 @@
1
+ \name{enve.TRIBS-class}
2
+ \Rdversion{1.1}
3
+ \docType{class}
4
+ \alias{enve.TRIBS-class}
5
+ \alias{enve.TRIBS}
6
+
7
+ \title{enve.TRIBS S4 class}
8
+ \description{Enve-omics representation of "Transformed-space Resampling In Biased Sets
9
+ (TRIBS)". This object represents sets of distances between objects,
10
+ sampled nearly-uniformly at random in "distance space". Subsampling
11
+ without selection is trivial, since both the distances space and the
12
+ selection occur in the same transformed space. However, it's useful to
13
+ compare randomly subsampled sets against a selected set of objects. This
14
+ is intended to identify overdispersion or overclustering (see
15
+ `enve.TRIBStest`) of a subset against the entire collection of objects
16
+ with minimum impact of sampling biases. This object can be produced by
17
+ `enve.tribs` and supports S4 methods `plot` and `summary`.}
18
+ \section{Objects from the Class}{Objects can be created by calls of the form \code{new(enve.TRIBS ...)}}
19
+ \section{Slots}{
20
+ \describe{
21
+ \item{\code{distance}:}{(\code{numeric}) Centrality measurement of the distances between the
22
+ << selected objects (without subsampling).}
23
+ \item{\code{points}:}{(\code{matrix}) Position of the different objects in distance
24
+ << space.}
25
+ \item{\code{distances}:}{(\code{matrix}) Subsampled distances, where the rows are replicates
26
+ << and the columns are subsampling levels.}
27
+ \item{\code{spaceSize}:}{(\code{numeric}) Number of objects.}
28
+ \item{\code{selSize}:}{(\code{numeric}) Number of selected objects.}
29
+ \item{\code{dimensions}:}{(\code{numeric}) Number of dimensions in the distance space.}
30
+ \item{\code{subsamples}:}{(\code{numeric}) Subsampling levels (as fractions, from 0 to 1).}
31
+ \item{\code{call}:}{(\code{call}) Call producing this object.}
32
+ }
33
+ }
34
+ \section{Methods}{
35
+ No methods defined with class "enve.TRIBS" in the signature.
36
+ }
37
+
38
+ \author{Luis M. Rodriguez-R [aut, cre]}
39
+
40
+
41
+
42
+
43
+
@@ -0,0 +1,19 @@
1
+ \name{enve.TRIBS.merge}
2
+ \alias{enve.TRIBS.merge}
3
+ \title{enve TRIBS merge}
4
+ \description{Merges two `enve.TRIBS` objects generated from the same objects at
5
+ different subsampling levels.}
6
+ \usage{enve.TRIBS.merge(x, y)}
7
+ \arguments{
8
+ \item{x}{First `enve.TRIBS` object.}
9
+ \item{y}{Second `enve.TRIBS` object.}
10
+ }
11
+
12
+ \value{Returns an `enve.TRIBS` object.}
13
+
14
+ \author{Luis M. Rodriguez-R [aut, cre]}
15
+
16
+
17
+
18
+
19
+
@@ -0,0 +1,37 @@
1
+ \name{enve.TRIBStest-class}
2
+ \Rdversion{1.1}
3
+ \docType{class}
4
+ \alias{enve.TRIBStest-class}
5
+ \alias{enve.TRIBStest}
6
+
7
+ \title{enve.TRIBStest S4 class}
8
+ \description{Test of significance of overclustering or overdispersion in a selected
9
+ set of objects with respect to the entire set (see `enve.TRIBS`). This
10
+ object can be produced by `enve.tribs.test` and supports S4 methods
11
+ `plot` and `summary`.}
12
+ \section{Objects from the Class}{Objects can be created by calls of the form \code{new(enve.TRIBStest ...)}}
13
+ \section{Slots}{
14
+ \describe{
15
+ \item{\code{pval.gt}:}{(\code{numeric}) P-value for the overdispersion test.}
16
+ \item{\code{pval.lt}:}{(\code{numeric}) P-value for the overclustering test.}
17
+ \item{\code{all.dist}:}{(\code{numeric}) Empiric PDF of distances for the entire dataset (subsampled at selection
18
+ size).}
19
+ \item{\code{sel.dist}:}{(\code{numeric}) Empiric PDF of distances for the selected objects (without subsampling).}
20
+ \item{\code{diff.dist}:}{(\code{numeric}) Empiric PDF of the difference between `all.dist` and `sel.dist`. The
21
+ p-values are estimating by comparing areas in this PDF greater than and
22
+ lesser than zero.}
23
+ \item{\code{dist.mids}:}{(\code{numeric}) Midpoints of the empiric PDFs of distances.}
24
+ \item{\code{diff.mids}:}{(\code{numeric}) Midpoints of the empiric PDF of difference of distances.}
25
+ \item{\code{call}:}{(\code{call}) Call producing this object.}
26
+ }
27
+ }
28
+ \section{Methods}{
29
+ No methods defined with class "enve.TRIBStest" in the signature.
30
+ }
31
+
32
+ \author{Luis M. Rodriguez-R [aut, cre]}
33
+
34
+
35
+
36
+
37
+
@@ -0,0 +1,24 @@
1
+ \name{enve.__prune.iter}
2
+ \alias{enve.__prune.iter}
3
+ \title{enve prune iter}
4
+ \description{Internal function for enve.prune.dist}
5
+ \usage{enve.__prune.iter(t, dist, min_dist, quiet)}
6
+ \arguments{
7
+ \item{t}{
8
+ }
9
+ \item{dist}{
10
+ }
11
+ \item{min_dist}{
12
+ }
13
+ \item{quiet}{
14
+ }
15
+ }
16
+
17
+
18
+
19
+ \author{Luis M. Rodriguez-R [aut, cre]}
20
+
21
+
22
+
23
+
24
+
@@ -0,0 +1,24 @@
1
+ \name{enve.__prune.reduce}
2
+ \alias{enve.__prune.reduce}
3
+ \title{enve prune reduce}
4
+ \description{Internal function for enve.prune.dist}
5
+ \usage{enve.__prune.reduce(t, nodes, min_dist, quiet)}
6
+ \arguments{
7
+ \item{t}{
8
+ }
9
+ \item{nodes}{
10
+ }
11
+ \item{min_dist}{
12
+ }
13
+ \item{quiet}{
14
+ }
15
+ }
16
+
17
+
18
+
19
+ \author{Luis M. Rodriguez-R [aut, cre]}
20
+
21
+
22
+
23
+
24
+
@@ -0,0 +1,33 @@
1
+ \name{enve.__tribs}
2
+ \alias{enve.__tribs}
3
+ \title{enve tribs}
4
+ \description{Internal ancilliary function (see `enve.tribs`).}
5
+ \usage{enve.__tribs(rep, frx, selection, dimensions, dots, dist.method,
6
+ summary.fx, dist)}
7
+ \arguments{
8
+ \item{rep}{
9
+ }
10
+ \item{frx}{
11
+ }
12
+ \item{selection}{
13
+ }
14
+ \item{dimensions}{
15
+ }
16
+ \item{dots}{
17
+ }
18
+ \item{dist.method}{
19
+ }
20
+ \item{summary.fx}{
21
+ }
22
+ \item{dist}{
23
+ }
24
+ }
25
+
26
+
27
+
28
+ \author{Luis M. Rodriguez-R [aut, cre]}
29
+
30
+
31
+
32
+
33
+
@@ -0,0 +1,64 @@
1
+ \name{enve.barplot}
2
+ \alias{enve.barplot}
3
+ \title{enve barplot}
4
+ \description{Creates nice barplots from tab-delimited tables}
5
+ \usage{enve.barplot(x, sizes, top = 25, colors.per.group = 9, bars.width = 4,
6
+ legend.ncol = 1, other.col = "#000000", add.trend = FALSE,
7
+ organic.trend = FALSE, sort.by = median, min.report = 101,
8
+ order = NULL, col, ...)}
9
+ \arguments{
10
+ \item{x}{Can be either the input data or the path to the file containing the
11
+ table. If it contains the data, it must be a data frame or an object
12
+ coercible to data frame. If it is a path, it must point to a
13
+ tab-delimited file containing header (first row) and row names (first
14
+ column).}
15
+ \item{sizes}{A numeric vector containing the real size of the samples (columns) in
16
+ the same order of the input table. If set, the values are assumed to be
17
+ 100\%, otherwise the sum of the columns is used.}
18
+ \item{top}{Maximum number of categories to display. Any additional categories will
19
+ be listed as "Others".}
20
+ \item{colors.per.group}{Number of categories in the first two saturation groups of colors. The
21
+ third group contains the remaining categories if needed.}
22
+ \item{bars.width}{Width of the barplot with respect to the legend.}
23
+ \item{legend.ncol}{Number of columns in the legend.}
24
+ \item{other.col}{Color of the "Others" category.}
25
+ \item{add.trend}{Controls if semi-transparent areas are to be plotted between the bars to
26
+ connect the regions (trend regions).}
27
+ \item{organic.trend}{Controls if the trend regions are to be smoothed (curves). By default,
28
+ trend regions have straight edges. If TRUE, forces add.trend=TRUE.}
29
+ \item{sort.by}{Any function that takes a numeric vector and returns a numeric scalar.
30
+ This function is applied to each row, and the resulting values are used
31
+ to sort the rows (decreasingly). Good options include: sd, min, max,
32
+ mean, median.}
33
+ \item{min.report}{Minimum percentage to report the value in the plot. Any value above 100
34
+ indicates that no values are to be reported.}
35
+ \item{order}{Controls how the rows should be ordered. If NULL (default), sort.by is
36
+ applied per row and the results are sorted decreasingly. If NA, no
37
+ sorting is performed, i.e., the original order is respected. If a vector
38
+ is provided, it is assumed to be the custom order to be used (either by
39
+ numeric index or by row names).}
40
+ \item{col}{Colors to use. If provided, overrides the variables `top` and
41
+ `colors.per.group`, but `other.col` is still used if the vector is
42
+ insufficient for all the rows. An additional palette is available with
43
+ col='coto' (contributed by Luis (Coto) Orellana).}
44
+ \item{\dots}{Any additional parameters to be passed to `barplot`.}
45
+ }
46
+
47
+
48
+
49
+ \author{Luis M. Rodriguez-R [aut, cre]}
50
+
51
+
52
+
53
+
54
+ \examples{
55
+ # Load data
56
+ data("phyla.counts", package="enveomics.R", envir=environment())
57
+ # Create a barplot sorted by variance with organic trends
58
+ enve.barplot(phyla.counts, # Counts of phyla in four sites
59
+ sizes=c(250,100,75,200), # Total sizes of the datasets of each site
60
+ bars.width=2, # Decrease from default, so the names are fully displayed
61
+ organic.trend=TRUE, # Nice curvy background
62
+ sort.by=var # Sort by variance across sites
63
+ )
64
+ }
@@ -0,0 +1,37 @@
1
+ \name{enve.cliopts}
2
+ \alias{enve.cliopts}
3
+ \title{enve cliopts}
4
+ \description{Generates nicely formatted command-line interfaces for
5
+ functions (_closures_ only).}
6
+ \usage{enve.cliopts(fx, rd_file, positional_arguments, usage, mandatory = c(),
7
+ vectorize = c(), ignore = c(), number = c(), defaults = list(),
8
+ o_desc = list(), p_desc = "")}
9
+ \arguments{
10
+ \item{fx}{Function for which the interface should be generated.}
11
+ \item{rd_file}{(Optional) .Rd file with the standard documentation of the function.}
12
+ \item{positional_arguments}{(Optional) Number of _positional_ arguments passed to parse_args
13
+ (package:optparse).}
14
+ \item{usage}{(Optional) Usage passed to OptionParser (package:optparse).}
15
+ \item{mandatory}{Mandatory arguments.}
16
+ \item{vectorize}{Arguments of the function to vectorize (comma-delimited). If numeric,
17
+ use also `number`.}
18
+ \item{ignore}{Arguments of the function to ignore.}
19
+ \item{number}{Force these arguments as numerics. Useful for numeric
20
+ vectors (see `vectorize`) or arguments with no defaults.}
21
+ \item{defaults}{Defaults to use instead of the ones provided by the formals.}
22
+ \item{o_desc}{Descriptions of the options. Help from `rd` is ignored for arguments
23
+ present in this list.}
24
+ \item{p_desc}{Description of the function. Help from `rd` is ignored for the
25
+ function description unless this value is an empty string.}
26
+ }
27
+
28
+ \value{Returns a `list` with keys: `options`, a named list with the values for
29
+ the function's arguments; and `args`, a vector with zero or more strings
30
+ containing the positional arguments.}
31
+
32
+ \author{Luis M. Rodriguez-R [aut, cre]}
33
+
34
+
35
+
36
+
37
+
@@ -0,0 +1,19 @@
1
+ \name{enve.col.alpha}
2
+ \alias{enve.col.alpha}
3
+ \title{enve col alpha}
4
+ \description{Modify alpha in a color (or vector of colors).}
5
+ \usage{enve.col.alpha(col, alpha = 1/2)}
6
+ \arguments{
7
+ \item{col}{Color or vector of colors. It can be any value supported by `col2rgb`,
8
+ such as 'darkred' or '#009988'.}
9
+ \item{alpha}{Alpha value to add to the color, from 0 to 1.}
10
+ }
11
+
12
+ \value{Returns a color or a vector of colors in hex notation including alpha.}
13
+
14
+ \author{Luis M. Rodriguez-R [aut, cre]}
15
+
16
+
17
+
18
+
19
+
@@ -0,0 +1,18 @@
1
+ \name{enve.col2alpha}
2
+ \alias{enve.col2alpha}
3
+ \title{enve col2alpha}
4
+ \description{Takes a vector of colors and sets the alpha.}
5
+ \usage{enve.col2alpha(x, alpha)}
6
+ \arguments{
7
+ \item{x}{A vector of any value base colors.}
8
+ \item{alpha}{Alpha level to set (in the 0-1 range).}
9
+ }
10
+
11
+
12
+
13
+ \author{Luis M. Rodriguez-R [aut, cre]}
14
+
15
+
16
+
17
+
18
+
@@ -0,0 +1,26 @@
1
+ \name{enve.df2dist}
2
+ \alias{enve.df2dist}
3
+ \title{enve df2dist}
4
+ \description{Transform a dataframe (or coercible object, like a table) into a `dist` object.}
5
+ \usage{enve.df2dist(x, obj1.index = 1, obj2.index = 2, dist.index = 3,
6
+ default.d = NA, max.sim = 0)}
7
+ \arguments{
8
+ \item{x}{A table (or coercible object) with at least three columns: (1) ID of the object 1,
9
+ (2) ID of the object 2, and (3) distance between the two objects.}
10
+ \item{obj1.index}{Index of the column containing the ID of the object 1.}
11
+ \item{obj2.index}{Index of the column containing the ID of the object 2.}
12
+ \item{dist.index}{Index of the column containing the distance.}
13
+ \item{default.d}{Default value (for missing values)}
14
+ \item{max.sim}{If not-zero, assumes that the values are similarity (not distance)
15
+ and this is the maximum similarity (corresponding to distance 0).
16
+ Applies transformation: distance = (max.sim - values)/max.sim.}
17
+ }
18
+
19
+ \value{Returns a `dist` object.}
20
+
21
+ \author{Luis M. Rodriguez-R [aut, cre]}
22
+
23
+
24
+
25
+
26
+
@@ -0,0 +1,25 @@
1
+ \name{enve.df2dist.group}
2
+ \alias{enve.df2dist.group}
3
+ \title{enve df2dist group}
4
+ \description{Transform a dataframe (or coercible object, like a table) into a `dist` object, where
5
+ there are 1 or more distances between each pair of objects.}
6
+ \usage{enve.df2dist.group(x, obj1.index = 1, obj2.index = 2, dist.index = 3,
7
+ summary = median, empty.rm = TRUE)}
8
+ \arguments{
9
+ \item{x}{A dataframe (or coercible object) with at least three columns: (1) ID of the object 1,
10
+ (2) ID of the object 2, and (3) distance between the two objects.}
11
+ \item{obj1.index}{Index of the column containing the ID of the object 1.}
12
+ \item{obj2.index}{Index of the column containing the ID of the object 2.}
13
+ \item{dist.index}{Index of the column containing the distance.}
14
+ \item{summary}{Function summarizing the different distances between the two objects.}
15
+ \item{empty.rm}{Remove rows with empty or NA groups}
16
+ }
17
+
18
+ \value{Returns a `dist` object.}
19
+
20
+ \author{Luis M. Rodriguez-R [aut, cre]}
21
+
22
+
23
+
24
+
25
+
@@ -0,0 +1,26 @@
1
+ \name{enve.df2dist.list}
2
+ \alias{enve.df2dist.list}
3
+ \title{enve df2dist list}
4
+ \description{Transform a dataframe (or coercible object, like a table) into a `dist` object.}
5
+ \usage{enve.df2dist.list(x, groups, obj1.index = 1, obj2.index = 2,
6
+ dist.index = 3, empty.rm = TRUE, ...)}
7
+ \arguments{
8
+ \item{x}{A dataframe (or coercible object) with at least three columns: (1) ID of the object 1,
9
+ (2) ID of the object 2, and (3) distance between the two objects.}
10
+ \item{groups}{Named array where the IDs correspond to the object IDs, and the values correspond to
11
+ the group.}
12
+ \item{obj1.index}{Index of the column containing the ID of the object 1.}
13
+ \item{obj2.index}{Index of the column containing the ID of the object 2.}
14
+ \item{dist.index}{Index of the column containing the distance.}
15
+ \item{empty.rm}{Remove incomplete matrices}
16
+ \item{\dots}{Any other parameters supported by `enve.df2dist.group`.}
17
+ }
18
+
19
+ \value{Returns a `list` of `dist` object.}
20
+
21
+ \author{Luis M. Rodriguez-R [aut, cre]}
22
+
23
+
24
+
25
+
26
+