gistfile1.txt

mymts<-dataset%>%
  dplyr::select(VO2i_methode1,DO2i_methode1,RQ_methode1)%>%
  ts(., frequency=31557600)
plot(ts(mymts))
mod.foreca <- foreca(mymts, n.comp = 3, plot = TRUE, spectrum.control=list(method="pspectrum"))
mod.foreca
plot(mod.foreca)
biplot(mod.foreca)
summary(mod.foreca)
mod.foreca$scores <- ts(mod.foreca$scores, start = start(mymts),freq = frequency(mymts))
plot(mod.foreca$scores)
round(cor(mod.foreca$scores), 3)
spec <- mvspectrum(mod.foreca$scores, "pspectrum");plot(spec)