[R] correlation between rows of different data frames

Bert Gunter gunter.berton at gene.com
Wed Aug 28 16:37:03 CEST 2013


You appear to have no idea what you are doing. I suggest you post on
the Bioconductor list, not here, for guidance. Better yet, get help
locally from someone who does know what they are doing.

Cheers,
Bert

On Wed, Aug 28, 2013 at 6:36 AM, Robin Mjelle <robinmjelle at gmail.com> wrote:
> Hi,
>
> I have two data frames with time serie datamatrix. I want to pick a row X
> from the first matrix and see if it correlates with row Y in the second
> matrix. These are gene expression values and I probaly need to do some
> scaling first, but I wonder if you have any suggestions on how to do the
> correlation. A graph for each correlation would also be nice.
>
>> head(dataH)
>              T0h T0.25h T0.5h T1h T2h T3h T6h T12h T24h T48h C0h C0.25h
> C0.5h C1h C2h C3h C6h
> NM_001001130  68     95    56  43  66  62  68   90   63   89  65     85
> 58  49  81  49  76
> NM_001001144   0      1     4   0   1   1   1    4    1    2   1      3
> 1   0   6   3   0
> NM_001001152  79    129    52  50  24  45 130  154  112  147  56     85
> 67  33  52  31  93
> NM_001001160   1      1     2   0   0   1   0    0    0    1   0      1
> 2   4   2   1   0
> NM_001001176   0      0     0   0   0   0   0    0    0    0   0      0
> 0   0   0   0   0
> NM_001001177   1      3     2   3   0   1   1    0    2    0   1      6
> 4   1   3   0   0
>              C12h C24h C48h
> NM_001001130   73   48   77
> NM_001001144    1    2    1
> NM_001001152   77   65  109
> NM_001001160    2    1    0
> NM_001001176    0    0    0
> NM_001001177    3    0    2
>> head(dataM)
>                   T0h T0.25h T0.5h   T1h   T2h   T3h   T6h  T12h  T24h
> T48h   C0h C0.25h
> mmu-let-7a-1-3p   450    329   395   300   469   320   245   281   300
> 376   166    242
> mmu-let-7a-2-3p     0      0     0     0     0     0     0     0     0
> 0     0      0
> mmu-let-7a-5p   93070  62871 65013 69714 73983 60836 49885 63788 62161
> 65595 48525  32111
> mmu-let-7b-3p      54     14    35    15    47     7     9    13    26
> 50    17     15
> mmu-let-7b-5p      19     16    18    25    18    16    19    16    18
> 42    14      8
> mmu-let-7c-1-3p     0      0     0     0     0     0     0     0     0
> 0     0      0
>                 C0.5h   C1h   C2h   C3h   C6h  C12h  C24h  C48h
> mmu-let-7a-1-3p   211   141   127    98   113   170   187   249
> mmu-let-7a-2-3p     0     0     0     0     0     0     0     0
> mmu-let-7a-5p   30148 17108 22988 27513 37940 57235 40395 38018
> mmu-let-7b-3p      15     3     7     9    12    23    20    14
> mmu-let-7b-5p       8     2    12    10     7    15     6     6
> mmu-let-7c-1-3p     0     0     0     0     0     0     0     0
>
>         [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.



-- 

Bert Gunter
Genentech Nonclinical Biostatistics

Internal Contact Info:
Phone: 467-7374
Website:
http://pharmadevelopment.roche.com/index/pdb/pdb-functional-groups/pdb-biostatistics/pdb-ncb-home.htm



More information about the R-help mailing list