Hi, everyone. I am analyzing the correlation of WGBS dataset and RNAseq dataset. Generally, more than one methylated CpGs could be detected in the transcriptional unit of a gene, I wanted to correlate the methylation level of every CpG site with the TPM of that gene. If the number of genes are not much, I could do this in excel or R, however, I have got almost 20 thousand genes, it is almost impossible to do the correlation analysis in excel. I was wondering how to achieve this in R? Any tips or help to solve my problems would be greatly appreciated. Thanks.