```march19 <- as.matrix(read.csv(“C:/Users/clair/OneDrive/Documents/AAA Spring 2025/Data/Final Project/031925.csv”)) graphics::image(march19, col=gray.colors(5000, start=0, end=1)) graphics::image(march19, col = hcl.colors(1000, “YlOrRd”, rev = TRUE)) graphics::image(march19, col = hcl.colors(1000, “YlOrRd”, rev = FALSE))

graphics::hist(march19, breaks=“FD”, freq=F, main=““, xlab=”Temperature (F)“, ylab=”Frequency”, col=“gray90”)

blkwhite <- as.matrix(read.csv(“C:/Users/clair/OneDrive/Documents/AAA Spring 2025/Data/Final Project/blkwhitecontrol.csv”)) brds_0402 <- as.matrix(read.csv(“C:/Users/clair/OneDrive/Documents/AAA Spring 2025/Data/Final Project/040225.csv”))

graphics::hist(blkwhite, breaks=“FD”, freq=F, main=““, xlab=”Temperature (F)“, ylab=”Frequency”, col=“gray90”) clsize <- 4 cluster_colors <- hcl.colors(clsize, “YlOrRd”) cluster_mat <- matrix(march19, nrow=dim()[1], ncol = dim()[2]) for (i in 1:clsize) { cluster_ordered[tt$cluster==center_order[i]] <- clsize+1-i }

graphics::hist(brds_0402, breaks=“FD”, freq=F, main=““, xlab=”Temperature (F)“, ylab=”Frequency”, col=“gray90”)

graphics::image(march19, col = hcl.colors(1000, “YlOrRd”, rev = TRUE)) ```