Hacktoberfest 2026: los issues que los mantenedores marcaron para octubre, abiertos y aptos para principiantes. Explorar issues de Hacktoberfest

Import of BigWig files with rtracklayer fails under MultiCoreParam

Abierto
#111 21 comentarios 0 reacciones 0 asignados Ver en GitHub

Nadie ha tomado este issue todavía.

Evaluación

Dificultad
4/5
Tiempo estimado
3-5 días
Aptitud para principiantes
35/100
Tipo de issue
Error
Claridad
Necesita aclaración
Estado de actividad
Estancado
Stack tecnológico
r
Área
distributed-systems

Línea de trabajo

Start by running the reproducible example with bplapply, comparing SerialParam and MulticoreParam around rtracklayer's import and empty BigWig tiles. Trace the interaction between MulticoreParam workers and import, then confirm the cause and define a regression check covering empty and non-empty regions.

Escrito por el modelo de indexación a partir del texto del issue.

Descripción

I'm seeing a weird bug in in my CAGEfightR package when trying to import BigWigFiles.

Here's a reproducible example:

# Load the example data
library(BiocParallel)
library(rtracklayer)
library(CAGEfightR)
data('exampleDesign')

# Example bigwig files
bw_plus <- system.file('extdata', exampleDesign$BigWigPlus,
                       package = 'CAGEfightR')
bw_plus <- BigWigFileList(bw_plus)

# Genome
mm9 <- SeqinfoForUCSCGenome("mm9")
grl <- tileGenome(mm9, ntile=3)

# Import, knowing that some areas will be empty
import(bw_plus[[1]], which=grl[[1]]) # empty
import(bw_plus[[1]], which=grl[[2]]) # empty
import(bw_plus[[1]], which=grl[[3]]) # has signal

# This works
lapply(bw_plus, import, which=grl[[1]])
lapply(bw_plus, import, which=grl[[2]])
lapply(bw_plus, import, which=grl[[3]])

# This works
register(SerialParam())
bplapply(bw_plus, import, which=grl[[1]])
bplapply(bw_plus, import, which=grl[[2]])
bplapply(bw_plus, import, which=grl[[3]])

# This works
register(MulticoreParam(3))
bplapply(bw_plus, import)

# This fails!
register(MulticoreParam(3))
bplapply(bw_plus, import, which=grl[[1]])
bplapply(bw_plus, import, which=grl[[2]])
bplapply(bw_plus, import, which=grl[[3]])

To summarize:

  • You can import entire BigWig-files using bplapply.
  • You can import part of a BigWig using bplapply and SerialParam
  • You cannot import an empty tile from a BigWig file with bplapply and MultiCoreParam.

What's causing this strange behavior?

Lenguaje dominante
R
Estrellas
69
Forks
32
Sin PR fusionados en 30 d

Guía de contribución

No hay ninguna guía de contribución indexada para este repositorio

Primeros pasos

  1. Lee el issue completo y luego la guía de contribución del proyecto.
  2. Comenta en el issue que vas a ocuparte — evita que dos personas hagan lo mismo.
  3. Haz un fork del repositorio y trabaja en una rama.
  4. Abre un pull request que haga referencia al número del issue.

Más de Bioconductor/BiocParallel

Todos los issues de Bioconductor/BiocParallel

Issues similares

Más issues de R

Recibe los nuevos issues en tu correo

Un resumen breve de issues de GitHub para principiantes.