add labels to “ggplot :: facet_wrap” chart from seconds table

Hermes · August 6, 2020, 4:14am

How to add the labels at the failure points to a "ggplot :: facet_wrap" graph. Given the need to show these in each of the scenarios graphed The points of interest can be found in the table: TINFAULT_lbl

  Fault          LastFault          Start_Time            End_Time Duration    labels
1: IN-FAULT HIGH-PRESS-PROTECT 2020-08-03 16:41:33 2020-08-03 16:41:36 00:00:03 InFault_1
2: IN-FAULT HIGH-PRESS-PROTECT 2020-08-03 16:41:43 2020-08-03 16:41:45 00:00:02 InFault_2

 And it is of interest to locate them ("Start_Time"),
mark them and label them according to the "labels" column

library(data.table)
library(ggplot2)
DTdata<-data.table(
                 Seq = c(1L,2L,3L,4L,5L,6L,7L,
                         8L,9L,10L,11L,12L,13L,14L,15L,16L,17L,18L,
                         19L,20L,21L,22L,23L,24L,25L,26L,27L,28L,29L,
                         30L,31L,32L,33L,34L,35L,36L),
                Time = c("03/08/2020 16:41:17",
                         "03/08/2020 16:41:18","03/08/2020 16:41:19",
                         "03/08/2020 16:41:20","03/08/2020 16:41:21","03/08/2020 16:41:22",
                         "03/08/2020 16:41:23","03/08/2020 16:41:24",
                         "03/08/2020 16:41:25","03/08/2020 16:41:26",
                         "03/08/2020 16:41:27","03/08/2020 16:41:28","03/08/2020 16:41:29",
                         "03/08/2020 16:41:30","03/08/2020 16:41:31",
                         "03/08/2020 16:41:32","03/08/2020 16:41:33","03/08/2020 16:41:34",
                         "03/08/2020 16:41:35","03/08/2020 16:41:36",
                         "03/08/2020 16:41:37","03/08/2020 16:41:38",
                         "03/08/2020 16:41:39","03/08/2020 16:41:40","03/08/2020 16:41:41",
                         "03/08/2020 16:41:42","03/08/2020 16:41:43",
                         "03/08/2020 16:41:44","03/08/2020 16:41:45",
                         "03/08/2020 16:41:46","03/08/2020 16:41:47","03/08/2020 16:41:48",
                         "03/08/2020 16:41:49","03/08/2020 16:41:50",
                         "03/08/2020 16:41:51","03/08/2020 16:41:52"),
                Temp = c(31L,31L,31L,31L,31L,
                         31L,33L,33L,33L,33L,33L,33L,33L,34L,34L,34L,
                         34L,34L,34L,35L,35L,35L,35L,35L,35L,36L,36L,
                         36L,36L,36L,36L,37L,37L,37L,37L,37L),
               Press = c(18.32,18.32,18.32,18.32,
                         18.32,18.32,19.32,19.32,19.32,19.32,19.32,19.32,
                         19.32,19.84,19.84,19.84,19.84,19.84,19.84,
                         20.36,20.36,20.36,20.36,20.36,20.36,20.89,20.89,
                         20.89,20.89,20.89,20.89,21.44,21.44,21.44,21.44,
                         21.44),
                ICT = c(30L,31L,31L,31L,31L,
                         31L,31L,31L,32L,32L,32L,32L,32L,32L,33L,33L,
                         33L,33L,33L,33L,34L,34L,34L,34L,34L,34L,35L,
                         35L,35L,35L,35L,35L,36L,36L,36L,36L),
              DCCurr = c(74L,74L,74L,74L,74L,
                         76L,76L,76L,76L,76L,76L,76L,76L,76L,76L,76L,
                         100L,100L,100L,100L,76L,76L,76L,76L,86L,86L,86L,
                         101L,101L,101L,86L,86L,86L,86L,90L,90L),
               Fault = c("NO-FAULT","NO-FAULT",
                         "NO-FAULT","NO-FAULT","NO-FAULT","NO-FAULT","NO-FAULT",
                         "NO-FAULT","NO-FAULT","NO-FAULT","NO-FAULT",
                         "NO-FAULT","NO-FAULT","NO-FAULT","NO-FAULT","NO-FAULT",
                         "IN-FAULT","IN-FAULT","IN-FAULT","IN-FAULT",
                         "NO-FAULT","NO-FAULT","NO-FAULT","NO-FAULT","NO-FAULT",
                         "NO-FAULT","IN-FAULT","IN-FAULT","IN-FAULT",
                         "NO-FAULT","NO-FAULT","NO-FAULT","NO-FAULT","NO-FAULT",
                         "NO-FAULT","NO-FAULT"),
           LastFault = c("N/A","N/A","N/A","N/A",
                         "N/A","N/A","N/A","N/A","N/A","N/A","N/A","N/A",
                         "N/A","N/A","N/A","N/A","HIGH-PRESS-PROTECT",
                         "HIGH-PRESS-PROTECT","HIGH-PRESS-PROTECT",
                         "HIGH-PRESS-PROTECT","N/A","N/A","N/A","N/A","N/A","N/A",
                         "HIGH-PRESS-PROTECT","HIGH-PRESS-PROTECT",
                         "HIGH-PRESS-PROTECT","N/A","N/A","N/A","N/A","N/A","N/A","N/A")
        )
