Not necessary. It is just an example. But, in general, more the replicates, better would be your statistical power in deducing biological function. One biological condition and one replicate? But text file supplied has 12 cel files. Does each cel file represent one time point? i.e 12 cel files =12 time points x 1 experiment?
[This is a follow-up to my previous post]
Hello world,
I have a set of genes whose RNA expression is catalogued in development (FlyBase temporal RNA-seq). I want to find which genes go up and down at time point A_07 and possibly cluster genes to their patterns also.
I put this into R and used the "Timecourse package" https://bioconductor.org/packages/release/bioc/vignettes/timecourse/inst/doc/timecourse.pdf
I've tried to replicate the script as close as I can including labelling the columns accordingly.
Here is the following matrix: 
I used the mb.long function and the result:

The results come out to be "NaN" with no error message by R. I've tried expanding the prior.df function from 1-5 but nothing changes this result. Whats going wrong?
1 answer
I only have 1 biological condition and only 1 replicate. Do you need multiple replicates?
Here is the data: http://m.uploadedit.com/bbtc/1512924884525.txt
I used your code and got the following. "out1" is giving me the exact same thing and when I try to plot it, nothing shows. The plot image is for "out2"

Thank you so much for your time
12 time points x 1 experiment. [I have a list of 100 genes, this is just a sample]
So I dont expect a plot from plotProfile with my results? How can I get a clustered image like this one:

I tried one replicate per time point using fruitfly data. mb.long fails to run. I tried the same using simulated data and again it fails. I am not sure if it is possible to run timecourse package function with 1 sample per timepoint and one condition .
Hmm must be so...
Do you know of another way of producing a similar plot as the one I recently posted? Im gonna start a new thread on this
Log in to answer this question.
Please post some example data (not image).
I created a dummy data set for this. Please go through the code below and try to compare with your code. Probably the data structure in R might be incorrect. Check your data types. I added comments to the code. Simulated data has 4 times points (1-4) and two conditions (Normal and tumor represented by A and B respectively). Gene names are 4 (gene1-4). Hence data frame is 4x8. Compare your data frame with final data frame in the example (name: total_exp)