[BioC] Diffbind error

jluis.lavin at unavarra.es jluis.lavin at unavarra.es
Wed Dec 19 16:17:07 CET 2012


Dear list,

--I'm trying to use Diffbind package to analyze some samples from a ChIP
seq experiment.
--I created a .csv samplesheet with the following fields and which samples
to analyze:

SampleID Tissue	Factor	Condition Replicate Reads Control	Peaks

Chip_K27.1+ Neural K27 Responsive 1 Chip_K27.BED Input_K.BED  Chip_K27.BED

Chip_K27.2+ Neural K27 Responsive 1 Chip_K27d.BED Input_K.BED Chip_K27d.BED

--All my BED files have 6 lines and have been retrieved from MACS peak
caller.

--When Diffbind reads the sheet everything seems correct:

> chippy = dba(sampleSheet="Peaksets_sample_sheet.csv")

Chip_K27.1+ Neural K27 Responsive  1 raw
Chip_K27.2+ Neural K27 Responsive  1 raw

#it plots the correlation heatmap

> plot(chippy)

#But fails in this step after a while
> chippy = dba.count(chippy, minOverlap=3)

Error in read.table(fn, skip = skipnum) : no lines available in input

#Alternatively if I try this command
> chippy = dba.count(chippy)

Error in cond$counts : $ operator is invalid for atomic vectors
In addition: Warning message:
In mclapply(arglist, fn, ..., mc.preschedule = TRUE, mc.allow.recursive =
TRUE) :
  all scheduled cores encountered errors in user code

Would any of you shed any light on this issue?

Best wishes and thanks in advance

JL


> sessionInfo()
R version 2.15.1 (2012-06-22)
Platform: x86_64-redhat-linux-gnu (64-bit)

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C
 [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8
 [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8
 [7] LC_PAPER=C                 LC_NAME=C
 [9] LC_ADDRESS=C               LC_TELEPHONE=C
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C

attached base packages:
[1] parallel  stats     graphics  grDevices utils     datasets  methods
[8] base

other attached packages:
[1] DiffBind_1.4.0       Biobase_2.18.0       GenomicRanges_1.10.5
[4] IRanges_1.16.4       BiocGenerics_0.4.0   BiocInstaller_1.8.3

loaded via a namespace (and not attached):
 [1] amap_0.8-7         edgeR_3.0.7        gdata_2.12.0       gplots_2.11.0
 [5] gtools_2.7.0       limma_3.14.3       RColorBrewer_1.0-5 stats4_2.15.1
 [9] tools_2.15.1       zlibbioc_1.4.0


-- 
Dr. José Luis Lavín Trueba

Dpto. de Producción Agraria
Grupo de Genética y Microbiología
Universidad Pública de Navarra
31006 Pamplona
Navarra
SPAIN



More information about the Bioconductor mailing list