new version

DESCRIPTION

@@ -1,10 +1,10 @@
 Package: MicEco
 Title: Various functions for microbial community data
-Version: 0.9.5
+Version: 0.9.6
 Authors@R: person("Jakob", "Russel", email = "[email protected]", role = c("aut", "cre"))
 Description: Collection of functions for microbiome analyses. E.g. fitting neutral models and standardized effect sizes of phylogenetic beta diversities, and much more.
 Depends: R (>= 3.2.5)
-Imports: stats, utils, phyloseq, foreach, doSNOW, picante, vegan, snow, bbmle, Hmisc, abind, reshape2
+Imports: stats, utils, phyloseq, foreach, doSNOW, picante, vegan, snow, bbmle, Hmisc, abind, reshape2, eulerr
 License: GPL (>= 3) | file LICENSE
 Encoding: UTF-8
 LazyData: true

NAMESPACE

@@ -7,6 +7,7 @@ export(comdistnt.par)
 export(community_rrna)
 export(neutral.fit)
 export(neutral.rand)
+export(phylo_euler)
 export(proportionality)
 export(ps_refactor)
 export(rarefy_rrna)
@@ -24,6 +25,7 @@ export(ses.permtest)
 import(Hmisc)
 import(bbmle)
 import(doSNOW)
+import(eulerr)
 import(foreach)
 import(parallel)
 import(phyloseq)
@@ -32,9 +34,13 @@ import(utils)
 import(vegan)
 importFrom(abind,abind)
 importFrom(reshape2,dcast)
+importFrom(stats,aggregate)
 importFrom(stats,as.dist)
+importFrom(stats,as.formula)
 importFrom(stats,cov)
+importFrom(stats,dist)
 importFrom(stats,dnorm)
+importFrom(stats,hclust)
 importFrom(stats,na.omit)
 importFrom(stats,p.adjust)
 importFrom(stats,pbeta)

R/phylo_euler.R

@@ -0,0 +1,40 @@
+#' Make Euler (Venn) diagram of shared taxa (ASVs, OTUs) across sample groups
+#'
+#' Make Euler (Venn) diagram of shared taxa (ASVs, OTUs) across sample groups from a phyloseq object. Overlap can be weighted by relative abundance
+#' @param ps A phyloseq object
+#' @param group The grouping factor. Should match variable in sample_data(ps)
+#' @param weight If TRUE, the overlaps are weighted by abundance
+#' @param type "percent" or "count"
+#' @param relative Should abundances be made relative
+#' @param ... Additional arguments
+#' @keywords venn euler diagram
+#' @return An euler plot
+#' @import phyloseq
+#' @import eulerr
+#' @importFrom stats aggregate as.formula
+#' @export
+
+phylo_euler <- function(ps, group, weight = FALSE, type = "percent", relative = TRUE, ...){
+
+    if(relative){
+        ps <- transform_sample_counts(ps, function(x) x/sum(x))
+    }
+
+    ps_melted <- psmelt(ps)
+
+    ps_agg <- aggregate(as.formula(paste("Abundance ~ OTU +",group)), data = ps_melted, function(x) mean(x))
+
+    ps_mat <- reshape2::dcast(as.formula(paste("OTU ~ ",group)), data = ps_agg, value.var = "Abundance")
+
+    ps_mat <- ps_mat[, -1]
+    ps_mat_bin <- (ps_mat>0)*1
+
+    if(weight){
+        df <- euler(ps_mat_bin, weights = rowMeans(ps_mat))
+    } else {
+        df <- euler(ps_mat_bin)
+    }
+
+    plot(df, quantities = list(type=type), ...)
+
+}

R/ps_refactor.R

@@ -6,6 +6,7 @@
 #' @param psmelted A phyloseq object melted into a data.frame with psmelt
 #' @param ... Arguments passed to hclust
 #' @importFrom reshape2 dcast
+#' @importFrom stats dist hclust
 #' @export
 ps_refactor <- function(psmelted, ...){
     if(!is(psmelted, "data.frame")){

README.md

@@ -12,6 +12,11 @@ MicEco: Various functions for analysis for microbial community data
 ### Citation
 [![DOI](https://zenodo.org/badge/83547545.svg)](https://zenodo.org/badge/latestdoi/83547545)
 
+
+#### phylo_euler
+
+Make Euler (Venn) diagram of shared taxa (ASVs, OTUs) across sample groups from a phyloseq object. Overlap can be weighted by relative abundance
+
 #### adonis_OmegaSq
 
 Calculate the unbiased effect size estimation (partial) omega-squared for adonis (PERMANOVA) models