From c1cd6e92460ae0e7221bc944447250ec0d91c932 Mon Sep 17 00:00:00 2001
From: Bridget McGivern <95941779+bmcgivern13@users.noreply.github.com>
Date: Mon, 25 Mar 2024 15:50:16 -0600
Subject: [PATCH] Update 02_GROWdb_Carbon_Curation.R

added code for CAZy SOM figure
---
 .../02_GROWdb_Carbon_Curation.R               | 26 +++++++++++++++++++
 1 file changed, 26 insertions(+)

diff --git a/USA_SurfaceWater/GenomeBasedCarbonAnalyses/02_GROWdb_Carbon_Curation.R b/USA_SurfaceWater/GenomeBasedCarbonAnalyses/02_GROWdb_Carbon_Curation.R
index cc5399a..fe421f3 100644
--- a/USA_SurfaceWater/GenomeBasedCarbonAnalyses/02_GROWdb_Carbon_Curation.R
+++ b/USA_SurfaceWater/GenomeBasedCarbonAnalyses/02_GROWdb_Carbon_Curation.R
@@ -157,3 +157,29 @@ dunnTest(abund ~ group,
          data=lifestyle,
          method="bh") # headwater and river are significantly different (p.adj=0.001703097)
 
+###########
+## breaking out CAZymes
+###########
+carbon_gene%>%
+  left_join(.,river_info,by=c("code"="Sample"))%>%
+  filter(is.na(StreamUpdated)==0)%>%
+  filter(Guild=="heterotroph-polymer")%>%
+  mutate(filter=ifelse(substr(final_ID,1,3)=="CBM","yes",
+                       ifelse(substr(final_ID,1,1)=="K","yes","no")))%>%
+  mutate(group=ifelse(StreamUpdated>=7,"river",ifelse(StreamUpdated>=4,"midsize","headwater")))%>%
+  filter(filter=="no")%>%
+  group_by(final_ID,code,group)%>%
+  summarise(sum=sum(geTMM))%>%
+  mutate(count=ifelse(sum>1,1/57,0))%>%
+  ggplot()+
+  geom_jitter(aes(x=reorder(final_ID,log(sum+1)),y=-log(sum+1)),size=0.7)+
+  geom_boxplot(aes(x=reorder(final_ID,log(sum+1)),y=-log(sum+1)),fill="NA",color="green",outlier.colour = "NA")+
+  geom_bar(aes(x=reorder(final_ID,log(sum+1)),y=count*100,fill=group),stat="identity",position="stack")+
+  geom_hline(yintercept = 50)+
+  geom_hline(yintercept = -5)+
+  geom_hline(yintercept = -3)+
+  theme_classic()+
+  coord_flip()+
+  ylim(c(-5,60))+
+  xlab("CAZy Family")
+