Question: chippeakanno file reading error
0
gravatar for saddamhusain77
4 months ago by
saddamhusain770 wrote:

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 • 117 views
ADD COMMENTlink modified 4 months ago • written 4 months ago by saddamhusain770
Answer: chippeakanno file reading error
0
gravatar for Ou, Jianhong
4 months ago by
Ou, Jianhong1.1k
United States
Ou, Jianhong1.1k wrote:

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

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

ADD COMMENTlink written 4 months ago by Ou, Jianhong1.1k
Answer: chippeakanno file reading error
0
gravatar for saddamhusain77
4 months ago by
saddamhusain770 wrote:

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 COMMENTlink modified 4 months ago • written 4 months ago by saddamhusain770

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 REPLYlink written 4 months ago by Ou, Jianhong1.1k
Answer: chippeakanno file reading error
0
gravatar for saddamhusain77
4 months ago by
saddamhusain770 wrote:

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 COMMENTlink written 4 months ago by saddamhusain770
Please log in to add an answer.

Help
Access

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