DTdata[,Time:=as.POSIXct(Time,tz="EST","%d/%m/%Y %H:%M:%S"),]
print(DTdata,topn=3)
#
DTgpl <- melt(DTdata[,.SD,.SDcols=is.numeric],id.vars=c(1))
#> *Warning in melt.data.table(DTdata[, .SD, .SDcols = is.numeric], id.vars = c(1)):*
#> *'measure.vars' [Temp, Press, ICT, DCCurr, ...] are not all of the same type. By*
#> *order of hierarchy, the molten data value column will be of type 'double'. All*
#> *measure variables not of type 'double' will be coerced too. Check DETAILS in ?*
#> *melt.data.table for more on coercion.*
DTgpl=DTgpl[DTdata[,1:2],on="Seq"]
setcolorder(DTgpl,c(1,4,2,3))
#plot
min <- DTgpl[1,Time] 
max <- DTgpl[.N,Time]
DTgpl[,ggplot(.SD, aes(x = Time,y = value,group = variable,colour=variable)) + 
        geom_line() + 
        theme(legend.position="none") +theme(text = element_text(size=9))+
        facet_wrap( ~ variable, ncol = 2,scales = "free")
      # Set axis limits c(min, max)
      + scale_x_datetime(limits = c(min, max),date_labels = "%b/%d %H:%M")
      ,]

###
#FAULTS
ODUcl<-c("Seq","Time","Fault","LastFault" )
DTFaults<-na.omit(DTdata[,c("Seq","Time","Fault","LastFault" )])
DTFaults[,grpF:= rleid(Fault=="IN-FAULT"),]
ptsFaults<-DTFaults[,
                        .(LastFault=head(LastFault,1),
                          Start_Time =head(Time,1),
                          End_Time =tail(Time, 1))
                        ,by=.(grpF,Fault)][,
                        Duration:=hms::as_hms(difftime(End_Time , Start_Time,"%H:%M:%S"))][,grpF:=NULL]
#table Faults
##labels IN FAULT
TINFAULT=ptsFaults[which(Fault =="IN-FAULT"),]
##labels IN FAULT
TINFAULT_lbl<-TINFAULT[,labels:=rleid(Start_Time , prefix="InFault_")]

Hermes · August 9, 2020, 6:56am

I have managed to mark the start of the failure, by introducing these changes:

p<-DTgpl[,ggplot(.SD, aes(x = Time,y = value,group = variable,colour=variable)) + 
        geom_line() + 
        theme(legend.position="none") +theme(text = element_text(size=9))+
        facet_wrap( ~ variable, ncol = 2,scales = "free")
      # Set axis limits c(min, max)
      + scale_x_datetime(limits = c(min, max),date_labels = "%b/%d %H:%M")
      ,]

TINFAULT_lbl=TINFAULT_lbl[DTgpl,on=.(Start_Time=Time),nomatch=NULL][,.(Start_Time,variable,value,labels)]

p + geom_point(data = TINFAULT_lbl, aes(x = Start_Time, y = value, colour = variable), size=2)

Rplot ggplot_labels2
But I get this alert when I try to show the labels:

p + geom_point(data = TINFAULT_lbl, aes(x = Start_Time, y = value, colour = variable), size=2)
> p + geom_point(data = TINFAULT_lbl, aes(x = Start_Time, y = value, colour = variable), size=2)+
+     geom_text(data = TINFAULT_lbl,aes(x = Start_Time, y = value, colour = variable),label=labels)
Error in rep(value[[k]], length.out = n) : 
  attempt to replicate an object of type 'closure'

How to proceed to display the labels?

Hermes · August 9, 2020, 11:37am

p + geom_point(data = TINFAULT_lbl, aes(x = Start_Time, y = value, colour = variable), size=2)+
    geom_text(data = TINFAULT_lbl,aes(x = Start_Time, y = value*1.01,label=labels),color = "red",
              size = 3,check_overlap = T)

Rplot ggplot_labels3

system · August 16, 2020, 11:37am

This topic was automatically closed 7 days after the last reply. New replies are no longer allowed.