Skip to content
Extraits de code Groupes Projets
Valider cb4a194c rédigé par Michel Crucifix's avatar Michel Crucifix
Parcourir les fichiers

plot routines

parent c2e1e7a6
Aucune branche associée trouvée
Aucune étiquette associée trouvée
Aucune requête de fusion associée trouvée
...@@ -8,7 +8,7 @@ Maintainer: Michel Crucifix <michel.crucifix@uclouvain.be> ...@@ -8,7 +8,7 @@ Maintainer: Michel Crucifix <michel.crucifix@uclouvain.be>
Description: Collection of time series analysis tools, focused on spectral decomposition and reconstruction, with an original focus on astronomical theory of paleoclimates and cyclostratigraphy. Description: Collection of time series analysis tools, focused on spectral decomposition and reconstruction, with an original focus on astronomical theory of paleoclimates and cyclostratigraphy.
License: MIT licence License: MIT licence
LazyLoad: yes LazyLoad: yes
RoxygenNote: 7.2.3 RoxygenNote: 7.3.1
Imports: Rdpack, sfsmisc, cmna, gsignal, pracma Imports: Rdpack, sfsmisc, cmna, gsignal, pracma
RdMacros: Rdpack RdMacros: Rdpack
Depends: Depends:
......
...@@ -10,11 +10,13 @@ export(cwt_morlet) ...@@ -10,11 +10,13 @@ export(cwt_morlet)
export(hilbert_extension) export(hilbert_extension)
export(mem) export(mem)
export(mfft) export(mfft)
export(mfft_anova)
export(mfft_real) export(mfft_real)
export(mfft_real_C) export(mfft_real_C)
export(mfft_real_ter) export(mfft_real_ter)
export(powerspectrum.wavelet) export(powerspectrum.wavelet)
export(reconstruct_mfft) export(reconstruct_mfft)
export(reconstruct_morlet)
export(ssa) export(ssa)
importFrom(Rdpack,reprompt) importFrom(Rdpack,reprompt)
importFrom(cmna,goldsectmax) importFrom(cmna,goldsectmax)
......
...@@ -93,6 +93,9 @@ enhance_ridge <- function(A,R,confidence=0.95,inter=0.08) ...@@ -93,6 +93,9 @@ enhance_ridge <- function(A,R,confidence=0.95,inter=0.08)
} }
#' @export reconstruct_morlet
reconstruct_morlet <- function(W,scales=c(-Inf,Inf), periods=NULL) reconstruct_morlet <- function(W,scales=c(-Inf,Inf), periods=NULL)
{ {
if (!(attr(W,"class")=="wavelet")) stop ("object is not a wavelet transform") if (!(attr(W,"class")=="wavelet")) stop ("object is not a wavelet transform")
......
...@@ -88,7 +88,7 @@ mfft <- function(xdata, minfreq=NULL, maxfreq=NULL, flag=1, nfreq=30) { ...@@ -88,7 +88,7 @@ mfft <- function(xdata, minfreq=NULL, maxfreq=NULL, flag=1, nfreq=30) {
# take care and verify that this actually works this way # take care and verify that this actually works this way
OUT = data.frame(Freq=Freq, Ampl=Ampl, Phase=Phase) OUT = data.frame(Freq=Freq, Ampl=Ampl, Phase=Phase)
attr(OUT,"class") = "mfft_deco" attr(OUT,"class") = "mfft_deco"
return(data.frame(Freq=Freq, Ampl=Ampl, Phase=Phase)) return(OUT)
} }
# #
......
...@@ -182,7 +182,7 @@ reconstruct_mfft <- function(M){ ...@@ -182,7 +182,7 @@ reconstruct_mfft <- function(M){
#' MFFT ANOVA #' MFFT ANOVA
#' @rdname mfft_deco #' @rdname mfft_deco
#' @export mfft_anova #' @export mfft_anova
reconstruct_mfft <- function(M){ mfft_anova <- function(M){
if (!(attr(M,"class")=="mfft_deco")) stop ("object is not a MFFT decomposition") if (!(attr(M,"class")=="mfft_deco")) stop ("object is not a MFFT decomposition")
xdata <- attr(M,"xdata") xdata <- attr(M,"xdata")
nfreq <- attr(M,"nfreq") nfreq <- attr(M,"nfreq")
......
...@@ -2,13 +2,18 @@ ...@@ -2,13 +2,18 @@
% Please edit documentation in R/mfft_real.R % Please edit documentation in R/mfft_real.R
\name{reconstruct_mfft} \name{reconstruct_mfft}
\alias{reconstruct_mfft} \alias{reconstruct_mfft}
\alias{mfft_anova}
\alias{plot.mfft_deco} \alias{plot.mfft_deco}
\title{MFFT reconstruction} \title{MFFT reconstruction}
\usage{ \usage{
reconstruct_mfft(M, times) reconstruct_mfft(M)
mfft_anova(M)
\method{plot}{mfft_deco}(M, ...) \method{plot}{mfft_deco}(M, ...)
} }
\description{ \description{
MFFT reconstruction MFFT reconstruction
MFFT ANOVA
} }
0% Chargement en cours ou .
You are about to add 0 people to the discussion. Proceed with caution.
Terminez d'abord l'édition de ce message.
Veuillez vous inscrire ou vous pour commenter