This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Error in DESeq analysis directory argument - No such file or directory

Hello,

I'm using DESeq2 for RNA-Seq data analysis. Everything looks great but I get an error about my directory at the final step even though I'm in the correct working directory

sampleCondition <- read.table(file=sampleFiles, header = TRUE, sep = "", row.names = 1)

sampleTable <- data.frame(sampleName = sampleFiles, fileName = sampleFiles, condition = samplecondition)

> getwd()
[1] "C:/Users/xyz/Downloads"

directory <- "/Users/xyz/Downloads/FeatureCounts.txt"
 ddsHTSeq <- DESeqDataSetFromHTSeqCount(sampleTable = sampleTable,directory = directory, design = ~condition)

Error in file(file, "rt") : cannot open the connection
In addition: Warning message:
In file(file, "rt") :
  cannot open file '/Users/xyz/Downloads/FeatureCounts.txt/0': No such file or directory

please help! I'm pretty sure the directory exists and that I am currently working in that directory (see getwd). Everything else in the analysis looks good

rna-seq deseq hts
• 2,292 views • link
written by AW • AW

If sampleFiles are in the directory: "C:/Users/xyz/Downloads", try directory <- "C:\\Users\\xyz\\Downloads"?

• 0 views • link
written by AK AK

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • Error in DESeqDataSetFromHTSeqCount
    written by abhay.kanodia • abhay.kanodia
    I am using HTseq pipeline for DESeq2: Directory = "/Users/abhaykanodia/Desktop/smallRNA/" condition = c("WT1", "WT2", "WT3", "NTC1", "NTC2", "NTC3) sampleFiles= c("AK1a_counts.txt","AK2a_counts.txt","AK3a_counts.txt"," AK4a_counts.txt","AK5a_counts.txt","AK6a_counts.txt") sampleName = c("AK1", "AK2", …
  • Specifying unequal replicates in deseq2 with htseq-count input
    written by kousi31 • kousi31
    I have two controls and four treatments. I calculated gene-counts using htseq. directory <- "~/d/new_de/gene_counts/" sampleFiles <- c("P-1.deseqfile", "P-2.deseqfile", "P-3.deseqfile", "P-4.deseqfile", "NP-3.deseqfile", "NP-4.deseqfile" ) sampleNames …
  • deseq2 replicates error
    written by Balatheskulter • Balatheskulter
    > setwd("C:/Users/SKULTER/Desktop/New folder/") > directory <- 'C:/Users/SKULTER/Desktop/New folder/' > samplefiles <- grep("SRR",list.files(directory),value = TRUE) > condition <- c('SRR5500529.txt','SRR5500530.txt') > sampletable <- data.frame(samplname = samplefiles,filename = …
  • DESeq2 compare all levels using wald and lrt test
    written by aiswaryabioinfo • aiswaryabioinfo
    I am currently working in plant pathogen interaction and is trying to make a gene co-expression network using RNA-seq data. I used HTseq-count to get …
  • DESeq2 - HTSeq-count txt error
    written by jsl • jsl
    Hi, I am trying to perform differential analysis on my HT-seq-counts .txt files, using a command below's command, sampleFiles <- grep("CD39",list.files(directory),value=TRUE) sampleCondition <- sub("(.*annotated...).*","",sampleFiles) sampleTable …
  • Error when running ddsHTSeq command in DESeq2
    written by Gabrielle.corso • Gabrielle.corso
    Hi there- I am trying to analyze RNA seq counts in R and am receiving an error code when running this chunk of code: # …
  • setting Coldata levels in Deseq2
    written by AP • AP
    Hello everyone, I am trying to see differential expression of a fungal genome with and without treatment of a chemical gossypol. I have one sample …
  • How to give directory path in DESeq2
    written by AP • AP
    Hello everyone, I did my transcripts count using Htseq and now I am trying to use DESeq2 package in R to see differential expression of …
  • Strange MA-plot using DESeq2
    written by stan stan
    Hello, I am using DESeq2 to perform some differential expression analysis between 3 levels of a disease condition. However, I am getting a strange looking …
  • DESeq2 model design: dose and time effect
    written by CuriusScientist CuriusScientist
    I have three samples each with 3 biological replicates i.e. one control and other at dose 12ug/ml and 25ug/ml at 24hr ``` > library('DESeq2') > …