GeneCycle (version 1.1.6)

dominant.freqs: Dominant Frequencies in Multiple (Genetic) Time Series

Description

dominant.freqs returns the m dominant frequencies (highest peaks) in each of the periodogram computed for the individual time series.

Usage

dominant.freqs(x, m=1, ...)

Value

A matrix (or vector, if only 1 time series is considered) with the dominant frequencies. In a matrix, each column corresponds to one time series.

Arguments

x

multivariate (genetic) time series (each column of this matrix corresponds to a separate variable/time series), or a vector with a single time series

m

number of dominant frequences

...

arguments passed to periodogram

Author

Konstantinos Fokianos and Korbinian Strimmer (https://strimmerlab.github.io/).

See Also

periodogram, spectrum.

Examples

# load GeneCycle library
library("GeneCycle")

# load data set
data(caulobacter)

# how many samples and how many genes?
dim(caulobacter)

# first three dominant frequencies for each gene
dominant.freqs(caulobacter, 3)

# first four dominant frequencies for gene no. 1000
dominant.freqs(caulobacter[,1000], 4)