From 2f676e763261f763c2e2fba0385b2504f63c0a39 Mon Sep 17 00:00:00 2001
From: marcelcosta <marcelcosta@mcosta.cat>
Date: Tue, 20 Sep 2022 17:17:15 +0200
Subject: [PATCH] =?UTF-8?q?Arreglar=20gr=C3=A1ficos=20para=20exportar=20cu?=
 =?UTF-8?q?ando=20es=20de=20vacunaci=C3=B3n.?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 invivos/app.R | 17 ++++++++++++++---
 1 file changed, 14 insertions(+), 3 deletions(-)

diff --git a/invivos/app.R b/invivos/app.R
index 315aabc..9f13173 100755
--- a/invivos/app.R
+++ b/invivos/app.R
@@ -586,10 +586,21 @@ server <- function(input, output) {
       if (input$fig_id %in% c("Cinética Grupo", "Cinética Individual")){
         if (input$fig_id == "Cinética Grupo"){
           if (input$vacc == "Sí"){
+            std<-function(x, na.rm=T){sd(x, na.rm=na.rm)/sqrt(length(x))}
+            errbar<-table %>% group_by(Group,`ID tumor`, Timepoint) %>% 
+              summarise(mean=mean(Volume, na.rm=T), std=std(Volume)) %>% 
+              mutate(Timepoint2=as.numeric(as.character(Timepoint))) %>% 
+              mutate(x=Timepoint2-input$`errorbar-width`, xend=Timepoint2+input$`errorbar-width`)
+            
             g<-ggplot(table, aes(as.numeric(as.character(Timepoint)), Volume, color=Group, group=Group))+
-              scale_x_continuous(expand = expansion(mult = c(0,0.05)), limits = c(0, (round(max(as.numeric(as.character(table$Timepoint))) / 5)+1)*5))+
+              scale_x_continuous(expand = expansion(mult = c(0,0.0)), 
+                                 breaks=sort(unique(errbar$Timepoint2)),
+                                 limits = c(0,max(as.numeric(as.character(table$Timepoint)))*1.1))+
               facet_grid(factor(`ID tumor`, labels = c("Vaccination", "Rechallenge"))~., scale="free_y")+
-              theme_bw()  
+              theme_bw()+
+              geom_segment(data=errbar, aes(y=mean, yend=mean+std, x=Timepoint2,  xend=Timepoint2))+
+              geom_segment(data=errbar, aes(y=mean+std, yend=mean+std,
+                                            x=x,xend=xend))
           }else{
             std<-function(x, na.rm=T){sd(x, na.rm=na.rm)/sqrt(length(x))}
             errbar<-table %>% group_by(Group, Timepoint) %>% 
@@ -612,7 +623,7 @@ server <- function(input, output) {
           if (input$vacc == "Sí"){
             g<-ggplot(table, aes(as.numeric(as.character(Timepoint)), Volume, color=Group, group=`ID animal`))+
               scale_x_continuous(expand = expansion(mult = c(0,0.0)), 
-                                 breaks=sort(unique(errbar$Timepoint2)),
+                                 breaks=sort(unique(as.numeric(as.character(table$Timepoint)))),
                                  limits = c(0, (round(max(as.numeric(as.character(table$Timepoint))) / 5)+1)*5))+
               facet_grid(factor(`ID tumor`, labels = c("Vaccination", "Rechallenge"))~Group, scale="free_y")+
               theme_bw()