Hi,
I would like to know if there is any Bioconductor R package that can deal with RNA-seq time-series allowing for time offset between individuals. That is, I have 10 individuals, each with variable sampling (variable in both sampling time and number of samplings) throughout their pregnancy and I would like to have those genes that show a significant difference in expression throughout pregnancy. A mock design example is below:
INDIVIDUAL Sampling_day
Individual1 41
Individual1 80
Individual1 120
Individual1 150
Individual1 200
Individual2 30
Individual2 70
Individual2 180
Thanks,
João
Thanks Aaron. I am indeed expecting non-linear trends but I was not able to find how to create a spline based-matrix.
See Section 9.6.2 of the limma user's guide for an example.