chippeakanno file reading error
3
0
Entering edit mode
@saddamhusain77-19420
Last seen 2.4 years ago

Hi, I want to check peak overlap between two chip-seq datasets using ChiPpeakAnno package. When I try to read the MACS output file with the code provided in the vignette it shows the following error.

brg1bed <- system.file("extdata", "brg1.bed", package="ChIPpeakAnno")

grbrg1 <- toGRanges(brg1bed, format="BED", header=FALSE)

Error in read.table(data, header = header, comment.char = comment.char, :

no lines available in input

In addition: Warning message:

In file(file, "rt") :

file("") only supports open = "w+" and open = "w+b": using the former

I have replaced the "extdata" to my working directory but it still showed error.

When I read MACS output with read.delim it showed the same error

sl1bed <- system.file("chip brg1", "sl1.txt", package="ChIPpeakAnno")

grsl1 <- toGRanges(sl1bed, format="MACS", header=FALSE)

Error in read.table(data, header = header, comment.char = comment.char, :

no lines available in input

In addition: Warning message:

In file(file, "rt") :

file("") only supports open = "w+" and open = "w+b": using the former.

Any suggestion would be appreciated.

ChiPpeakAnno • 388 views
ADD COMMENT
0
Entering edit mode
Ou, Jianhong ★ 1.2k
@ou-jianhong-4539
Last seen 9 weeks ago
United States

suppose your brg1.bed is in current working directory, please try to run:

grbrg1 <- toGRanges("brg1.bed", format="BED", header=FALSE)

ADD COMMENT
0
Entering edit mode
@saddamhusain77-19420
Last seen 2.4 years ago

Thank you for your response, Jianhong. I tried what you suggested and got this error, even thought my working directory does contain brg1.bed file

brg1bed <- system.file("extdata", "brg1.bed", package="ChIPpeakAnno")

grbrg1 <- toGRanges("brg1bed", format="BED", header=FALSE)

Error in file(file, "rt") : cannot open the connection

In addition: Warning message:

In file(file, "rt") : cannot open file 'brg1bed': No such file or directory.

Could you please tell me what is the "extdata" in the code? is it a directory within the package where I am importing the data?

ADD COMMENT
0
Entering edit mode

you have typo, change brg1bed to brg1.bed for more information please refer:https://cran.r-project.org/doc/contrib/Paradis-rdebuts_en.pdf

ADD REPLY
0
Entering edit mode
@saddamhusain77-19420
Last seen 2.4 years ago

Sure, I will go through it. Indeed it was a typo error as changing brg1bed to brg1.bed worked out. Thank you for your time.

ADD COMMENT

Login before adding your answer.

Traffic: 321 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6