I'm working with a circadian timeseries dataset generated from repeated blood draws from a group of subjects. So for each subject, I have gene expression data from multiple timepoints throughout the day. Looking at the documentation for the RAIN package (https://www.bioconductor.org/packages/release/bioc/html/rain.html), it looks like I can perform a repeated measures analysis for the data from each subject by setting the 'longitudinal' argument to TRUE.

What I want to know is if this option (or RAIN in general) supports pooling data from multiple subjects (i.e. biological replicates), while accounting for the repeated measures within each subject?

In linear modeling, this would be similar to using a mixed effects model that incorporates all of the data, with time as a fixed effect and subject as a random effect.