Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
DESeq2 error in data.frame (multiple treatments and multiple replicates) KYR Bioinformatics 3 03-19-2019 10:05 AM
ExomeCNV - R - data.frame error sophiespo Bioinformatics 4 02-25-2016 01:16 AM
Matching a column from a data frame on the columns of another data frame zisis86 Bioinformatics 2 05-29-2014 05:32 AM
easyRNAseq error sylvian Bioinformatics 0 06-27-2013 06:33 PM
Solution for easyRNASeq Error in .doBasicCount(obj) weasteam Bioinformatics 0 05-02-2013 03:28 PM

Thread Tools
Old 11-18-2014, 02:41 PM   #1
Junior Member
Location: Penn

Join Date: Jun 2013
Posts: 4
Default easyRNAseq problem: Error in

easyRNAseq terminates with the error:

Error in, ...) :
no rows to aggregate

I use easyRNAseq to enter my data for DESeq:

> library(easyRNASeq)
> library(BSgenome.Hsapiens.UCSC.hg19)
> setwd('/home/sylvian/SerumFedStarved/DESeq/')
> chr.sizes=seqlengths(Hsapiens)
> bamfiles=dir(getwd(), pattern='*.bam$')
> conditions=c('fed','starved','fed','starved','fed','starved')
> names(conditions)<-c('DRSNpilot1_STARhg19_paired_sorted.bam','DRSNpilot2_STARhg19_paired_sorted.bam','TQ24_STARhg19_paired_sorted.bam','TQ25_STARhg19_paired_sorted.bam','TQ28_STARhg19_paired_sorted.bam','TQ29_STARhg19_paired_sorted.bam')
> rnaseqgenes<-easyRNASeq(filesDirectory=getwd(),
+ organism='Hsapiens',
+ chr.sizes=chr.sizes,
+ readLength=100L,
+ annotationMethod='biomaRt',
+ format='bam',
+ count='genes',
+ summarization='geneModels',
+ filenames=bamfiles,
+ conditions=conditions,
+ normalization=TRUE,
+ fitType='local',
+ outputFormat='DESeq')
Checking arguments...
Fetching annotations...
Computing gene models...
Summarizing counts...
Processing DRSNpilot1_STARhg19_paired_sorted.bam
Error in, ...) :
no rows to aggregate

Has anybody by any chance had the same problem and knows how to fix this?

I used to do
but this now produces an error:
Error in easyRNASeq(filesDirectory = getwd(), organism = "Hsapiens", chr.sizes = chr.sizes, :
chr.sizes should be a named list of integers. 'Use 'as.integer' to convert from numeric.

chr.sizes=seqlengths(Hsapiens) produces

chr1 chr2 chr3
249250621 243199373 198022430
chr4 chr5 chr6
191154276 180915260 171115067
chr7 chr8 chr9
159138663 146364022 141213431

Thanks a lot
sylvian is offline   Reply With Quote

deseq, easyrnaseq

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 09:32 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO