This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to get binary data from gene set data?

Hi, I want to melt a geneset data to get binary interactions. How can I get this? I have mentioned below the input and output.

Thank you,

Input:
AANWWTGC_UNKNOWN    4208    481 6095
AAAYRNCTG_UNKNOWN   4052    9842    23047
AAAYWAACM_HFH4_01   4208    161753  3151

Output:
AANWWTGC_UNKNOWN    4208
AANWWTGC_UNKNOWN    481
AANWWTGC_UNKNOWN    6095
AAAYRNCTG_UNKNOWN   4052
AAAYRNCTG_UNKNOWN   9842
AAAYRNCTG_UNKNOWN   23047
AAAYWAACM_HFH4_01   4208
AAAYWAACM_HFH4_02   161753
AAAYWAACM_HFH4_03   3151
rna-seq r
• 1,251 views • link
updated by Kevin Blighe Kevin Blighe • written by Neu • Neu

Do all the rows of your input have four columns? And are there column names?

• 0 views • link
written by Sean Davis Sean Davis

1 answer

Try the melt() function in reshape:

MyDataFrame
                 V1   V2     V3    V4
1  AANWWTGC_UNKNOWN 4208    481  6095
2 AAAYRNCTG_UNKNOWN 4052   9842 23047
3 AAAYWAACM_HFH4_01 4208 161753  3151

require(reshape)

melt(MyDataFrame, id.vars=colnames(MyDataFrame)[1])
                 V1 variable  value
1  AANWWTGC_UNKNOWN       V2   4208
2 AAAYRNCTG_UNKNOWN       V2   4052
3 AAAYWAACM_HFH4_01       V2   4208
4  AANWWTGC_UNKNOWN       V3    481
5 AAAYRNCTG_UNKNOWN       V3   9842
6 AAAYWAACM_HFH4_01       V3 161753
7  AANWWTGC_UNKNOWN       V4   6095
8 AAAYRNCTG_UNKNOWN       V4  23047
9 AAAYWAACM_HFH4_01       V4   3151
• 0 views • link
written by Kevin Blighe Kevin Blighe

Log in to answer this question.

More posts like this