install.packages("gridExtra")
install.packages("devtools")
devtools::install_github("IPS-LMU/emuR")

library(emuR)
library(wrassp)
library(tidyverse)

path2ae = file.path(mypath, "emuR_demoData", "ae_emuDB")
# load database
ae = load_emuDB(path2ae, verbose = F)
list_ssffTrackDefinitions(ae) #has dft

#add_ssffTrackDefinition(ae,"dft",onTheFlyFunctionName = dftSpectrum)

sS.sl = query(ae, "[Phonetic==s|S]")

sS.dft = get_trackdata(ae, sS.sl, "dft", resultType = "tibble",cut=0.5) 

View(sS.dft)

sS.dft.long = convert_wideToLong(sS.dft,calcFreqs = TRUE)

View(sS.dft.long)

ggplot(sS.dft.long)+
  aes(x=freq,y=track_value)+
  geom_line() +
  facet_wrap(~sl_rowIdx+labels)

sS.dft.long.mean = sS.dft.long %>%
  group_by(labels,freq) %>%
  summarise(track_value=mean(track_value))

ggplot(sS.dft.long.mean)+
  aes(x=freq,y=track_value,col=labels)+
  geom_line()

sS.dft.long.mean[sS.dft.long.mean$freq>2000 & sS.dft.long.mean$freq<3000,]

sS.dft.2000to3000 = sS.dft.long %>%
  filter(freq>2000 & freq <3000) %>%
  group_by(labels,sl_rowIdx) %>%
  summarise(track_value=mean(track_value))

ggplot(sS.dft.2000to3000) +
  aes(x=labels,y=track_value)+
  geom_boxplot()


sS.dft.long.mean = sS.dft.long.mean %>%
  group_by(labels) %>%
  mutate(reconstructed = emuR::dct(track_value,fit=T))

sS.dft.long.mean = sS.dft.long.mean%>%
  group_by(labels)%>%
  mutate(coefficients = dct(track_value),
         reconstructed = dct(track_value,fit=T),
         smoothed = dct(track_value,m=5,fit=T))

ggplot(sS.dft.long.mean)+
  geom_line(aes(x=freq,y=track_value,col=labels)) +
  geom_line(aes(x=freq,y=smoothed,col=labels))


ggplot(sS.dft.long.mean)+
  aes(x=freq,y=track_value,col=labels)+
  geom_line() +
  geom_line(aes(y=smoothed),lwd=2)


ggplot(sS.dft.long) +
  aes(x = freq, y = track_value,col=labels) +
  geom_line() +
  facet_wrap( ~ sl_rowIdx + labels)

sS.dft.long = sS.dft.long %>%
  group_by(sl_rowIdx) %>%
  mutate(smoothed = emuR::dct(track_value,m=4,fit=T))

ggplot(sS.dft.long) +
  aes(x = freq, y = track_value,col=labels) +
  geom_line() +
  geom_line(aes(y=smoothed),lwd=2) +
  facet_wrap( ~ sl_rowIdx + labels)

sS.dctCoefficients =
  sS.dft.long %>%
  group_by(sl_rowIdx,labels)%>%
  do(data_frame(DCT=dct(.$track_value,m=3,fit=FALSE)))%>%
  mutate(Kodierung=c("k0","k1","k2","k3"))%>%
  spread(Kodierung,DCT)