09.11(Fri)

[R code]

library(lubridate)

data<-read.csv("/Users/amy/Desktop/BIgCon/data processing/modeldata_nooutlier2.csv") data<-data[,c("date","ratings")] wday <- factor(wday(date),labels=c("Sun","Mon","Tue","Wed","Thu","Fri","Sat")) month <- month(date)

X = list() date <- unique(data$date) for ( i in 1:365){ x1 <- data %>% filter(date==date[i]) x<-as.vector(x1$ratings) X = append(X,list(x)) } cluster = tsclust(X, k=5, distance="dtw_basic", type="hierarchical") cl = slot(cluster, "cluster") plot(cluster, type="sc") table(month,cluster@cluster)