[BioC] edgeR: Problem with subsetting a DGEList in latest package version
Steve Lianoglou
lianoglou.steve at gene.com
Wed Aug 13 23:58:48 CEST 2014
Hi,
On Wed, Aug 13, 2014 at 2:37 PM, Katja Hebestreit <katjah at stanford.edu> wrote:
> Hello edgeR maintainers,
>
> I just ran into an error with code that is several weeks old. I cannot subset my DEGList object anymore. I created it from the HTseq output:
>
>
>> files <- list.files("Gene_counts",
> recursive = FALSE, pattern = "counts",
> full.names=TRUE)
>
>> rc <- readDGE(files, labels = labels)
>> dim(rc)
> [1] 23429 40
>
>> out1 <- which(colnames(rc) == "PGE_old_E")
>> rc <- rc[,-out1]
> Error in `$<-.data.frame`(`*tmp*`, "group", value = integer(0)) :
> replacement has 0 rows, data has 39
Can please you show us the output of:
R> lapply(rc, dim)
Thanks,
-steve
--
Steve Lianoglou
Computational Biologist
Genentech
More information about the Bioconductor
mailing list