###################################################
### chunk number 1: 
###################################################
  library(cycle) 


###################################################
### chunk number 2: 
###################################################
set.seed(1)


###################################################
### chunk number 3: 
###################################################
data(yeast)
  yeast <- yeast[1:200,]


###################################################
### chunk number 4: 
###################################################
yeast <- filter.NA(yeast, thres=0.25)


###################################################
### chunk number 5: 
###################################################
yeast <- fill.NA(yeast,mode="mean")


###################################################
### chunk number 6: 
###################################################
yeast <- standardise(yeast)


###################################################
### chunk number 7: 
###################################################
auto.corr <- 0
for (i in 2:dim(exprs(yeast))[[2]]){
 auto.corr[i] <- cor(exprs(yeast)[,i-1],exprs(yeast)[,i])
 }
auto.corr 


###################################################
### chunk number 8: 
###################################################
T.yeast <- 85  


###################################################
### chunk number 9: 
###################################################
times.yeast <-  pData(yeast)$time 
times.yeast


###################################################
### chunk number 10: 
###################################################
NN <- 100 


###################################################
### chunk number 11: 
###################################################
fdr.rr <- fdrfourier(eset=yeast,T=T.yeast,times=times.yeast,background.model="rr",N=NN,progress=FALSE)


###################################################
### chunk number 12: 
###################################################
fdr.ar1 <- fdrfourier(eset=yeast,T=T.yeast,times=times.yeast,background.model="ar1",N=NN,progress=FALSE)


###################################################
### chunk number 13: 
###################################################
sum(fdr.rr$fdr < 0.25) 
sum(fdr.ar1$fdr < 0.25)


###################################################
### chunk number 14: 
###################################################
fdr.ar1$fdr[which(fdr.ar1$fdr < 0.25)]