Question: Selecting the highest Pearson correlation coefficient of two matrices
1
gravatar for jivarajivaraj
15 months ago by
jivarajivaraj50 wrote:

Hi,

I have two matrices of normalised gene expression values from single cell data from two time points. In rows I have the same genes but cells in each matrix is not same(different time point). In figure 2 C of this paper https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5668937/ authors say they have calculated Pearson correlation and give one Pearson correlation coefficient for each comparison as the highest correlation coefficient.

I am using many Pearson calculation like corrplot.mixed(cor(cbind(h_160,h_140)), order="hclust", tl.col="black") that gives me one Pearson correlation coefficient for each pair of cells (A matrix of Pearson correlation coefficient). So how I could select the highest Pearson correlation coefficient for these two matrices??

ADD COMMENTlink modified 14 months ago by Kevin Blighe51k • written 15 months ago by jivarajivaraj50
2
gravatar for Kevin Blighe
14 months ago by
Kevin Blighe51k
Kevin Blighe51k wrote:

The maximum would just be:

cormatrix <- cor(cbind(h_160,h_140))
max(cormatrix, na.rm=TRUE)

Note that I have 'parallellised' the cor() function in R, here (see section 2 cor function, parallellised): R functions edited for parallel processing

In addition, it looks like they are just doing a simple network plot in Figure 2C. I also use Pearson or Spearmann correlation coefficients in my tutorial, here: Network plot from expression data in R using igraph

Kevin

ADD COMMENTlink written 14 months ago by Kevin Blighe51k
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: 1989 users visited in the last hour