geneARMA-package {geneARMA} | R Documentation |
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.
Package: | geneARMA |
Type: | Package |
Version: | 1.0 |
Date: | 2009-10-07 |
License: | GPL-3 |
LazyLoad: | yes |
Timothy McMurry and Arthur Berg
Maintainer: Timothy McMurry <tmcmurry@depaul.edu>
Ning Li, et al. Functional clustering of periodic transcriptional profiles through ARMA(p,q)
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)