'Killed' error during RNAseq alignment
Entering edit mode
Sukeerthi • 0
Last seen 12 months ago

I am trying to align RNA sequencing data using the following code. However, it shows an error killed at the end. It does not specify the reason as well. I initially assumed memory could have been an issue, but the error seems to continue even when the memory issue was sorted in a new system. This is my code. Any help would be greatly appreciated.

myExons <- exons(TxDb.Hsapiens.UCSC.hg19.knownGene, columns = c("tx_id", "gene_id"))
myExons <- myExons[lengths(myExons$gene_id) == 1]

dfExons <- as.data.frame(myExons)
SAF <- data.frame(GeneID = dfExons$gene_id, Chr = dfExons$seqnames, Start = dfExons$start,
    End = dfExons$end, Strand = dfExons$strand)

This is the error message I am receiving

subjunc data.frame myExons Rsubread • 723 views
Entering edit mode

As you can see in the output, some columns in this GRanges object are character vectors, others are lists, so this cannot readily be parsed into a meaningful data.frame, that is why it takes long and crashes. Your code does not do any alignment, so why don't you state what the end goal is so one can suggest better code to avoid these issues. Do you want to make a tx2gene map for tximport?

Entering edit mode
Last seen 12 hours ago
WEHI, Melbourne, Australia

Neither your code nor the image show any use of Rsubread commands. You don't appear to be doing any alignment, just making a data.frame.

Anyway, there is no R or Bioconductor error called "Killed". "Killed" is an error arising from your operating system, perhaps because you've exceeded your memory quota. It cannot be caused by R commands, except indirectly by using memory.

The hg19 genome is very out of date. It would be far better and more efficient to use Rsubread's inbuilt hg38 annotation.


Login before adding your answer.

Traffic: 804 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6