server.R 134 KB
Newer Older
stevenn's avatar
stevenn committed
1
shinyServer(function(input, output,session) {
2
  
svolant's avatar
svolant committed
3
4
  hide(id = "loading-content", anim = TRUE, animType = "fade",time=1.5)
  hide(id = "loading-content-bar", anim = TRUE, animType = "fade",time=1.5)
5
6
7
8
9
  #####################################################
  ##
  ##                    LOAD FILES
  ##
  #####################################################
stevenn's avatar
stevenn committed
10
  
11
12
  
  
stevenn's avatar
stevenn committed
13
14
  ## Create base for contrast
  rand = floor(runif(1,0,1e9))
Amine  GHOZLANE's avatar
Amine GHOZLANE committed
15
  namesfile = tempfile(pattern = "BaseContrast", tmpdir = tempdir(), fileext = "")
stevenn's avatar
stevenn committed
16
  file.create(namesfile,showWarnings=FALSE)
17
  target = NULL
18
  
svolant's avatar
svolant committed
19
20
21
22
  
  ## JSON name for masque
  curdir  = getwd()
  json_name = tempfile(pattern = "file", tmpdir = paste(curdir,"www","masque","todo",sep= .Platform$file.sep),  fileext = ".json")
23
  
svolant's avatar
svolant committed
24
25
26
  ## Pass for MASQUE
  pass = gsub("file","",basename(file_path_sans_ext(json_name)))
  
svolant's avatar
svolant committed
27
  ## Popup messages
svolant's avatar
svolant committed
28
  observe(if(input$AddRegScatter) info("By adding the regression line, you will lose interactivity."))
29
  
30
  ## Reactive target
Amine  GHOZLANE's avatar
Amine GHOZLANE committed
31
  values <- reactiveValues(TargetWorking = target,labeled=NULL,fastq_names_only=NULL,R1fastQ=NULL,R2fastQ=NULL,
svolant's avatar
svolant committed
32
                           json_name=json_name,num=0,pass=pass,login_email = NULL,is.valid =NULL,
Amine  GHOZLANE's avatar
Amine GHOZLANE committed
33
34
35
                           biom_masque = NULL,tree_masque=NULL, masque_key = NULL, count_table_masque = NULL, 
                           rdp_annot_masque = NULL, rdp_thres_masque = NULL,
                           paths_fastq_tmp=NULL,curdir=curdir, error_progress=FALSE)
36
  
stevenn's avatar
stevenn committed
37
  ## Counts file
stevenn's avatar
stevenn committed
38
  dataInputCounts <-reactive({ 
39
    
svolant's avatar
svolant committed
40
    data = NULL
stevenn's avatar
stevenn committed
41
    inFile <- input$fileCounts
Amine  GHOZLANE's avatar
Amine GHOZLANE committed
42
43
    if (is.null(inFile) && is.null(values$count_table_masque)) return(NULL)
    #if (is.null(inFile)) return(NULL)
44
    
Amine  GHOZLANE's avatar
Amine GHOZLANE committed
45
46
47
48
49
50
51
52
53
54
55
56
57
    if (!is.null(values$count_table_masque) && file.exists(values$count_table_masque)){
      tryCatch(read.csv(values$count_table_masque,sep="\t",header=TRUE,check.names=FALSE)->data,
               error=function(e) sendSweetAlert(messageId="ErrorCounts",
                                                title = "Oops",
                                                text=paste("The count file can not be read in SHAMAN.\n \n",e),type ="error"))
    }
    else{
      tryCatch(read.csv(inFile$datapath,sep=input$sepcount,header=TRUE,check.names=FALSE)->data,
               error=function(e) sendSweetAlert(messageId="ErrorCounts",
                                                title = "Oops",
                                                text=paste("Your file can not be read in SHAMAN.\n \n",e),type ="error"))
    }
    #print(data)
svolant's avatar
svolant committed
58
    if(!is.null(data)){
Amine  GHOZLANE's avatar
Amine GHOZLANE committed
59
      colnames(data) = gsub("-",".",colnames(data))
svolant's avatar
svolant committed
60
61
      ## Rownames
      if(!TRUE%in%duplicated(data[,1])) rownames(data)=data[,1];data=data[,-1]
Amine  GHOZLANE's avatar
Amine GHOZLANE committed
62
      try(round(data, 0)->data, silent=T)
svolant's avatar
svolant committed
63
    }
64
    
stevenn's avatar
stevenn committed
65
66
67
    return(as.data.frame(data))
  })
  
68
69
  
  
stevenn's avatar
stevenn committed
70
71
  ## Taxo File
  dataInputTaxo <-reactive({ 
stevenn's avatar
stevenn committed
72
    
stevenn's avatar
stevenn committed
73
    inFile <- input$fileTaxo
74
    
Amine  GHOZLANE's avatar
Amine GHOZLANE committed
75
76
    if (is.null(inFile) && is.null(values$rdp_annot_masque)) return(NULL)
    #if (is.null(inFile)) return(NULL)
77
    
Amine  GHOZLANE's avatar
Amine GHOZLANE committed
78
    if(input$TypeTaxo=="Table" && !is.null(inFile)) 
Stevenn Volant's avatar
Stevenn Volant committed
79
    {
svolant's avatar
svolant committed
80
81
82
83
      tryCatch(read.csv(inFile$datapath,sep=input$septaxo,header=TRUE)->data,
               error=function(e) sendSweetAlert(messageId="ErrorTaxo",
                                                title = "Oops",
                                                text=paste("Your file can not be read in SHAMAN.\n \n",e),type ="error"))
84
      
Stevenn Volant's avatar
Stevenn Volant committed
85
      ## Rownames
svolant's avatar
svolant committed
86
87
88
89
90
91
92
93
94
      if(!is.null(data))
      {
        if(!TRUE%in%duplicated(data[,1])){ 
          DataNames=data[,1]
          colNames=colnames(data)[-1]
          data=as.matrix(data[,-1])
          rownames(data)=DataNames
          colnames(data) = colNames
        }
95
      }
Stevenn Volant's avatar
Stevenn Volant committed
96
    }
stevenn's avatar
stevenn committed
97
    
Amine  GHOZLANE's avatar
Amine GHOZLANE committed
98
    if(input$TypeTaxo=="RDP" && !is.null(inFile) || !is.null(values$rdp_annot_masque)) 
Amine  GHOZLANE's avatar
Amine GHOZLANE committed
99
    {
Amine  GHOZLANE's avatar
Amine GHOZLANE committed
100
101
102
103
104
105
106
107
108
109
110
111
      if (!is.null(values$rdp_annot_masque) && file.exists(values$rdp_annot_masque)){
        tryCatch(read_rdp(values$rdp_annot_masque,values$rdp_thres_masque)->data,
                 error=function(e) sendSweetAlert(messageId="ErrorRDP",
                                                  title = "Oops",
                                                  text=paste("The annotation file can not be read in SHAMAN.\n \n",e),type ="error"))
      }
      else{
        tryCatch(read_rdp(inFile$datapath,input$RDP_th)->data,
                 error=function(e) sendSweetAlert(messageId="ErrorRDP",
                                                  title = "Oops",
                                                  text=paste("Your file can not be read in SHAMAN.\n \n",e),type ="error"))
      }
svolant's avatar
svolant committed
112
      
Amine  GHOZLANE's avatar
Amine GHOZLANE committed
113
    }
stevenn's avatar
stevenn committed
114
115
116
    
    ## Add NA
    data=as.matrix(data)
Stevenn Volant's avatar
Stevenn Volant committed
117
    indNa = which(data=="")
stevenn's avatar
stevenn committed
118
119
    data[indNa]=NA
    
stevenn's avatar
stevenn committed
120
    return(as.data.frame(data))
stevenn's avatar
stevenn committed
121
  })
122
123
  
  
stevenn's avatar
stevenn committed
124
125
  ## BIOM File
  dataInputBiom <-reactive({ 
stevenn's avatar
stevenn committed
126
    
svolant's avatar
svolant committed
127
    data = NULL
stevenn's avatar
stevenn committed
128
    inFile <- input$fileBiom
129
    
Amine  GHOZLANE's avatar
Amine GHOZLANE committed
130
    if (!is.null(inFile) && is.null(values$biom_masque)){
svolant's avatar
svolant committed
131
132
133
134
      tryCatch(read_biom(inFile$datapath)->data,
               error=function(e) sendSweetAlert(messageId="ErrorBiom1",
                                                title = "Oops",
                                                text=paste("Your file can not be read in SHAMAN.\n \n",e),type ="error"))
135
136
      
    }
svolant's avatar
svolant committed
137
138
139
140
141
142
    if (!is.null(values$biom_masque) && file.exists(values$biom_masque)){ 
      tryCatch(read_biom(values$biom_masque)->data,
               error=function(e) sendSweetAlert(messageId="ErrorBiom2",
                                                title = "Oops",
                                                text=paste("Your file can not be read in SHAMAN.\n \n",e),type ="error"))
    }
stevenn's avatar
stevenn committed
143
    return(data)
stevenn's avatar
stevenn committed
144
  })
145
  
Amine  GHOZLANE's avatar
Amine GHOZLANE committed
146
147
148
149
150
151
152
  observeEvent(input$fileCounts,{
    values$count_table_masque=NULL;
  })
  observeEvent(input$fileTaxo,{
    values$rdp_annot_masque=NULL;
    values$rdp_thres_masque=NULL;
  })
svolant's avatar
svolant committed
153
154
155
  observeEvent(input$fileBiom,{
    values$biom_masque=NULL;
  })
156
  
157
  
svolant's avatar
svolant committed
158
159
160
161
162
  ## Unifrac File (tree)
  dataInputTree <-reactive({ 
    
    data = NULL
    inFile <- input$fileTree
163
    
svolant's avatar
svolant committed
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
    if (!is.null(inFile) && is.null(values$tree_masque)) {
      try(read.tree(inFile$datapath)->data, silent=T)
      CheckTree = CheckTreeFile(data)
      data = CheckTree$tree
      try(readLines(inFile$datapath)->treeseq, silent=T)
      return(list(data=data, Error=CheckTree$Error, Warning=CheckTree$Warning, treeseq=treeseq))
    }
    
    if (!is.null(values$tree_masque) && file.exists(values$tree_masque)) {
      try(read.tree(values$tree_masque)->data, silent=T)
      CheckTree = CheckTreeFile(data)
      data = CheckTree$tree
      try(readLines(values$tree_masque)->treeseq, silent=T)
      return(list(data=data, Error=CheckTree$Error, Warning=CheckTree$Warning, treeseq=treeseq))
    }
179
    
svolant's avatar
svolant committed
180
181
182
183
184
  })
  
  
  observeEvent(input$fileTree,{
    values$tree_masque=NULL;
svolant's avatar
svolant committed
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
  })
  
  
  # Infobox Tree (Unifrac)
  output$InfoTreePhylo_box <- renderInfoBox({
    input$fileTree
    tree_tmp = isolate(dataInputTree())
    tree = tree_tmp$data
    
    res = infoBox(h6(strong("Phylogenetic tree")), subtitle = h6(strong("Load the phylogenetic tree (optional)")) ,color = "light-blue",width=NULL,fill=TRUE, icon = icon("upload"))
    if(!is.null(tree)){
      if(!is.null(isolate(input$fileTree))){
        res = infoBox(h6(strong("Phylogenetic tree")), subtitle = h6("The phylogenetic has been loaded") ,color = "green",width=NULL,fill=TRUE, icon = icon("thumbs-o-up"))
        if(!is.null(tree_tmp$Warning)){      
          res = infoBox(h6(strong("Phylogenetic tree")), subtitle = h6(tree_tmp$Warning) ,color = "orange",width=NULL,fill=TRUE, icon = icon("warning"))
        }
        if(!is.null(tree_tmp$Error)){      
          res = infoBox(h6(strong("Phylogenetic tree")), subtitle = h6(tree_tmp$Error),color = "red",width=NULL,fill=TRUE, icon = icon("thumbs-o-down"))
        }
      }
    } 
    return(res)
  })
208
  
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
  observe({
    val <- input$annotationKingdomthreshold
    # Control the value, min, max, and step.
    # Step size is 2 when input value is even; 1 when value is odd.
    updateSliderInput(session, "annotationPhylumthreshold", value = input$annotationPhylumthreshold,
                      min = val, max = 1, step = 0.005)
  })
  observe({
    val <- input$annotationPhylumthreshold[2]
    # Control the value, min, max, and step.
    # Step size is 2 when input value is even; 1 when value is odd.
    updateSliderInput(session, "annotationClassthreshold", value = input$annotationClassthreshold,
                      min = val, max = 1, step = 0.005)
  })
  observe({
    val <- input$annotationClassthreshold[2]
    # Control the value, min, max, and step.
    # Step size is 2 when input value is even; 1 when value is odd.
    updateSliderInput(session, "annotationOrderthreshold", value = input$annotationOrderthreshold,
                      min = val, max = 1, step = 0.005)
  })
  observe({
    val <- input$annotationOrderthreshold[2]
    # Control the value, min, max, and step.
    # Step size is 2 when input value is even; 1 when value is odd.
    updateSliderInput(session, "annotationFamilythreshold", value = input$annotationFamilythreshold,
                      min = val, max = 1, step = 0.005)
  })
  observe({
    val <- input$annotationFamilythreshold[2]
    # Control the value, min, max, and step.
    # Step size is 2 when input value is even; 1 when value is odd.
    updateSliderInput(session, "annotationGenusthreshold", value = input$annotationGenusthreshold,
                      min = val, max = 1, step = 0.005)
  })
  observe({
    val <- input$annotationGenusthreshold[2]
    # Control the value, min, max, and step.
    # Step size is 2 when input value is even; 1 when value is odd.
    updateSliderInput(session, "annotationSpeciethreshold", value = input$annotationSpeciethreshold,
                      min = val, max = 1, step = 0.005)
  })
stevenn's avatar
stevenn committed
251
252
  ## Input data
  dataInput <-reactive({ 
stevenn's avatar
stevenn committed
253
    
stevenn's avatar
stevenn committed
254
    data = NULL
Stevenn Volant's avatar
Stevenn Volant committed
255
256
    check = NULL
    percent = NULL
svolant's avatar
svolant committed
257
258
    Taxo = NULL
    Counts = NULL
svolant's avatar
svolant committed
259
260
    inputData = NULL
    
stevenn's avatar
stevenn committed
261
    if(input$FileFormat=="fileCounts")
stevenn's avatar
stevenn committed
262
    {
stevenn's avatar
stevenn committed
263
      Counts = dataInputCounts()
svolant's avatar
svolant committed
264
265
266
      if(!input$NoTaxoFile) Taxo = dataInputTaxo()
      if(!is.null(Counts) && input$NoTaxoFile) {Taxo = data.frame(rownames(Counts),row.names = rownames(Counts));names(Taxo)=NA}
      
Stevenn Volant's avatar
Stevenn Volant committed
267
268
      if(!is.null(Counts) && !is.null(Taxo))
      { 
Amine  GHOZLANE's avatar
Amine GHOZLANE committed
269
        tmp = GetDataFromCT(Counts,Taxo, ifelse(input$TypeTable=="MGS" && input$FileFormat!="fileBiom", TRUE, FALSE))
Stevenn Volant's avatar
Stevenn Volant committed
270
        data = list(counts=tmp$counts,taxo=tmp$taxo)
svolant's avatar
svolant committed
271
272
273
        ## Remove row with only O
        # data[["counts"]] = data[["counts"]][rowSums(data[["counts"]])>1,]
        
Stevenn Volant's avatar
Stevenn Volant committed
274
275
276
        check = list(CheckCounts=tmp$CheckCounts,CheckTaxo=tmp$CheckTaxo,CheckPercent=tmp$CheckPercent)
        percent = tmp$Percent
      }    
stevenn's avatar
stevenn committed
277
278
    }
    
stevenn's avatar
stevenn committed
279
    if(input$FileFormat=="fileBiom")
stevenn's avatar
stevenn committed
280
    {
stevenn's avatar
stevenn committed
281
      tmpBIOM = dataInputBiom()
Stevenn Volant's avatar
Stevenn Volant committed
282
283
284
285
      if(!is.null(tmpBIOM))
      {
        tmp = GetDataFromBIOM(tmpBIOM)
        data = list(counts=tmp$counts,taxo=tmp$taxo)
svolant's avatar
svolant committed
286
287
288
        ## Remove row with only O
        # data[["counts"]] = data[["counts"]][rowSums(data[["counts"]])>1,]
        
Stevenn Volant's avatar
Stevenn Volant committed
289
290
291
        check = list(CheckCounts=tmp$CheckCounts,CheckTaxo=tmp$CheckTaxo,CheckPercent=tmp$CheckPercent)
        percent = tmp$Percent
      }
stevenn's avatar
stevenn committed
292
293
    }
    
svolant's avatar
svolant committed
294
    
295
296
297
298
299
300
301
302
303
304
305
    
    #     if(input$FileFormat=="fileRData")
    #     {
    #       inFile <- input$fileRData
    #       load(inFile)
    #       if(!is.null(inputData)){
    #         data = inputData$data
    #         check = inputData$check
    #         percent = inputData$percent
    #       }
    #     }
svolant's avatar
svolant committed
306
    
Stevenn Volant's avatar
Stevenn Volant committed
307
    return(list(data=data,check=check,percent=percent))
stevenn's avatar
stevenn committed
308
309
310
  })
  
  
311
312
313
314
315
316
317
  ## Size factor file (optional)
  dataSizeFactors <-reactive({ 
    
    inFile <- input$fileSizeFactors
    
    if (is.null(inFile)) return(NULL)
    
svolant's avatar
svolant committed
318
319
320
321
    tryCatch(read.csv(inFile$datapath,sep=input$sepsize,header=TRUE)->data,
             error=function(e) sendSweetAlert(messageId="ErrorSizeFactor",
                                              title = "Oops",
                                              text=paste("Your file can not be read in SHAMAN.\n \n",e),type ="error"))
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
    return(as.data.frame(data))
  })
  
  
  ## Size factor file (optional)
  SizeFactors_fromFile <-reactive({ 
    
    Error = NULL
    Check = TRUE
    
    data = dataSizeFactors()
    normFactors = dataMergeCounts()$normFactors
    
    if(!is.null(data)){
      ## Check the format
      
      tmp = as.numeric(data)
      names(tmp) = colnames(data)
      
      if(length(tmp)!=length(normFactors)){Error = "The number of samples is not the same than in the target file, size factors will be estimated"; Check = FALSE}
      if(!identical(names(tmp),names(normFactors))){Error = "The names are not the same or in the same order than in the target file, size factors will be estimated"; Check = FALSE}
343
      
344
345
346
347
348
349
350
351
352
353
354
      if(Check) normFactors = tmp
    }
    
    return(list(Check = Check,Error = Error,normFactors=normFactors))
  })
  
  
  
  
  
  
Amine  GHOZLANE's avatar
Amine GHOZLANE committed
355
  ## Merge counts data
356
  dataMergeCounts <-reactive({
Stevenn Volant's avatar
Stevenn Volant committed
357
    input$RunDESeq
Amine  GHOZLANE's avatar
Amine GHOZLANE committed
358
359
360
    
    counts = NULL
    CheckTarget = FALSE
361
    normFactors = NULL
Stevenn Volant's avatar
Stevenn Volant committed
362
    CT_noNorm = NULL
363
    CT_Norm = NULL
364
    ChTM = NULL
365
    ChMC = NULL
366
    data = isolate(dataInput()$data)
367
368
    target = isolate(values$TargetWorking)
    labeled= isolate(values$labeled)
369
    taxo = isolate(input$TaxoSelect)
370
    withProgress(
371
      if(!is.null(data$counts) && !is.null(data$taxo) && nrow(data$counts)>0 && nrow(data$taxo)>0 && !is.null(taxo) && taxo!="..." && !is.null(target)) 
372
      {
373
        design = GetDesign(isolate(input),target)
374
        print(design)
375
376
        ChTM = CheckTargetModel(input,target,labeled,data$counts)$Error
        if(!is.null(design) && is.null(ChTM))
377
        {
378
          tmp = isolate(GetCountsMerge(input,data,taxo,target,design))
379
380
381
          #ChMC = tmp$Error
          #if (!is.null(ChMC))
          #{
382
383
384
385
386
387
388
389
390
391
392
393
394
395
            counts = tmp$counts
            ## Filtering the counts
            if(isolate(input$AddFilter) && !is.null(isolate(input$SliderThSamp)) && !is.null(isolate(input$SliderThAb)))
            {
              ind.filter =Filtered_feature(counts,isolate(input$SliderThSamp),isolate(input$SliderThAb))$ind
              counts = counts[-ind.filter,]
            }
            CheckTarget = tmp$CheckTarget
            #target = tmp$target
            #labeled = tmp$labeled
            normFactors = tmp$normFactors
            ## OTU table, norm and no norm
            CT_noNorm = tmp$CT_noNorm
            CT_Norm = tmp$CT_Norm
396
          #}
397
        }
398
      }
399
400
      ,message="Merging the counts ...")
    return(list(counts=counts,CheckTarget=CheckTarget,normFactors=normFactors,CT_noNorm=CT_noNorm, CT_Norm=CT_Norm, Error = ChMC))
401
    #return(list(counts=counts,target=target,labeled=labeled,normFactors=normFactors,CT_noNorm=CT_noNorm))
Amine  GHOZLANE's avatar
Amine GHOZLANE committed
402
  })
403
404
  
  
Stevenn Volant's avatar
Stevenn Volant committed
405
406
407
408
409
410
411
  # Infobox Error counts
  output$InfoErrorCounts <- renderInfoBox({
    
    tmp = dataInput()
    data = tmp$data
    check = tmp$check
    cond = (!is.null(data$counts) && nrow(data$counts)>0)
412
    res =infoBox(h6(strong("Count table")), subtitle = h6("Load the count table") ,color = "light-blue",width=NULL,fill=TRUE, icon = icon("upload"))
Stevenn Volant's avatar
Stevenn Volant committed
413
414
415
    
    if(cond)
    {
416
417
418
      if(!is.null(check$CheckCounts$Warning)) res = infoBox(h6(strong("Count table")), subtitle = h6(check$CheckCounts$Warning), icon = icon("warning"),color = "orange",width=NULL,fill=TRUE)
      if(!is.null(check$CheckCounts$Error)) res = infoBox(h6(strong("Count table")), subtitle = h6(check$CheckCounts$Error), icon = icon("thumbs-o-down"),color = "red",width=NULL,fill=TRUE)
      if(is.null(check$CheckCounts$Error) && is.null(check$CheckCounts$Warning)) res = infoBox(h6(strong("Count table")), subtitle = h6(paste("Format of the count table seems to be OK")), icon = icon("thumbs-o-up"),color = "green",width=NULL,fill=TRUE)
Stevenn Volant's avatar
Stevenn Volant committed
419
420
421
422
    }
    
    return(res)
  })
423
  
Stevenn Volant's avatar
Stevenn Volant committed
424
425
426
427
428
429
430
  # Infobox Error counts
  output$InfoErrorTaxo <- renderInfoBox({
    
    tmp = dataInput()
    data = tmp$data
    check = tmp$check
    cond = (!is.null(data$taxo) && nrow(data$taxo)>0)
svolant's avatar
svolant committed
431
    res = infoBox(h6(strong("Taxonomy table")), subtitle = h6("Load the taxonomy table") ,color = "light-blue",width=NULL,fill=TRUE, icon = icon("upload"))
Stevenn Volant's avatar
Stevenn Volant committed
432
433
434
435
436
437
438
439
    
    if(cond)
    {
      if(!is.null(check$CheckTaxo$Warning)) res = infoBox(h6(strong("Taxonomy table")), subtitle = h6(check$CheckTaxo$Warning), icon = icon("warning"),color = "orange",width=NULL,fill=TRUE)
      if(!is.null(check$CheckTaxo$Error)) res = infoBox(h6(strong("Taxonomy table")), subtitle = h6(check$CheckTaxo$Error), icon = icon("thumbs-o-down"),color = "red",width=NULL,fill=TRUE)
      if(is.null(check$CheckTaxo$Error) && is.null(check$CheckTaxo$Warning)) res = infoBox(h6(strong("Taxonomy table")), subtitle = h6(paste("Format of the taxonomy table seems to be OK")), icon = icon("thumbs-o-up"),color = "green",width=NULL,fill=TRUE)
    }
    
svolant's avatar
svolant committed
440
    if(input$NoTaxoFile && input$FileFormat=="fileCounts") res = infoBox(h6(strong("Taxonomy table")), subtitle = h6("No taxonomy table has been uploaded, the analysis can only be done at the OTU/gene level"), icon = icon("warning"),color = "orange",width=NULL,fill=TRUE)
Stevenn Volant's avatar
Stevenn Volant committed
441
442
    return(res)
  })
443
444
  
  
Stevenn Volant's avatar
Stevenn Volant committed
445
446
447
448
449
450
451
  # Infobox Error counts
  output$valueErrorPercent <- renderInfoBox({
    
    tmp = dataInput()
    data = tmp$data
    check = tmp$check
    cond = (!is.null(data$counts) && nrow(data$counts)>0 && !is.null(data$taxo) && nrow(data$taxo)>0)
452
    res = shinydashboard::valueBox(paste0(0, "%"),h6(strong("Annotated features")), color = "light-blue",width=NULL,icon = icon("list"))
453
    
Stevenn Volant's avatar
Stevenn Volant committed
454
455
456
    if(cond)
    {
      percent = round(100*tmp$percent,2)
457
458
      if(percent==0) res = shinydashboard::valueBox(paste0(percent, "%"),h6(strong("Annotated features")), color = "red",width=NULL,icon = icon("list"))  
      if(percent!=0) res = shinydashboard::valueBox(paste0(percent, "%"),h6(strong("Annotated features")), color = "green",width=NULL,icon = icon("list"))  
Stevenn Volant's avatar
Stevenn Volant committed
459
460
461
462
463
      
    }
    
    return(res)
  })
stevenn's avatar
stevenn committed
464
  
465
  
466
467
468
469
470
471
472
473
  ####### Filtering the counts (sliders)
  
  output$ThAb <- renderUI({
    input$AddFilter
    
    res = NULL
    counts = isolate(dataMergeCounts()$counts)
    tot = rowSums(counts)
474
    #save(counts,tot,file="testFilter.RData")
475
    withProgress({tmp = SelectThreshAb(counts,lambda=max(round(sum(counts)/nrow(counts)*0.05),min(tot)+1),graph=FALSE)},message="Loading...")
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
    
    res = sliderInput("SliderThAb","Threshold on the total abundance (in log)",min=0,max=round(max(log(tot+1)),1),value = log(tmp+1))
    return(res)
  })
  
  
  output$ThSamp <- renderUI({
    input$AddFilter
    
    res = NULL
    counts = isolate(dataMergeCounts()$counts)
    counts.bin = as.matrix(counts)
    counts.bin[which(counts>0)] = 1
    nbSampByFeat = rowSums(counts.bin)
    
    ## Default value
    val = round(max(nbSampByFeat)*0.2)
    
    res = sliderInput("SliderThSamp","Threshold on the minimal number of samples",min=0,max=max(nbSampByFeat),value = val)
    return(res)
  })
  
  
  ## Plot for the filtering step$
  
  # plot_filter(counts,th.samp,th.abund,type="Scatter")
  
503
  output$Plot_ThAb <- renderPlot({
504
505
506
507
508
    counts = dataMergeCounts()$counts
    ## output of plot_filter is ggplot class
    plot_filter(counts,input$SliderThSamp,input$SliderThAb,type="Abundance")
    
  })
509
  
510
  output$Plot_ThSamp <- renderPlot({
511
512
513
514
    counts = dataMergeCounts()$counts
    ## output of plot_filter is ggplot class
    plot_filter(counts,input$SliderThSamp,input$SliderThAb,type="Samples")
  })
515
  
516
517
518
519
520
  output$Plot_Scatter_Filter <- renderScatterD3({
    counts = dataMergeCounts()$counts
    ## output of plot_filter is ggplot class
    plot_filter(counts,input$SliderThSamp,input$SliderThAb,type="Scatter")
  })
521
522
523
524
525
  #####################################################
  ##
  ##                DYNAMIC MENU
  ##
  #####################################################
stevenn's avatar
stevenn committed
526
527
528
529
530
  
  
  
  output$dymMenu <- renderMenu({
    
Stevenn Volant's avatar
Stevenn Volant committed
531
532
533
    tmp = dataInput()
    data = tmp$data
    check = tmp$check
stevenn's avatar
stevenn committed
534
    
Stevenn Volant's avatar
Stevenn Volant committed
535
536
537
    ## Check error in the counts and taxonomy table 
    CheckOK = (is.null(check$CheckCounts$Error) && is.null(check$CheckTaxo$Error)  && is.null(check$CheckPercent))
    if(!is.null(data$counts) && !is.null(data$taxo) && nrow(data$counts)>0 && nrow(data$taxo)>0 && CheckOK)
Amine  GHOZLANE's avatar
Amine GHOZLANE committed
538
    {
539
      
svolant's avatar
svolant committed
540
      sidebarMenu(id = "side",
541
542
543
544
545
546
547
548
549
550
551
                  menuItem("Statistical analysis",
                           menuSubItem("Run differential analysis",tabName="RunDiff"),
                           menuSubItem("Diagnostic plots",tabName="DiagPlotTab"),
                           menuSubItem("Tables",tabName="TableDiff"),
                           icon = icon("bar-chart-o"), tabName = "AnaStat"
                  ),
                  menuItem("Visualization",icon = icon("area-chart"),
                           menuSubItem("Global views",tabName="GlobVisu"),
                           menuSubItem("Comparison plots",tabName="CompPlot"),
                           tabName = "Visu")
                  #menuItem("Perspective plots", icon = icon("pie-chart"), tabName = "Krona")
stevenn's avatar
stevenn committed
552
      )
svolant's avatar
svolant committed
553
    } else{ sidebarMenu(id = "side",NULL)}
554
    
stevenn's avatar
stevenn committed
555
556
557
  })
  
  
558
  
stevenn's avatar
stevenn committed
559
560
561
562
563
564
565
  #####################################################
  ##
  ##                DATA TABLE
  ##
  #####################################################
  
  ## Counts Table
566
  output$DataCounts <- DT::renderDataTable(
Stevenn Volant's avatar
Stevenn Volant committed
567
    dataInput()$data$counts, 
stevenn's avatar
stevenn committed
568
    options = list(lengthMenu = list(c(10, 50, -1), c('10', '50', 'All')),
569
                   pageLength = 10,scrollX=TRUE, processing=FALSE
stevenn's avatar
stevenn committed
570
571
    ))
  
572
  ## Counts Table
573
574
575
576
577
  output$DataVenn<- DT::renderDataTable(#{
    #SelContrast = input$ContrastList_table_FC
    #resDiff = ResDiffAnal()
    #BaseContrast = read.table(namesfile,header=TRUE)
    GetData_venn(input,input$ContrastList_table_FC,read.table(namesfile,header=TRUE),ResDiffAnal())$df.tot,
578
579
580
581
    #}
    options = list(lengthMenu = list(c(10, 50, -1), c('10', '50', 'All')),
                   pageLength = 10,scrollX=TRUE, processing=FALSE
    ))
582
583
  
  
stevenn's avatar
stevenn committed
584
  ## Taxonomy table
585
  output$DataTaxo <- DT::renderDataTable(
Stevenn Volant's avatar
Stevenn Volant committed
586
    dataInput()$data$taxo, 
stevenn's avatar
stevenn committed
587
    options = list(lengthMenu = list(c(10, 50, -1), c('10', '50', 'All')),
svolant's avatar
svolant committed
588
                   pageLength = 10,scrollX=TRUE, processing=FALSE
stevenn's avatar
stevenn committed
589
    ))
590
  
stevenn's avatar
stevenn committed
591
  
stevenn's avatar
stevenn committed
592
  ## Tab box for data visualisation
stevenn's avatar
stevenn committed
593
  output$TabBoxData <- renderUI({
594
    tree = dataInputTree()$data
Stevenn Volant's avatar
Stevenn Volant committed
595
    data=dataInput()$data
svolant's avatar
svolant committed
596
597
    res=NULL
    if(!is.null(tree))
598
    {
599
600
601
602
603
604
605
606
      res = tabBox(width = NULL, selected = "Count table",
                   tabPanel("Count table",DT::dataTableOutput("DataCounts")),
                   tabPanel("Taxonomy",DT::dataTableOutput("DataTaxo")),
                   tabPanel("Summary",h5(strong("Percentage of annotation")),htmlOutput("SummaryView"),
                            br(),h5(strong("Number of features by level:")),plotOutput("SummaryViewBarplot",width = 1200,height=500)),
                   tabPanel("Phylogeny", PhyloTreeMetaROutput('PhyloTreeMetaR'))
      )
      
607
    }
svolant's avatar
svolant committed
608
    else if(is.null(tree))
stevenn's avatar
stevenn committed
609
    {
610
611
612
613
614
615
      res = tabBox(width = NULL,selected = "Count table",
                   tabPanel("Count table",DT::dataTableOutput("DataCounts")),
                   tabPanel("Taxonomy",DT::dataTableOutput("DataTaxo")),
                   tabPanel("Summary",h5(strong("Percentage of annotation")),htmlOutput("SummaryView"),
                            br(),h5(strong("Number of features by level:")),plotOutput("SummaryViewBarplot",width = 1200,height=500))
      )
svolant's avatar
svolant committed
616
      
stevenn's avatar
stevenn committed
617
    }
svolant's avatar
svolant committed
618
619
620
621
622
623
624
    return(res)
  })
  
  observe({
    data=dataInput()$data
    if(!is.null(data$counts) && !is.null(data$taxo) && nrow(data$counts)>0 && nrow(data$taxo)>0)
    {
625
      showElement("tabboxdata_col",anim=TRUE)
svolant's avatar
svolant committed
626
    } else hideElement("tabboxdata_col",anim=TRUE)
stevenn's avatar
stevenn committed
627
    
stevenn's avatar
stevenn committed
628
  })
629
  
630
631
632
633
  output$PhyloTreeMetaR <- renderPhyloTreeMetaR({
    PhyloTreeMetaR(dataInputTree()$treeseq,NULL)
  })
  
svolant's avatar
svolant committed
634
  output$SummaryView <- renderGvis({
svolant's avatar
svolant committed
635
636
    tmp = dataInput()
    data = tmp$data
svolant's avatar
svolant committed
637
638
    taxo = data$taxo
    counts = data$counts
svolant's avatar
svolant committed
639
640
    check = tmp$check
    cond = (!is.null(data$counts) && nrow(data$counts)>0 && !is.null(data$taxo) && nrow(data$taxo)>0 && is.null(check$CheckTaxo$Error) && is.null(check$CheckCounts$Error))
641
    
svolant's avatar
svolant committed
642
    res = NULL
svolant's avatar
svolant committed
643
    if(cond)
svolant's avatar
svolant committed
644
    {
svolant's avatar
svolant committed
645
      taxo = rbind(taxo,rep(NA,ncol(taxo)))
646
647
648
649
650
651
652
653
654
655
656
657
658
659
      #tmpPercent = round(apply(is.na(taxo),2,table)["FALSE",]/(nrow(taxo)-1)*100,2)
      
      tmp = apply(is.na(taxo),2,table)
      
      if (class(tmp) == "list") {
        tmp2 = sapply(tmp, function (x) {if (! "FALSE" %in% names(x)) {x["FALSE"] = 0} ; return(x["FALSE"])})
      }
      else
      {
        tmp2 = tmp["FALSE",]
      }
      
      tmpPercent = round(tmp2/(nrow(taxo)-1)*100,2)
      
660
      
svolant's avatar
svolant committed
661
      df <- data.frame(Label = colnames(taxo),Value = tmpPercent)
662
      
Amine  GHOZLANE's avatar
Amine GHOZLANE committed
663
664
665
      # res = gvisGauge(df,options=list(min=0, max=100, greenFrom=80,
      #                                 greenTo=100, yellowFrom=60, yellowTo=80,
      #                                 redFrom=0, redTo=60, width=1200, height=300))
svolant's avatar
svolant committed
666
      res = gvisGauge(df,options=list(min=0, max=100, greenFrom=80,
667
                                      greenTo=100, yellowFrom=60, yellowTo=80,
Amine  GHOZLANE's avatar
Amine GHOZLANE committed
668
                                      redFrom=0, redTo=60, width=800, height=200))
svolant's avatar
svolant committed
669
670
671
672
673
674
    }
    return(res)
  })
  
  
  output$SummaryViewBarplot <- renderPlot({
svolant's avatar
svolant committed
675
676
    tmp = dataInput()
    data = tmp$data
svolant's avatar
svolant committed
677
678
    taxo = data$taxo
    counts = data$counts
svolant's avatar
svolant committed
679
680
    check = tmp$check
    cond = (!is.null(data$counts) && nrow(data$counts)>0 && !is.null(data$taxo) && nrow(data$taxo)>0 && is.null(check$CheckTaxo$Error) && is.null(check$CheckCounts$Error))
681
    
svolant's avatar
svolant committed
682
    res = NULL
svolant's avatar
svolant committed
683
    if(cond)
svolant's avatar
svolant committed
684
685
686
687
688
689
690
691
692
693
694
695
696
    {
      colors=rep(c("#1f77b4","#aec7e8","#ff7f0e","#ffbb78", "#2ca02c","#98df8a","#d62728","#ff9896","#9467bd","#c5b0d5","#8c564b",
                   "#c49c94","#e377c2","#f7b6d2","#7f7f7f", "#c7c7c7","#bcbd22","#dbdb8d","#17becf","#9edae5"),ceiling(ncol(taxo)/20))
      tmp = apply(taxo,2,unique)
      nbfeatures = as.numeric(lapply(tmp,length)) -as.numeric(lapply(lapply(tmp,is.na),any))
      df <- data.frame(Label = colnames(taxo),Count = nbfeatures)
      df$Label = factor(df$Label,levels =colnames(taxo) )
      res = ggplot(df,aes(x=Label,y=Count,fill=Label))+geom_bar(stat="identity")
      res = res + theme_bw() + xlab("Taxonomy") + scale_fill_manual(values=colors) + guides(fill=FALSE)
    }
    return(res)
  })
  
stevenn's avatar
stevenn committed
697
698
699
700
701
  #####################################################
  ##
  ##                TARGET FILE
  ##
  #####################################################
stevenn's avatar
stevenn committed
702
  
703
  
stevenn's avatar
stevenn committed
704
  ## Load target file
705
  observe({ 
stevenn's avatar
stevenn committed
706
707
    
    inFile <- input$fileTarget
Stevenn Volant's avatar
Stevenn Volant committed
708
    counts = dataInput()$data$counts
svolant's avatar
svolant committed
709
    labeled = 0
svolant's avatar
svolant committed
710
    data = NULL
711
    
stevenn's avatar
stevenn committed
712
713
    if (is.null(inFile)) return(NULL)
    
svolant's avatar
svolant committed
714
715
    ## Read the data
    try(read.csv(inFile$datapath,sep=input$septarget,header=TRUE)->data,silent=TRUE)
716

svolant's avatar
svolant committed
717
718
719
720
721
722
    if(!is.null(data))
    {
      data = as.data.frame(data)
      names = colnames(data)
      
      ## Change the rownames
svolant's avatar
svolant committed
723
      if(!TRUE%in%duplicated(data[,1])) rownames(data)=gsub(pattern = "-",replacement = ".",as.character(data[,1]))
svolant's avatar
svolant committed
724
725
726
727
728
      
      ## Keep only the row which are in the count table
      ind = which(rownames(data)%in%colnames(counts))
      data = as.data.frame(data[ind,])
      colnames(data) = names
svolant's avatar
svolant committed
729
730
      
      
svolant's avatar
svolant committed
731
732
733
734
735
      ## Replace "-" by "."
      if(ncol(data)>1 && nrow(data)>1){
        ind_num = which(sapply(as.data.frame(data[,-1]),is.numeric)) + 1
        if(length(ind_num)>0){
          data_tmp =cbind( as.data.frame(apply(as.data.frame(data[,-ind_num]),2,gsub,pattern = "-",replacement = ".")),data[,ind_num])
Amine  GHOZLANE's avatar
Amine GHOZLANE committed
736
          #data_tmp =cbind( as.data.frame(as.data.frame(data[,-ind_num])),data[,ind_num])
svolant's avatar
svolant committed
737
738
739
740
          colnames(data_tmp) = c(colnames(data)[-ind_num],colnames(data)[ind_num])
          data = data_tmp
        }
        if(length(ind_num)==0){data = as.data.frame(apply(data,2,gsub,pattern = "-",replacement = "."))}
741
      }
svolant's avatar
svolant committed
742
      
743
      values$TargetWorking = as.data.frame(data)
744
745
      #       ind_sel = Target_selection()
      #       if(length(ind))
svolant's avatar
svolant committed
746
747
748
749
750
751
752
753
754
755
756
      # target = as.data.frame(apply(target,2,gsub,pattern = "-",replacement = "."))
      
      #ord = order(rownames(data))
      #data = data[ord,]
      ### A SUPPRIMER 
      #rownames(data) <- colnames(counts)
      
      # Percent annotated
      #     print(ind)
      #     print(colnames(counts))
      #     print(rownames(data))
757
      values$labeled = length(ind)/length(colnames(counts))*100.0
svolant's avatar
svolant committed
758
    }
759
    
760
761
762
    # return(list(target = target, labeled=labeled))
  })
  
763
  
svolant's avatar
svolant committed
764
765
766
767
768
769
770
771
772
  
  
  #############################################################
  ##
  ##                        MASQUE
  ##
  #############################################################
  
  
svolant's avatar
svolant committed
773
  
svolant's avatar
svolant committed
774
775
776
  observeEvent(input$dir,{
    
    inFiles <- input$dir
Amine  GHOZLANE's avatar
Amine GHOZLANE committed
777
    
svolant's avatar
svolant committed
778
    if (!is.null(inFiles)){
779
780
781
      # values$fastq_names_only = unique(paste(values$fastq_names_only,inFiles$name))
      values$paths_fastq_tmp = rbind(isolate(values$paths_fastq_tmp),inFiles)
      values$fastq_names_only = isolate(unique(values$paths_fastq_tmp[,"name"]))
svolant's avatar
svolant committed
782
    }
svolant's avatar
svolant committed
783
784
  })
  
785
  
svolant's avatar
svolant committed
786
  
svolant's avatar
svolant committed
787
  ## Create a fasta file containing the contaminant
svolant's avatar
svolant committed
788
789
790
791
792
793
794
795
796
797
798
  CreateFasta <- reactive({
    seq = NULL
    tmp = tempdir()
    fastaName = paste(tmp,paste(basename(file_path_sans_ext(json_name)),"_contaminant.fasta",sep=""),sep = .Platform$file.sep)
    
    if(!file.exists(fastaName)) file.create(fastaName,showWarnings=FALSE)
    if(input$PairedOrNot=="y"){seq =paste("#Seq1\n",input$R1primer,"\n \n","#Seq2\n",input$R2primer,sep="")}
    if(input$PairedOrNot=="n"){seq =input$primerSingle}
    if(!is.null(seq))  write(seq, file=fastaName)
    
  })
svolant's avatar
svolant committed
799
  
svolant's avatar
svolant committed
800
801
802
  
  ## Action with submit button
  MasqueSubmit <- eventReactive(input$submit,{
Amine  GHOZLANE's avatar
Amine GHOZLANE committed
803
804
    #activate check_mail
    CMP = CheckMasque(input, values,check_mail = TRUE)
svolant's avatar
svolant committed
805
    Error = CMP$Error
806
    
svolant's avatar
svolant committed
807
    isJSONalreadyExist = file.exists(paste(values$curdir,"www","masque","doing",basename(json_name),sep= .Platform$file.sep))
808
    
svolant's avatar
json    
svolant committed
809
    if(is.null(Error) && !isJSONalreadyExist)
svolant's avatar
svolant committed
810
    {
svolant's avatar
svolant committed
811
812
      CreateFasta()
      values$num = 1
svolant's avatar
svolant committed
813
      tmp = tempdir()
svolant's avatar
svolant committed
814
      # home <- normalizePath("~")
svolant's avatar
svolant committed
815
      home <- ""
svolant's avatar
svolant committed
816
      # path_glob = file.path(home, paste(unlist(dir()$path[-1]), collapse = .Platform$file.sep))
817
818
      
      
svolant's avatar
svolant committed
819
820
821
822
823
824
825
826
827
828
829
      ## Paired-end
      if(input$PairedOrNot=="y"){
        cmp = 0
        nfiles = length(values$R1fastQ)+length(values$R2fastQ)
        
        withProgress(message = 'Uploading files...', value = 0, {
          pathToR1 = paste(tmp,"Masque_files_R1",sep=.Platform$file.sep)
          pathToR2 = paste(tmp,"Masque_files_R2",sep=.Platform$file.sep)
          
          if(dir.exists(pathToR1)){file.remove(list.files(pathToR1,full.names =TRUE))} else dir.create(pathToR1)
          if(dir.exists(pathToR2)){file.remove(list.files(pathToR2,full.names =TRUE))} else dir.create(pathToR2)
830
831
832
833
834
835
          for(i in values$R1fastQ){
            ind=which(i==values$paths_fastq_tmp[,"name"])[1]
            file.rename(from=values$paths_fastq_tmp[,"datapath"][ind], to=paste(tmp,"Masque_files_R1",i,sep= .Platform$file.sep))
            cmp = cmp +1
            incProgress(cmp/nfiles, detail = "Forward fastq files...")
          }
svolant's avatar
svolant committed
836
837
          for(i in values$R2fastQ){
            ind=which(i==values$paths_fastq_tmp[,"name"])[1]
Amine  GHOZLANE's avatar
Amine GHOZLANE committed
838
            file.rename(from=values$paths_fastq_tmp[,"datapath"][ind], to=paste(tmp,"Masque_files_R2",i,sep= .Platform$file.sep))
svolant's avatar
svolant committed
839
840
841
            cmp = cmp +1
            incProgress(cmp/nfiles, detail = "Reverse fastq files...")
          }
svolant's avatar
svolant committed
842
        })
843
        
svolant's avatar
svolant committed
844
      } else{
845
        
svolant's avatar
svolant committed
846
847
848
849
        cmp = 0
        nfiles = length(values$fastq_names_only)
        
        withProgress(message = 'Uploading files...', value = 0, {
850
          
svolant's avatar
svolant committed
851
852
853
854
          pathTo = paste(tmp,"Masque_files",sep=.Platform$file.sep)
          
          if(dir.exists(pathTo)){file.remove(list.files(pathTo,full.names =TRUE))} else dir.create(pathTo)
          
svolant's avatar
svolant committed
855
856
          for(i in values$fastq_names_only){
            ind=which(i==values$paths_fastq_tmp[,"name"])[1]
Amine  GHOZLANE's avatar
Amine GHOZLANE committed
857
            file.rename(from=values$paths_fastq_tmp[,"datapath"][ind], to=paste(tmp,"Masque_files",i,sep= .Platform$file.sep));cmp = cmp +1;incProgress(cmp/nfiles)}
svolant's avatar
svolant committed
858
        })
859
      }
svolant's avatar
svolant committed
860
861
      
      ## Create JSON file
svolant's avatar
svolant committed
862
863
      withProgress(message = 'Creating JSON file...',{CreateJSON(input,values)})
      if(file.exists(values$json_name)) values$num = 1
svolant's avatar
svolant committed
864
865
866
867
868
869
870
871
872
873
      sendSweetAlert(messageId="SuccessMasque",
                     title = "Success",
                     text = paste("Your data have been submitted. You will receive an e-mail once the computation over. <br /> This can take few hours.
                                  <br /> 
                                  <br /> 
                                  <br /> 
                                  <em> Remind: You can close shaman and use your key to check the progression and get your results: </em>",values$pass),
                     type = "success",
                     html=TRUE
                     )
svolant's avatar
svolant committed
874
875
    }
    
svolant's avatar
svolant committed
876
877
878
879
880
881
882
883
884
885
886
  })
  
  
  observeEvent(input$submit,{  
    
    tryCatch(MasqueSubmit(),
             error=function(e) sendSweetAlert(messageId="ErrorMasque",
                                              title = "Oops",
                                              text=paste("Something wrong when submitting.\n \n",e),type ="error"))
    
    
svolant's avatar
svolant committed
887
  },priority = 1)
svolant's avatar
svolant committed
888
  
svolant's avatar
svolant committed
889
  
svolant's avatar
svolant committed
890
891
  ## FastQ list
  output$FastQList_out <- renderUI({
svolant's avatar
svolant committed
892
893
894
895
896
897
898
899
900
901
902
    res = NULL
    if(!is.null(input$dir)){
      NullBox = h3(strong("0 FastQ file detected"),style="color:red;  text-align: center")
      res = NullBox
      
      if(length(values$fastq_names_only)>0)
      {
        res =list(selectInput("FastQList",label = "List of the fastq files in the selected directory",isolate(values$fastq_names_only),multiple =TRUE,selectize=FALSE,size = 6),
                  actionButton("RemoveFastQbut",'Remove file(s)',icon=icon("remove")))
      } else res = NullBox
    }
svolant's avatar
svolant committed
903
904
905
    return(res)
  })
  
906
907
908
  
  
  
svolant's avatar
svolant committed
909
910
911
912
913
914
915
  ## Remove FastQ function
  RemoveFastQ <-eventReactive(input$RemoveFastQbut,{
    
    if(length(input$FastQList)>0)
    {
      ind = which(values$fastq_names_only%in% input$FastQList)
      values$fastq_names_only = values$fastq_names_only[-ind]
svolant's avatar
svolant committed
916
      values$paths_fastq_tmp = values$paths_fastq_tmp[-ind,]
svolant's avatar
svolant committed
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
      updateSelectInput(session, "FastQList","List of the fastq files in the selected directory",values$fastq_names_only)
    }
  })
  
  
  
  
  
  ## Remove FastQ
  observeEvent(input$RemoveFastQbut,{  
    
    RemoveFastQ()
    if(input$MatchFiles_button>=1) MatchFiles()
    
  },priority=1)
  
  
  
  ## Update R1 and R2 lists
  MatchFiles <-reactive({
    
    if(length(values$fastq_names_only)>0 && input$PairedOrNot=='y')
    {
      indR1 = grep(input$R1files,values$fastq_names_only)
      indR2 = grep(input$R2files,values$fastq_names_only)
      
      ## If some are R1 and R2, removed from both list
      b12 = intersect(indR1,indR2)
      if(length(b12)>0) {indR1 = indR1[-which(indR1%in%b12)]; indR2 = indR2[-which(indR2%in%b12)]}
      
      
      if(length(indR1)>0 && length(indR2)>0){
        
        values$R1fastQ = values$fastq_names_only[indR1]
        values$R2fastQ = values$fastq_names_only[indR2]
        
        ## If some are only R1 or R2
        tmpR1 = gsub(input$R1files,x=values$R1fastQ,""); tmpR2 = gsub(input$R2files,x=values$R2fastQ,"")
        values$R1fastQ = values$R1fastQ[tmpR1%in%tmpR2];values$R2fastQ = values$R2fastQ[tmpR2%in%tmpR1]
956
        
svolant's avatar
svolant committed
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
        ## Update the files lists
        updateSelectInput(session, "R1filesList","",values$R1fastQ); updateSelectInput(session, "R2filesList","",values$R2fastQ)
      } else{updateSelectInput(session, "R1filesList","","");updateSelectInput(session, "R2filesList","","")}
    } else{updateSelectInput(session, "R1filesList","","");updateSelectInput(session, "R2filesList","","")}
  })
  
  
  observeEvent(input$MatchFiles_button,{  
    
    MatchFiles()
    
  },priority=1)
  
  
  
  ## Remove FastQ function directly from R1, R2
  RemoveFastQ_R1R2 <-eventReactive(input$RemoveFastQbut_R1R2,{
    
    if(length(input$R1filesList)>0)
    {
      ind = which(values$R1fastQ%in% input$R1filesList)
      values$R1fastQ = values$R1fastQ[-ind]
      updateSelectInput(session, "R1filesList","",values$R1fastQ)
    }
    
    if(length(input$R2filesList)>0)
    {
      ind = which(values$R2fastQ%in% input$R2filesList)
      values$R2fastQ = values$R2fastQ[-ind]
      updateSelectInput(session, "R2filesList","",values$R2fastQ)
987
    }
svolant's avatar
svolant committed
988
989
990
991
    
    
  })
  
svolant's avatar
svolant committed
992
993
994
995
  
  
  RemoveFastQ_R1R2_all <-eventReactive(input$LoadFiles,{
    
996
997
998
999
    values$R1fastQ = NULL
    updateSelectInput(session, "R1filesList","","")
    values$R2fastQ = NULL
    updateSelectInput(session, "R2filesList","","")
svolant's avatar
svolant committed
1000
    
For faster browsing, not all history is shown. View entire blame