Add Krona plot

inbo · Feb 12, 2024 · 535f195 · 535f195
1 parent 879049f
commit 535f195
Showing 1 changed file with 15 additions and 0 deletions.
diff --git a/source/rmarkdown/annelida_data_analyse.Rmd b/source/rmarkdown/annelida_data_analyse.Rmd
@@ -236,6 +236,21 @@ tidy_physeq_rarefied %>%
   tidytacos::tacoplot_stack()
 ```
 
+```{r}
+# Create a new combined column
+combined_column <- paste(sample_data(physeq_rarefied)$Landgebruik_MBAG, sample_data(physeq_rarefied)$Diepte, sep = "_")
+
+# Add the new combined column to the sample metadata with the desired name
+sample_data(physeq_rarefied)$Landgebruik_MBAG_diepte <- combined_column
+
+# Update the phyloseq object with the modified sample data
+physeq_rarefied <- merge_phyloseq(physeq_rarefied, sample_data(physeq_rarefied))
+
+
+psadd::plot_krona(physeq_rarefied, "MBAG_Olig01_Annelida_rar_species_alle_stalen_Landgebruik_MBAG_per_diepte", "Landgebruik_MBAG_diepte", trim = T)
+
+```
+
 ```{r}
 tidy_physeq_rarefied %>%
   tidytacos::tacoplot_stack(x = Landgebruik_MBAG)