library(data.table)
library(tidyverse)
library(magrittr)
library(fda)
library(gridExtra)
library(emuR)
fpcaRootDir <- "C:/Users/49178/Documents/FPCA35_all/scripts/"
source(file.path(fpcaRootDir, "defint.fd.R"))
source(file.path(fpcaRootDir, "fdPar.R"))
source(file.path(fpcaRootDir, "getPCscores.R"))
source(file.path(fpcaRootDir, "landmarkreg.nocurve.R"))
library(data.table)
library(tidyverse)
library(magrittr)
library(fda)
library(gridExtra)
library(emuR)
fpcaRootDir <- "C:/Users/49178/Documents/FPCA35_all/scripts/"
#pfad1 = "//phonetik.uni-muenchen.de/vdata/ERC2/Italian/Feb20"
pfad2 = "C:/Users/PGreca/Desktop/FPCA35_all"
db = load_emuDB(file.path(pfad2, "db_emuDB"))
setwd("//phonetik.uni-muenchen.de/vdata/ERC2/Italian/CueTradingPaper"
db = load_emuDB("db_emuDB"))
setwd("//phonetik.uni-muenchen.de/vdata/ERC2/Italian/CueTradingPaper")
db = load_emuDB("db_emuDB")
setwd("//phonetik.uni-muenchen.de/vdata/ERC2/Italian/CueTradingPaper/POST-REVIEW-SUBMITS/OpenData-LMU/Secondary_dataset")
e.s = query(db, "[MAU=e|je & End(Word, MAU)==0]", resultType = "emuRsegs")%>%sort()
e.w.s = requery_hier(db, e.s, "Word", resultType = "emuRsegs")%>%sort()
# remove cervelli, cervello, tenete, prete
temp = label(e.w.s) %in% c("cervelli", "cervello", "tenete", "prete", "preti")
e.s = e.s[!temp,]
e.w.s = e.w.s[!temp,]
# get /e/ in these words
e.cer = query(db, "[[#MAU = e ^ Word = cervelli|cervello] -> MAU = ddZ]", resultType = "emuRsegs")%>%sort()
e.cer.w = requery_hier(db, e.cer, "Word", resultType = "emuRsegs")%>%sort()
e.ten = query(db, "[[#MAU = e ^ Word = tenete] -> MAU = t]", resultType = "emuRsegs")%>%sort()
e.ten.w = requery_hier(db, e.ten, "Word", resultType = "emuRsegs")%>%sort()
e.preti = query(db, "[[#MAU = e ^ Word = prete|preti] -> MAU = v]", resultType = "emuRsegs")%>%sort()
e.preti.w = requery_hier(db, e.preti, "Word", resultType = "emuRsegs")%>%sort()
e.s = rbind(e.s, e.cer, e.ten, e.preti)
e.w.s = rbind(e.w.s, e.cer.w, e.ten.w, e.preti.w)
# check same number of rows
nrow(e.s) == nrow(e.w.s)
# check start time of e.s is always bigger than start time of word
# should all be F
any(start(e.s) < start(e.w.s))
# /E/
E.s = query(db, "[MAU=E & End(Word, MAU)==0]",resultType = "emuRsegs")%>%sort()
E.w.s = requery_hier(db, E.s, "Word", resultType = "emuRsegs")%>%sort()
nrow(E.s) == nrow(E.w.s)
any(start(E.s) < start(E.w.s))
# /o/
o.s = query(db, "[MAU=o|O|jo & End(Word, MAU)==0]",resultType = "emuRsegs")%>%sort()
o.w.s = requery_hier(db, o.s, "Word", resultType = "emuRsegs")%>%sort()
nrow(o.s) == nrow(o.w.s)
any(start(o.s) < start(o.w.s))
# /a/
a.s = query(db, "[MAU=a & End(Word, MAU)==0]", resultType = "emuRsegs")%>%sort()
