Question: Selecting the highest Pearson correlation coefficient of two matrices
1
gravatar for jivarajivaraj
8 months ago by
jivarajivaraj40 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 7 months ago by Kevin Blighe41k • written 8 months ago by jivarajivaraj40
2
gravatar for Kevin Blighe
7 months ago by
Kevin Blighe41k
Guy's Hospital, London
Kevin Blighe41k 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 7 months ago by Kevin Blighe41k
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: 1461 users visited in the last hour