Hacktoberfest 2026: as issues que os mantenedores marcaram para outubro, abertas e boas para iniciantes. Ver issues do Hacktoberfest

`expand()` and non-standard rowRanges columns

Aberta
  • Algo errado?
#85 1 comentário 0 reações 0 responsáveis Ver no GitHub

@jmg421 já está trabalhando nisso.

Desde 30/6/2026.

  • #95 de @jmg421 — aberto

Avaliação

Dificuldade
3/5
Tempo estimado
1-2 dias
Facilidade para iniciantes
45/100
Tipo de issue
Funcionalidade
Clareza
Razoavelmente clara
Status de atividade
Estagnada
Stack de tecnologia
r
Domínio
bioinformatics

Direção de pesquisa

Reproduza o exemplo em R fornecido com VariantAnnotation e compare rowRanges(vcf) antes e depois de expand(vcf). Comece localizando a implementação de expand() para CollapsedVCF; o trabalho estará concluído quando o objeto expandido mantiver colunas não padrão de rowRanges, como SNP_name e num_alts, para cada linha resultante.

Escrita pelo modelo de indexação a partir do texto da issue.

Descrição

hi there,

is there a way to retain non-standard rowRanges columns when I use expand() on a CollapsedVCF object? Some demo code below.

If not, please can you consider this an enhancement request? thanks!

all the best,

Janet

library(VariantAnnotation)

vcf <- VCF(rowRanges = GRanges("chr1", IRanges(1:4*3, width=c(1, 2, 1, 1))))
alt(vcf) <- DNAStringSetList("A", c("TT"), c("G", "A"), c("TT", "C"))
ref(vcf) <- DNAStringSet(c("G", c("AA"), "T", "G"))

## add some non-standard columns to rowRanges
mcols(rowRanges(vcf))$SNP_name <- paste("SNP_", 1:length(vcf), sep="")
mcols(rowRanges(vcf))$num_alts <- elementNROWS(alt(vcf))

## take a look
rowRanges(vcf)

GRanges object with 4 ranges and 6 metadata columns:
      seqnames    ranges strand |    SNP_name  num_alts            REF                ALT      QUAL
         <Rle> <IRanges>  <Rle> | <character> <integer> <DNAStringSet> <DNAStringSetList> <numeric>
  [1]     chr1         3      * |       SNP_1         1              G                  A        NA
  [2]     chr1       6-7      * |       SNP_2         1             AA                 TT        NA
  [3]     chr1         9      * |       SNP_3         2              T                G,A        NA
  [4]     chr1        12      * |       SNP_4         2              G               TT,C        NA
           FILTER
      <character>
  [1]        <NA>
  [2]        <NA>
  [3]        <NA>
  [4]        <NA>
  -------
  seqinfo: 1 sequence from an unspecified genome; no seqlengths


## those extra columns are lost when we expand
vcfLong <- expand(vcf)
rowRanges(vcfLong)
rowRanges(vcfLong)

GRanges object with 6 ranges and 4 metadata columns:
      seqnames    ranges strand |            REF            ALT      QUAL      FILTER
         <Rle> <IRanges>  <Rle> | <DNAStringSet> <DNAStringSet> <numeric> <character>
  [1]     chr1         3      * |              G              A        NA        <NA>
  [2]     chr1       6-7      * |             AA             TT        NA        <NA>
  [3]     chr1         9      * |              T              G        NA        <NA>
  [4]     chr1         9      * |              T              A        NA        <NA>
  [5]     chr1        12      * |              G             TT        NA        <NA>
  [6]     chr1        12      * |              G              C        NA        <NA>
  -------
  seqinfo: 1 sequence from an unspecified genome; no seqlengths
Linguagem predominante
R
Estrelas
32
Forks
21
Nenhum PR com merge em 30d

Preparar o ambiente

Este projeto não oferece contêiner de desenvolvimento, Dockerfile nem guia de contribuição, então a configuração fica por sua conta: comece pelo README e veja nosso guia da primeira contribuição para os passos gerais.

Primeiros passos

  1. Leia a issue inteira e depois o guia de contribuição do projeto.
  2. Comente na issue dizendo que vai assumir — evita que duas pessoas façam o mesmo trabalho.
  3. Faça um fork do repositório e trabalhe em uma branch.
  4. Abra um pull request que referencie o número da issue.

Mais de Bioconductor/VariantAnnotation

Todas as issues de Bioconductor/VariantAnnotation

Issues semelhantes

Mais issues de R

Receba novas issues na sua caixa de entrada

Um resumo curto de issues do GitHub para quem está começando.