Transition Analysis

library(solarf)
library(knitr) #For pretty tables on the report using kable 
library(plyr)
library(ggplot2)
library(GGally)
library(caret)
library(extrafont)

#Theme for Plot
science_theme = theme(panel.background=element_blank(), panel.grid.major = element_line(size = 0.5, color = "grey"), 
    axis.line = element_line(size = 0.7, color = "black"), legend.position = c(0.9, 
        0.7), text = element_text(size = 14))

I(Solar Centroid ID, Solar Transition)

solar.results <- readRDS(file="/Users/carlos/Google Drive/Research/SolarForecasting/datasets/mesowest/carlos/solarf/data/predictions.rds")
solar <- solar.results[,c("groundtruth.YEAR","groundtruth.MON","groundtruth.DAY","groundtruth.clusterid","istransition")]
solar$groundtruth.clusterid[is.na(solar$groundtruth.clusterid)] <- -9999
solar$istransition[solar$istransition == TRUE] <- 1
solar$istransition[solar$istransition == FALSE] <- 0
solar$istransition[is.na(solar$istransition)] <- -9999

write.table(solar[,c("groundtruth.clusterid","istransition")],file="~/Desktop/OCI/data/solarCid_solarTransition.csv",col.names=FALSE,row.names=FALSE,sep=",")

entropy <- read.table("~/Desktop/OCI/out/entropy-solarCid_solarTransition.txt",header=TRUE)
ggplot(entropy, aes(x=Lag, y=Ixy)) + geom_line(colour="#94003C") + science_theme + ggtitle("Point-to-Point Mutual Information") + xlab(expression(tau)) + ylab(expression(I(x[t],x[t+tau]))) +
  theme(
    axis.line=element_blank(),
  text = element_text(size=20, family="Interstate"),
  axis.text.x = element_text(hjust = 1, size=20,color="black"), 
  axis.text.y = element_text(hjust = 1, size=20,color="black"), 
  axis.title.y=element_text(size=20), 
  axis.title.x=element_text(size=20),
  axis.ticks = element_line(colour=NA),
  panel.grid.major = element_line(colour = 'gray', 
                                  linetype = 'dashed'),
  panel.grid.minor = element_line(colour = NA),  
  panel.background = element_rect(colour = 'white'))

I(Humidity Transition, Solar Transition)

humidity.results <- readRDS(file="/Users/carlos/Google Drive/Research/SolarForecasting/datasets/mesowest/carlos/solarf/data/predictions_Humidity.rds")
humidity <- humidity.results[,c("groundtruth.YEAR","groundtruth.MON","groundtruth.DAY","istransition")]
humidity$istransition[humidity$istransition == TRUE] <- 1
humidity$istransition[humidity$istransition == FALSE] <- 0
humidity$istransition[is.na(humidity$istransition)] <- -9999

humidity$sistransition <- solar$istransition[(solar$MON == humidity$MON) & (solar$DAY == humidity$DAY)]

pair <- merge(humidity,solar,by=c("groundtruth.YEAR","groundtruth.MON","groundtruth.DAY"))[,c("istransition.x","istransition.y")]
colnames(pair) <- c("istransition","sistransition")


write.table(pair[,c("istransition","sistransition")],file="~/Desktop/OCI/data/humidityTransition_SolarTransition.csv",col.names=FALSE,row.names=FALSE,sep=",")

entropy <- read.table("~/Desktop/OCI/out/entropy-humidityTransition_SolarTransition.txt",header=TRUE)
ggplot(entropy, aes(x=Lag, y=Ixy)) + geom_line(colour="#94003C") + science_theme + ggtitle("Point-to-Point Mutual Information") + xlab(expression(tau)) + ylab(expression(I(x[t],x[t+tau]))) +
  theme(
    axis.line=element_blank(),
  text = element_text(size=20, family="Interstate"),
  axis.text.x = element_text(hjust = 1, size=20,color="black"), 
  axis.text.y = element_text(hjust = 1, size=20,color="black"), 
  axis.title.y=element_text(size=20), 
  axis.title.x=element_text(size=20),
  axis.ticks = element_line(colour=NA),
  panel.grid.major = element_line(colour = 'gray', 
                                  linetype = 'dashed'),
  panel.grid.minor = element_line(colour = NA),  
  panel.background = element_rect(colour = 'white'))

I(Gust Transition, Solar Transition)

gust.results <- readRDS(file="/Users/carlos/Google Drive/Research/SolarForecasting/datasets/mesowest/carlos/solarf/data/predictions_Gust.rds")
gust <- gust.results[,c("groundtruth.YEAR","groundtruth.MON","groundtruth.DAY","istransition")]
gust$istransition[gust$istransition == TRUE] <- 1
gust$istransition[gust$istransition == FALSE] <- 0
gust$istransition[is.na(gust$istransition)] <- -9999

pair <- merge(gust,solar,by=c("groundtruth.YEAR","groundtruth.MON","groundtruth.DAY"))[,c("istransition.x","istransition.y")]
colnames(pair) <- c("istransition","sistransition")


write.table(pair[,c("istransition","sistransition")],file="~/Desktop/OCI/data/gustTransition_SolarTransition.csv",col.names=FALSE,row.names=FALSE,sep=",")

entropy <- read.table("~/Desktop/OCI/out/entropy-gustTransition_SolarTransition.txt",header=TRUE)
ggplot(entropy, aes(x=Lag, y=Ixy)) + geom_line(colour="#94003C") + science_theme + ggtitle("Point-to-Point Mutual Information") + xlab(expression(tau)) + ylab(expression(I(x[t],x[t+tau]))) +
  theme(
    axis.line=element_blank(),
  text = element_text(size=20, family="Interstate"),
  axis.text.x = element_text(hjust = 1, size=20,color="black"), 
  axis.text.y = element_text(hjust = 1, size=20,color="black"), 
  axis.title.y=element_text(size=20), 
  axis.title.x=element_text(size=20),
  axis.ticks = element_line(colour=NA),
  panel.grid.major = element_line(colour = 'gray', 
                                  linetype = 'dashed'),
  panel.grid.minor = element_line(colour = NA),  
  panel.background = element_rect(colour = 'white'))