0
$\begingroup$

So I have a file, let's call it pcawg.tsv. It's formatted like the one down below.

Tumor sample ID  Normal sample ID  Mechanism
a                1                 deletion
b                2                 ecDNA
c                3                 ecDNA

I have a file of tumor sample ID's with the ecDNA mechanism called ecDNA_samples.tsv. It's formatted like this:

Tumor sample ID
a
b

To filter and extract the lines of pcawg.tsv that contain the tumor sample ID's in ecDNA_samples.tsv into a new file, I used this line of code below.

grep -f ecDNA_samples.tsv pcawg.tsv > ecDNA.vcf.tsv

Would this be the right code to use?

$\endgroup$
3
  • 2
    $\begingroup$ This has been asked on unix.stackexchange.com/questions/110645/… $\endgroup$ Aug 9, 2021 at 6:46
  • $\begingroup$ Is there a particular reason why you want to use grep to do this? It would be possible to make a more robust solution using R. $\endgroup$
    – gringer
    Aug 9, 2021 at 23:50
  • $\begingroup$ I kind of have to use this with Linux because I'm working with a Linux-based secure computing portal, since we have patient data which is confidential. $\endgroup$ Aug 10, 2021 at 0:16

2 Answers 2

1
$\begingroup$

This may help:

cat ecDNA_samples.tsv | xargs -I % grep '^%\s' pcawg.tsv

The command above executes grep '^<line>\s' pcawg.tsv, for each line of ecDNA_samples.tsv. Which matches the lines taken from ecDNA_samples.tsv to the sequences of characters between line start and a white-space in pcawg.tsv.

$\endgroup$
4
  • $\begingroup$ So I ran this code as cat ecDNA_samples.tsv | xargs -I % grep '^%\s' pcawg.tsv>ecDNA_normal.tsv, but when I opened the file ecDNA_normal.tsv in my directory, the file was empty. $\endgroup$ Aug 10, 2021 at 1:55
  • $\begingroup$ Use >> instead of > so an empty result at the end doesn't end up as the result instead of everything you've captured. $\endgroup$
    – Ram RS
    Aug 10, 2021 at 16:19
  • $\begingroup$ I see! Thank you! $\endgroup$ Aug 11, 2021 at 3:17
  • $\begingroup$ @RamRS It worked. $\endgroup$ Aug 13, 2021 at 1:40
0
$\begingroup$

Just in case an R solution is acceptable as well, this should be fairly robust code because it properly parses the TSV files and matches up column names:

#!/usr/bin/env Rscript

library(tidyverse)
read_tsv("pcawg.tsv") %>%
   inner_join(read_tsv("ecDNA_samples.tsv")) %>%
   write_tsv("filtered_pcawg.tsv")
$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.