dominant.freqs {GeneTS}R Documentation

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 (=matrix where each column corresponds to one 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 (http://www.ucy.ac.cy/~fokianos/) and Korbinian Strimmer (http://www.stat.uni-muenchen.de/~strimmer/).

See Also

periodogram, spectrum.

Examples

# load GeneTS library
library(GeneTS)

# 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)

[Package GeneTS version 2.3 Index]