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

In GeneCycle: Identification of Periodically Expressed Genes

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, ...)

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

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.

Author(s)

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)

GeneCycle documentation built on April 3, 2025, 10:53 p.m.