Question: qCount in QuasR for RNAseq quantification
gravatar for tarek.mohamed
4.7 years ago by
tarek.mohamed270 wrote:

Hi All,

I am using QuasR to analyze mt paired end RNAseq data, now I have bam files for my 12 samples and the next step is to generate a count table for Differential gene expression analysis. I follwed the QuasR vignette and I used qCount() but I had a problem (R could not find my genome file)

sampleFile <- "DOX_RNAseq.txt"
proj_RNASeq<- qAlign(sampleFile, genomeFile,cacheDir = "E:/Doxorubicin Project RNA-seq Data")

annotationFile <- "BSgenome.Hsapiens.NCBI.GRCh38_annotation.gtf"
chrLen <- scanFaIndex(genomeFile)

Error in scanFaIndex(open(FaFile(file))) : 
  error in evaluating the argument 'file' in selecting a method for function 'scanFaIndex': Error in .io_check_exists(path(con)) : file(s) do not exist:


myposts • 1.3k views
ADD COMMENTlink written 4.7 years ago by tarek.mohamed270
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1387 users visited in the last hour