From f9f62052e4e574afe0af7d3912684969007b1cf6 Mon Sep 17 00:00:00 2001
From: Matteo Vecchi <57488469+tardipede@users.noreply.github.com>
Date: Wed, 16 Feb 2022 16:41:30 +0200
Subject: [PATCH] Update script.R

---
 script.R | 10 ++++++++++
 1 file changed, 10 insertions(+)

diff --git a/script.R b/script.R
index 286a785..a0f7b30 100644
--- a/script.R
+++ b/script.R
@@ -126,3 +126,13 @@ plot_final = (plot_wild|plot_experiment|plot_reproduction) + plot_layout(widths
 ggsave(plot_final,filename = "plot.jpg",width = 40,height = 20, units = "cm", dpi = 300)
 ggsave(plot_final,filename = "plot.svg",width = 40,height = 20, units = "cm", dpi = 300)
 
+
+### ADDITIONS AFTER REVIEW ###
+# After the first review round we added a statistical comparison between different survival rates between tardigrades in wild vs. experiment snails
+# And comparison with the same results, but for miter sor Tukce et al., 2018 (we used data from this paper supplementary material)
+
+prop.test(c(5,120),c(10,172)) # alive/total tardigrades vs. mites in feces of wild animals
+prop.test(c(218,40),c(269,69)) # alive/ingested tardigrades vs. mites in feces of experiment animals
+prop.test(c(5,218),c(10,296)) # alive/total expelled tardigrades in wild vs. lab snails
+
+