Question: error for toGRanges
0
gravatar for jolin0701-dy
3.0 years ago by
jolin0701-dy70
jolin0701-dy70 wrote:

I'd like to use ChIPpeakAnno to do downstream analysis after peak calling by macs2.

But I got some error for toGRanges to upload my results from macs2.

macs <- system.file("/Users/macs2", "/Users/macs2/cz3Rip_vs_cz1Inp_peaks.xls", package="ChIPpeakAnno")

macsOutput <- toGRanges(macs, format="MACS")

Error in read.table(data, header = header, ..., nrows = 5) : 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 also tried broadPeak, but still got the error.

Any suggestions?

Thanks so much~~~

chip-seq • 1.2k views
ADD COMMENTlink modified 3.0 years ago by igor7.6k • written 3.0 years ago by jolin0701-dy70
0
gravatar for igor
3.0 years ago by
igor7.6k
United States
igor7.6k wrote:

You don't need system.file if you know the path. Try:

macs <- "/Users/macs2/cz3Rip_vs_cz1Inp_peaks.xls"

Also, give a broadPeak, not xls file.

ADD COMMENTlink modified 3.0 years ago • written 3.0 years ago by igor7.6k
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 2.3.0
Traffic: 2106 users visited in the last hour