Hi all, I am using wgcna for coexpression analysis of time series data. I have 48 samples in total (16 time points with 3 replicates). I have identified different modules based on gene correlation across all samples. After this I want to associate these modules to trait of interest. I have categorized these 16 time points into 5 time zones so I want to see which are the different modules which are associated with these 5 time zones. My trait data is:
samples timezone1 timezone2 timezone3 timezone4 timezone5
dark11 1 0 0 0 0
dark11_1 1 0 0 0 0
dark9 1 0 0 0 0
dark7 0 1 0 0 0
dark7_1 0 1 0 0 0
dark5 0 0 1 0 0
dark5_1 0 0 1 0 0
and so on
I want to ask if I am doing it right? or is there any other way to associate the traits to modules in categorical data? After identifying the modules related to every time zone, I want to identify intramodular connectivity and hub genes.