Reshape data in R -


I have to be a matrix from this data frame. On the basis of values, if there is a connection between a pair of 1 gene, and if not then 0. Therefore, the value of ADRA1D and ADK will be 1, and such other pairs will also be there. But there is no pair of ADK and AR, so in that matrix this should be 0. 

  tab & lt; - read.table (text = "ID gene1 gene2 1 ADRA1D ADK 2 ADRA1B ADK 3 ADRA1A ADK 4 ADRB 1 ASIC 1 5 ADRB 1 ADK 6 ADRB 2 ASIC 1 ADRB 2 ADK 8 AGTR 1 ACH 9 AGTR 1 ADK 10 ELLOS 5 ADRB 1 11 ALOX 5 ADRB 2 12 LPPL 2 ADRB 1 13 ALPPL 2 ADRB 2 14 A2 A2 ART 1 15 AR ADARA 1 AR ADRAID 17 AR ADARA 1 B 18 AR ADARA 1A AR ADR 2 A 20 Mainly, I Want to build phylogenetic tree, so a matrix like this I do not have any price column, so how can I revive the library for this?  
 The code below does not work:  
 < Code> library (reshape) ct = cast (tab, gene1 ~ gene2)

if Using the reshape is not mandatory, so I had a look at the igraph. using the igraph package to get the symmetric matrix There is a way we First, change your data frame (the relevant 2 column) into a igraph object, and then get_adjacency is required.

Library (igraph) g < - graph.data.frame (tab [, c (2,3)] get.adjacency (g)

This gives you the proximity matrix for you to analyze the rest of the igraph Should be used. [16 x 16 spectrum matrix of class "dgCMatrix" [[16 column names 'ADRA1D', 'ADRA1B', 'suppressing ADRA1A]' ...]] ADRA1D. . . . . . . . . 1. . . . ADRA1B . . . . . . . . 1. . . . ADRA1A . . . . . . . . 1. . . . ADRB1. . . . . . . . . 1 1. . . . ADRB2. . . . . . . . . 1 1. . . . AGTR1. . . . . . . . . 1 1. . AOLX5 . 1 1. . . . . . . . . . . ALPPL2. . 1 1. . . . . . . . . . . AMY2A . . . 1. . . . . . . . . AR1 1 1. . . . . . . . . 1 1 1 ADK . . . . . . . . . . . . . . ASIC 1. . . . . . . . . . . . . . . HE . . . . . . . . . . . . . . ADORA1. . . . . . . . . . . . . . . ADRA2A . . . . . . . . . . . . . . ADRA2B . . . . . . . . . . . . . . One advantage of using

igraph is that many graph-based distance calculation methods are available to you shortest.paths < / P>


Comments

Post a Comment

Popular posts from this blog

ios - How do I use CFArrayRef in Swift? -

I am using Objective-class in my Swift Project through a bridging header. The method signature looks like this: - Some cement (some type) some parameters; I started by getting an example of class, calling method, and storing the value: var myInstance = MyClassWithThatMethod (); Var cfArr = myInstance.someMethod (some value); Then try to get the value in the array: var valueInArrayThatIWant = CFArrayGetValueAtIndex (cfArr, 0); However, I get the error unmanaged & lt; Cfarray & gt; Not like 'CFArray' . Unmanaged & lt; Cfarray & gt; also means? I looked, but I do not need to change the array in a fast array (though it would be good). I need to be able to get value from the array. I also tried CFArray method to pass in a function: func doSomeStuffOnArray (myArray: NSArray) {} Although I get the same error when using: doSomeStuffOnArray (cfArr); // unmanaged & lt; CFArray & gt; I am using CFArray because I want to sto...
Read more

eclipse plugin - Run java code error: Workspace is closed -

To create an automated project, I created a plug-in project with the following dependencies: Org.eclipse.core.resources org.eclipse.equinox.registry org.eclipse.core.runtime and the following The Java class is located in the src folder: Package Examiner; Import org.eclipse.core.resources.IProject; Import org.eclipse.core.resources.IWorkspaceRoot; Import org.eclipse.core.resources.ResourcesPlugin; Import org.eclipse.core.runtime.CoreException; Import org.eclipse.core.runtime.IProgressMonitor; Import org.eclipse.core.runtime.NullProgressMonitor; Public class tes {public static zero main (string [] args) {// TODO auto generated method stub IProgressMonitor progress monitor = new NullProgressMonitor (); IWorkspaceRoot root = ResourcesPlugin.getWorkspace (). GetRoot (); Ipoject project = root.jetproject ("desired projectname"); Try {Project.create (progress monitor); Project.open (progressMonitor); } Grip (CoreException E) {// TODO Auto generated blocking block e....
Read more

c - Error on building source code in VC 6 -

When I create my project in VC6, the BScMKE reports the following error Linking .. Creating Library ... .. / ... / PB / Bean / Debub / ICEERP.Lib and Object ........ P.B. / Bean / DBube / ICREPlay.APP Browse information file ... BSMMK: ​​error can not be read from BK 1505: file '../../PB/bin/TmpOut/ICReplay/Debug/ICReplay_PB.bsc' Error executing bscmake.exe Can someone help me? 2 Reason for reference from: "ICReplay_PB .bsc "file will be corrupted The trunk of the file will be happening with the compiler running out of disk space or is being interrupted when creating .sbr file I think this answer will help you!
Read more