Hi, I have a FPKM matrix of 9644 columns and 15 rows from a RNA-seq expression experiment. The columns are the gene and the rows represents expression of the gene at different time points. My basic idea is to identify transcription factor binding site (TFBS) upstream of target gene. I want to do de-novo motif discovery based on over-represented sequence search in regulatory regions of target gene along with a group of co-expressed gene. Therefore, my first idea is to search genes co-expressed with my target gene. How can I do the co-expression analysis in R. How can I visualize the results? Please help me with R scripts. Let us consider the data name is 'Rdata.csv'