# SetupExample.R, efg, 18 July 2005
# Create sample input File showing how microarray data can be processed

SampleCount  <- 11
TimeCount    <- 25  # Assume 25 points over 24 hours

set.seed(19)
Time <- sort(round( runif(TimeCount, 0, TimeCount), 5))

Expression <- matrix(0, SampleCount, TimeCount)

colnames(Expression) <- paste("Hour", 0:(TimeCount-1), sep="")
rownames(Expression) <- paste("Series", 1:SampleCount, sep="")

Expression[ 1,] <- 1  # constant:  not periodic
Expression[ 2,] <- cos(2*pi*Time/8  + pi/6)
Expression[ 3,] <- cos(2*pi*Time/12 + pi/2)

Expression[ 4,] <- cos(2*pi*Time/24 + 0)
Expression[ 5,] <- cos(2*pi*Time/24 + pi/4)
Expression[ 6,] <- cos(2*pi*Time/24 + pi/2)
Expression[ 7,] <- cos(2*pi*Time/24 + 3*pi/4)
Expression[ 8,] <- cos(2*pi*Time/24 + pi)
Expression[ 9,] <- cos(2*pi*Time/24 + 3*pi/2)

Expression[10,] <- cos(2*pi*Time/8  - pi/6) + cos(2*pi*Time/40)
Expression[11,] <- cos(2*pi*Time/12 + pi/2) + cos(2*pi*Time/18 + pi)

Expression <- round(Expression, 5)
write.csv(Expression, file="LombScargleExampleSet.csv")