This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Sorting RNA-seq data

Hello,

I have some RNA-seq data on an Excel spreadsheet. My gene/s of interest follow a particular expression pattern. I would like to know if it is possible to sort the data/identify genes that match 3 conditions set by me. For example, I'd like to be able to see all the genes that match this expression pattern :

Cell line A < Cell Line B < Cell line C > Cell line D
rna-seq
• 2,620 views • link
updated by zx8754 zx8754 • written by agrisimo2 • agrisimo2

Please provide example input and expected output.

• 1 views • link
written by zx8754 zx8754

For example-

           Cell A         Cell B       Cell C     Cell D
Gene A    0.115175459   3.484635909 6.571842857 4.349035833
Gene B    0.021664012   2.939972182 3.448264286 3.8535915
Gene C    0.014484529   3.347903818 5.250840143 4.148886458
Gene D    0.0749899     33.82436091 52.07118571 30.74083333

The command should exclude gene B from the data set, since it does not follow the pattern A < B < C > D, and provide me with a list of genes that does so Gene A,C and D.

• 1 views • link
written by agrisimo2 • agrisimo2

I have some RNA-seq data on an Excel spreadsheet.

enter image description here

• 1 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

3 answers

If you are using Microsoft Excel,

=IF(AND(A1<B1,B1<C1,C1>D1),"YES","NO")
• 1 views • link
updated by zx8754 zx8754 • written by EagleEye EagleEye

With awk that would be: 

awk '($2 < $3 && $3 < $4 && $4>$5 )'
• 1 views • link
updated by zx8754 zx8754 • written by Pierre Lindenbaum Pierre Lindenbaum

Using R:

# read the file, something like:
df1 <- read.table("myFile.txt")

# then filter
df1[ df1$CellA < df1$CellB & df1$CellB < df1$CellC & df1$CellC > df1$CellD, ]
• 1 views • link
written by zx8754 zx8754

Log in to answer this question.

More posts like this