diff --git a/source/rmarkdown/annelida_data_analyse.Rmd b/source/rmarkdown/annelida_data_analyse.Rmd
index ea1625f..7e7d927 100644
--- 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)