Question: Calculating correlations using data from different experimental methods for a co-expression network
0
gravatar for pixie@bioinfo
5 months ago by
pixie@bioinfo1.4k
pixie@bioinfo1.4k wrote:

Hello, I have data from RNA-seq as raw read counts from 3 replicates from control and 3 replicates from treatment. I have to calculate Pearson correlations (PCCs) among the genes after using usual normalization procedures. One of the important transcription factor which is important for the biological hypothesis does not have any reads mapped to it, except in one of the replicates.

The biologists want me to use the expression data from RT-PCR experiment for the particular TF for calculating the PCCs, while using raw read counts for other genes. Is this okay to do or is there any solution to this ? Thanks.

pcr rna-seq • 178 views
ADD COMMENTlink modified 5 months ago • written 5 months ago by pixie@bioinfo1.4k
2

No that doesn't sound okay. What biologists want is not always possible, I have this problem as bioinformatician also all the time.

ADD REPLYlink written 5 months ago by Benn7.5k
1

Unfortunately, Its becoming very hard to explain them the logic behind it -_-

ADD REPLYlink written 5 months ago by pixie@bioinfo1.4k

I am afraid that it is part of the job.

ADD REPLYlink written 5 months ago by Benn7.5k

Sounds like you have a RNA-seq expression matrix but that it does not include the key gene of interest? So, your colleagues asked you to, instead, use RT-PCR data for this gene of interest?

ADD REPLYlink written 5 months ago by Kevin Blighe47k

Correct, zero read counts for the given gene, so plug in RT-PCR data

ADD REPLYlink written 5 months ago by pixie@bioinfo1.4k

I mean, you could just correlate this gene's RT-PCR data to all other genes. As it is correlation, it is not sooo affected by the different distributions. I would use Spearman, though, not Pearson. Your sample n is low, so, non-parametric correlation is a better fit.

I imagine that RT-PCR data is normalised and log2 transformed?

ADD REPLYlink written 5 months ago by Kevin Blighe47k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1003 users visited in the last hour