matrix - Apply a function to all column pairs of two matrices in R -


i have 2 matrices same columns , rows names:

> metilacion[1:5,1:5]             a2bp1       a2m     a2ml1     a4galt       aaas paciente1   0.2804884 0.5816559 1.1814702 -0.6234276 -0.3997400 paciente2   0.5122471 1.2944264 0.5673766  0.4490407 -0.6045510 paciente3  -0.3116356 1.6085049 0.9970350  0.3379215 -0.4787046 paciente4  -0.7220941 0.8771948 2.1445474 -0.5837802 -0.4848246 paciente5  -0.3369999 1.5885716 0.8185654  0.2374583 -0.5698858 > expresion[1:5,1:5]            a2bp1         a2m     a2ml1      a4galt       aaas paciente1 -0.9082274 -0.17736185 0.8846485 -0.36059775 -0.5624139 paciente2 -1.7152290  1.62368019 0.3292617  1.35968899 -0.9220157 paciente3 -1.0581859  0.33028098 1.1020073  0.01870851 -0.9669236 paciente4 -0.8389615  1.33754885 0.5122861 -0.14583960 -0.8196533 paciente5 -1.5273835  0.06418637 0.2695209  0.03381359 -0.4461490

i want compute correlation coefficient between pair of columns between 2 matrices , generate object correlation values every column pairs.

for example, correlation coefficient between first column be:

> cor(metilacion[,1],expresion[,1]) [1] -0.09351992

so, want generate object contains correlation values.

thanks!

you can cor

 cor(metilacion,expression)  #          a2bp1         a2m      a2ml1     a4galt        aaas  #a2bp1  -0.4887051  0.03682951 -0.0404260  0.5795882 -0.03534625  #a2m    -0.5909642  0.01572799 -0.1469085  0.3503903 -0.19412101  #a2ml1   0.8006633  0.17242226  0.1294179 -0.5827062 -0.05502329  #a4galt -0.8036390  0.18066923 -0.2026173  0.6824085 -0.32097886  #aaas    0.9033514 -0.54378874  0.7694163 -0.7995712  0.13676285

if need cor of corresponding columns

 diag(cor(metilacion,expresion))

or

 mapply(cor, as.data.frame(metilacion), as.data.frame(expresion))  #    a2bp1         a2m       a2ml1      a4galt        aaas   #-0.48870510  0.01572799  0.12941787  0.68240850  0.13676285

Comments

Post a Comment

Popular posts from this blog

powershell Start-Process exit code -1073741502 when used with Credential from a windows service environment -

i'm running strange behavior powershell start-process call. here call: $process = start-process ` "c:\somepath\mybinary.exe" ` -passthru ` -credential $defaultcredential ` -wait ` -workingdirectory "c:\somepath" ` -loaduserprofile if ($process.exitcode -ne 0) { #do } this call return exit code - 1073741502 . after quick search, exit code seems related generic error when program not load required dll (aka. status_dll_init_failed ). when run without -credential $credential program runs correctly. in order isolate problem, manually launched some.exe in prompt target credential , runs smoothly. so problem seems come way start-process cmdlet launch process. i found potential solutions problem tried apply no luck : link , link . would have idea of what's going on here ? edit 1: run proc mon monitoring program activities when launched directly or via powershell script. problem seems occur when loading kernelbase
Read more

twig - Using Twigbridge in a Laravel 5.1 Package -

i'm trying create laravel 5.1 package uses twigbridge instead of .php views. can normal .php views display fine, when try make use of .twig templates i'm met error. possible use twigbridge templates when creating views laravel package or pretty stuck blade , php ? edit 1: loading twig template without extending works. it's when start using {% extends 'file' %} error comes out. i error error loading /virtual/laravel/packages/username/mypackage/src/views/sample.twig: template "templates.master" not defined () in "/virtual/laravel/packages/username/mypackage/src/views/sample.twig" @ line 1. my folders username |_ mypackage |_ src |_ controllers - samplecontroller.php |_ views |_ templates - master.twig - sample.twig twig template sample.twig {% extends 'templates.master' %} {% block content %} {{ parent()
Read more

c# - LINQ join Entities from HashSet's, Join vs Dictionary vs HashSet performance -

i have hashset's each stores t, have written test application compares different relation algorithms can think of, i'm not pleased results getting. does there exist more efficient ways of achieving object relations, once testing? using system; using system.collections.generic; using system.linq; using system.text; using system.threading.tasks; using system.diagnostics; namespace linqtests { class program { static void main(string[] args) { hashset<user> usertable = new hashset<user>(); hashset<userproperty> userpropertytable = new hashset<userproperty>(); #region lets create dummy data. console.writeline("please wait while creating dummy data, can take while..."); console.writeline(""); int rows = 1000000; for(int x = 0; x < rows; x++) { random rnd = new random(); // ad
Read more