ONR-Spectral-Sensitivity-Function-/kssAnalysis.R at master · LightingResearchCenter/ONR-Spectral-Sensitivity-Function- · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
library(Rmisc)
library(readr)
library(ggplot2)
library(readxl)
library(reshape2)
source("C:/Users/roohac/Documents/GitHub/ONR-Spectral-Sensitivity-Function-/kssAnalysisFunctions.R")


nighttime_KSS_normalized_time2 <- read_excel("//root/projects/ONR-EEG-BAA16_001/PROCESSED_DATA/KSS/nighttime_KSS_normalized_time2.xlsx")


nighttime_KSS_normalized_time2_long <- melt(nighttime_KSS_normalized_time2, id.vars = c("Subject_id", "session", "date", "color","condition","group","month"))


nighttime_KSS_normalized_time2_long2 <- subset(nighttime_KSS_normalized_time2_long, Subject_id != 122)
#1
kss <- nighttime_KSS_normalized_time2_long2
analyzeKSS_normDim(kss, TRUE)

#2
kss <- nighttime_KSS_normalized_time2_long2
analyzeKSS(kss, TRUE)

#3
red <- subset(nighttime_KSS_normalized_time2_long2, color == "red")
analyzeKSS(red, TRUE)

#4
blue <- subset(nighttime_KSS_normalized_time2_long2, color == "blue")
analyzeKSS(blue, FALSE)

#5
green <- subset(nighttime_KSS_normalized_time2_long2, color == "green")
analyzeKSS(green, FALSE)

#6
amber<- subset(nighttime_KSS_normalized_time2_long2, color == "amber")
analyzeKSS(amber, FALSE)

#7
cyan <- subset(nighttime_KSS_normalized_time2_long2, color == "cyan")
analyzeKSS(cyan, FALSE)