From 1e0c476568b0303a421ef47afc8ede8cdcc87fc4 Mon Sep 17 00:00:00 2001
From: marcelcosta <marcelcosta@protonmail.com>
Date: Tue, 3 Nov 2020 18:41:18 +0100
Subject: [PATCH] =?UTF-8?q?a=C3=B1adir=20randomizado=201?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 invivos/app.R | 48 ++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 48 insertions(+)

diff --git a/invivos/app.R b/invivos/app.R
index 90a1b95..89d0ff1 100755
--- a/invivos/app.R
+++ b/invivos/app.R
@@ -3,6 +3,8 @@ library(ggplot2)
 library(reshape2)
 library(openxlsx)
 library(dplyr)
+library(car)
+library(ggbeeswarm)
 source("../../funcions.R")
 
 # Define UI for application
@@ -25,6 +27,52 @@ ui <- fluidPage(
 
 # Define server logic required to draw a histogram
 server <- function(input, output) {
+  
+  output$distPlot <- renderPlot({
+    df<-read.xlsx(input$file1$datapath)
+    up_cuttof<-400
+    low_cuttof<-50
+    df<-df[df$Volumen < up_cuttof & df$Volumen > low_cuttof,]
+    df["Mouse"]<-gsub("[a-zA-Z]", "", df$MouseID)
+    
+    s<-shapiro.test(df$Volumen)[[2]]
+    
+    ngroup<-5
+    
+    ind.list<-list()
+    pval.list<-list()
+    lvn.list<-list()
+    test.list<-list()
+    for (data in 1:2000){
+      interr=T
+      while(interr == T){
+        ind<-sample(rep(1:ngroup, each=7), length(unique(df$Mouse)))
+        df_temp<-merge(df, data.frame("Mouse"=unique(df$Mouse), "group"=as.factor(ind)))
+        interr<-any(table(df_temp$group) < floor(nrow(df_temp)/5) | table(df_temp$group) > ceiling(nrow(df_temp)/5))
+      }
+      ind.list[[data]]<-df_temp[,c("MouseID","group")]
+      lvn.list[data]<-leveneTest(Volumen ~ group, data = df_temp[,3:4])[[2]][1]
+      if (s < 0.05){
+        k<-kruskal.test(df_temp$Volumen,df_temp$group)
+        test.list[data]<-k[[1]][1]
+        pval.list[data]<-k[[3]][1]
+      }else{
+        res.aov<-aov(Volumen~group, data=df_temp)
+        pval.list[data]<-summary(res.aov)[[1]][[5]][1]
+        test.list[data]<-summary(res.aov)[[1]][[4]][1]
+      }
+    }
+    
+    index<-which(unlist(lvn.list) == min(unlist(lvn.list)[which(unlist(pval.list) %in% sort(unlist(pval.list), decreasing = T)[1:20])]))
+    df_def<-merge(df, ind.list[[index]])
+    
+    ggplot(df_def, aes(group, Volumen))+
+      geom_boxplot(outlier.alpha = F)+
+      geom_jitter(width=0.25)+
+      geom_point(stat="summary", color="blue", size=3)+
+      lims(y=c(0,max(df_def$Volumen)+10))
+  })
+  
   output$downloadData <- downloadHandler(
     
     filename = function() {