geneARMA-package {geneARMA}R Documentation

Simulate, model, and display data from a time-course microarray experiment with periodic gene expression

Description

Fit models for periodic gene expression to time-course microarray data in a normal mixture model framework with mean approximated by a truncated Fourier series and covariance structure modeled by an ARMA(p,q) process. Estimation is performed with the EM algorithm.

Details

Package: geneARMA
Type: Package
Version: 1.0
Date: 2009-10-07
License: GPL-3
LazyLoad: yes

Author(s)

Timothy McMurry and Arthur Berg

Maintainer: Timothy McMurry <tmcmurry@depaul.edu>

References

Ning Li, et al. Functional clustering of periodic transcriptional profiles through ARMA(p,q)

Examples

set.seed(100)
Data <- geneARMAsim(400, ars=c(.5, .1))
f1 <- geneARMAfit(Data$Y, Data$tm, 2, 2, 2, 0, eps.conv = .001, max.iter = 15, tau.init=c(.25, .45))
plot(f1, y=NULL, "all.means")
plot(f1, y=NULL, "single.cluster", j=2)

[Package geneARMA version 1.0 Index]