R/Radiation_helper_functions.R

In rhlee12/Noble: Work With NEON Project Instrumented Data

    #
    # #subset to first and last weeks
    # first.pop=var.data[frst.per[2]>var.data$startDateTime&var.data$startDateTime>=frst.per[1],]
    # last.pop=var.data[last.per[2]>var.data$startDateTime&var.data$startDateTime>=last.per[1],]
    #
    # #Subset to nightime conditions
    #
    #
    #
    #
    #
    # first.pop=first.pop[]
    # last.pop=last.pop[which(strftime(last.pop$startDateTime, format="%H:%M:%S", tz=time.zone) %in% test.time),]
    #
    # ## Remove all NAs, do this after subsetting to nightime conditions
    # first.pop=first.pop[,-which(colSums(is.na(first.pop))==length(first.pop[,1]))]
    # last.pop=last.pop[,-which(colSums(is.na(last.pop))==length(last.pop[,1]))]
    #
    # f.test=stats::var.test(unlist(as.list(first.pop[,(2:length(colnames(first.pop)))])),# ------>f.test results####
    #                        unlist(as.list(last.pop[,(2:length(colnames(last.pop)))])))
    #
    # if(0.01>f.test$p.value){f.test.result="Fail"}else{f.test.result="Pass"} ################################################################
    #
    # message(paste0("Variance Stability Test: ", f.test.result))
    # utils::write.csv(x = data.frame(value=unlist(f.test)),file = paste0(raw.dir, "variance_stats.csv"))
    #

    #If variance has failed, make a plot of all variances for the full testg period
    # if(f.test.result=="Fail"&plot){
    #
    #     night.data=var.data[which(strftime(var.data$startDateTime, format="%H:%M:%S", tz=time.zone) %in% test.time),]
    #     m.night=reshape2::melt(night.data, id.vars="startDateTime")
    #     m.night$startDateTime=as.POSIXct(m.night$startDateTime, tz = time.zone)
    #     m.night$testPeriod=0
    #     m.night$testPeriod[frst.per[2]>m.night$startDateTime&m.night$startDateTime>=frst.per[1]]=1
    #     m.night$testPeriod[last.per[2]>m.night$startDateTime&m.night$startDateTime>=last.per[1]]=1
    #
    #     var.plot=ggplot2::ggplot(data=m.night, ggplot2::aes(x=startDateTime, y=value, color=as.factor(variable)))+
    #         ggplot2::geom_rect(
    #             ggplot2::aes(xmin=as.POSIXct(frst.per[1]),
    #                 xmax = as.POSIXct(frst.per[2]),
    #                 ymin = -Inf,
    #                 ymax = Inf),
    #             fill = 'gray',
    #             alpha = 0.01)+
    #         ggplot2::geom_rect(
    #             ggplot2::aes(xmin=as.POSIXct(last.per[2]),
    #                 xmax = as.POSIXct(last.per[1]),
    #                 ymin = -Inf,
    #                 ymax = Inf),
    #             fill = 'gray',
    #             alpha = 0.01)+
    #         ggplot2::geom_point()+
    #         ggplot2::theme_light()+
    #         ggplot2::ggtitle(paste0(site, " overnight variance values for radiation data products"), subtitle = "Gray boxes indicate test periods")+
    #         ggplot2::xlab("Date")+
    #         ggplot2::geom_smooth()
    #
    #     ggplot2::ggsave(plot = var.plot, device = "png", width = 8, height = 5, filename = paste0(site, "_night_vars.png"), path = raw.dir, units = "in")
    #
    #     #lm(m.night$startDateTime ~ m.night$value)
    # }

    #return(f.test.result)