cycle (version 1.26.0)

ar1analysis: Performs AR1 fitting

Description

Calculation of the autocorrelation coefficients genes and variance of corresponding random variables to fit gene expression time series by AR1 processes

Usage

ar1analysis(eset)

Arguments

eset

object of the class “ExpressionSet”

Value

alpha) for ExpressionSet features and variance (sigma2) of corresponding random variables obtained using the ar function of the stats package.

See Also

ar

Examples

data(yeast) # loading the reduced CDC28 yeast set (from the Mfuzz package)

# Data preprocessing 
if (interactive()){
data(yeast) 

yeast <- filter.NA(yeast) 
# filters genes with more than 25% of the expression values missing 

yeast  <- fill.NA(yeast) 
# for illustration only; rather use knn method for replacing missing values

tmp <- ar1analysis(yeast) 
# fits AR1 process autocorrelation coefficients  

plot(density(tmp$alpha),main="Autocorrelation")

}