display_repeatability.r

library(dplyr)
library(ggplot2)


# File paths to Input and output data
path.rpt          <- ("./rpt_no_public.csv")
path.plot.overall <- ("./rpt_overall.png")
path.plot.sideall <- ("./rpt_sideall.png")





# Load repeatability results
df.rpt <- read.csv(path.rpt)

# Exclude results with NA
df.rpt <- df.rpt[which(!is.na(df.rpt$rpt) & is.finite(df.rpt$rpt)), ]

# Select data with View==SideAll
df.rpt.sideall <- df.rpt[which(df.rpt$View=="SideAll"),]

# Simplify attributes under df.rpt$View
df.rpt$View <- replace(df.rpt$View, 
                       df.rpt$View %in% c("SideAll", "SideAverage"), 
                       "Side")

# Plot results (overall data)
p1 <- ggplot(df.rpt, aes(x=View, y=rpt)) + geom_boxplot() +
       facet_grid(rows=vars(GROWTH_STAGE), cols=vars(TREATMENT)) +
       labs(y = "Repeatability, Input: all RGB and  HSI data")
png(path.plot.overall)
print(p1)
dev.off()

# Plot results (sideall data)
df.rpt.sideall$frame_nr <- as.factor(df.rpt.sideall$frame_nr)
p2 <- ggplot(df.rpt.sideall, aes(x=frame_nr, y=rpt)) + geom_boxplot() +
       facet_grid(rows=vars(GROWTH_STAGE), cols=vars(TREATMENT)) +
       labs(y = "Repeatability, Input: RGB data from different vewing angles")
png(path.plot.sideall)
print(p2)
dev.off()



# EOF