panel.R

# Usefull library
library(ggplot2)
library(scales)
library(qpdf)
library(gridExtra)
library(gridtext)
library(dplyr)
library(grid)
library(ggh4x)
library(RColorBrewer)


time_panel = function (df_data_code, df_trend_code, type, p_threshold=0.1, missRect=FALSE, unit2day=365.25, period=NULL, last=FALSE, color=NULL) {


    if (type == 'sqrt(Q)') {
        df_data_code$Qm3s = sqrt(df_data_code$Qm3s)
    }
    
    maxQ = max(df_data_code$Qm3s, na.rm=TRUE)
    
    power = get_power(maxQ)
    
    dbrk = 10^power

    df_data_code$Qm3sN = df_data_code$Qm3s / dbrk

    if (!is.null(df_trend_code)) {
        
        df_trend_code$trendN = df_trend_code$trend / dbrk
        df_trend_code$interceptN = df_trend_code$intercept / dbrk
    }

    maxQN = max(df_data_code$Qm3sN, na.rm=TRUE)
    
    if (maxQN >= 5) {
        dbrk = 1.0
        accuracy = 0.1
    } else if (maxQN < 5 & maxQN >= 3) {
        dbrk = 0.5
        accuracy = 0.1
    } else if (maxQN < 3 & maxQN >= 2) {
        dbrk = 0.4
        accuracy = 0.1
    } else if (maxQN < 2 & maxQN >= 1) {
        dbrk = 0.2
        accuracy = 0.1
    } else if (maxQN < 1) {
        dbrk = 0.1
        accuracy = 0.1
    }
    
    dDate = as.numeric(df_data_code$Date[length(df_data_code$Date)] -
                       df_data_code$Date[1]) / unit2day

    # datebreak = round(as.numeric(dDate) / unit2day / 11 , 0)
    
    if (dDate >= 100) {
        datebreak = 25
        dateminbreak = 5
    } else if (dDate < 100 & dDate >= 50) {
        datebreak = 10
        dateminbreak = 1
    } else if (dDate < 50) {
        datebreak = 5
        dateminbreak = 1
    }
    
    p = ggplot() +