Problem importing 5C data
Entering edit mode
gugpuccio ▴ 10
Last seen 21 months ago

Dear All

I'm sorry if this question was answered other times, and if it's a stupid one, but i can't understand how to use the search bar in this site. Every time i look for a word it gives out empy search.

The problem:

I'm trying to import a 5C matrix with the HiTC package in R, but i'm unable to do it. This is the line of code i used:

M=importC(con = "./Desktop/ENR-HeLa-R2.csv", xgi.bed = "./Desktop/ENr.hg19.reverses", ygi.bed = "./Desktop/ENr.hg19.forwards")

And the error message is:

Genomic intervals loaded
Error in HTClist(unlist(obj)) : 
  all elements in '...' must be HTCexp objects
In addition: Warning message:
In mclapply(chromPair, function(chr) { :
  all scheduled cores encountered errors in user code

Then i tried to follow the import tutorial on the HiTC manual, trying to export and then re-importing the data:

exportC(E14.bin, file="E14")
importC(con="E14.mat", xgi.bed ="E14_xgi.bed", ygi.bed ="E14_ygi.bed")

Getting this error:

Genomic intervals loaded
Loading chrX-chrX...
Error in isN0(as(m, "matrix")) : 
  (list) object cannot be coerced to type 'double'

in this case i used "E14.mat" instead of "E14.count" as written in the tutorial, because it was the only output of the export function, together with the .bed files.

I don't understand what i'm doing wrong. I'm sorry if my question looks stupid, but i'm a newbie in R programming.
Thanks a lot in advance.

importing HiTC • 768 views

Login before adding your answer.

Traffic: 370 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6