SGFeatureCounts function in SGSeq package failed to generate splice graphs
1
0
Entering edit mode
mgdrnl • 0
@mgdrnl-14180
Last seen 4 months ago
Spain

Hello

I think this question is for the package developers, but maybe somebody had the same problem. I am traying to analyze 59 bam files in a computer with 110GB of available memory usage, all bam files occupy 149GB. I can't obtain the splice graph and I don't get any error message.. I monitorized the process and it looks like the memory is full for this process, so maybe that is the reason. I parallelize the process with 12 CPU, but still it is not enough. To obtain the SGFeatureCounts class, all bam files are downloaded to R first? If it is so, I guess I would need to increase the memory.. Could you give me any advice how can I complete the process in my computer?

Thank you,

Magda

sgseq • 1.1k views
ADD COMMENT
0
Entering edit mode
@leonard-goldstein-6845
Last seen 5 months ago
Australia

Hi Magda,

Can you please provide more details -- what commands were used to generate the SGFeatureCounts object and what package versions were used i.e. output of sessionInfo(). For larger analyses I strongly recommend not to use wrapper functions like analyzeFeatures() and instead run individual analysis steps separately (see section 12 of the vignette). This helps with troubleshooting and allows saving and inspecting of intermediate results. Please also see section 13 for recommendations on using multiple cores and typical memory requirements. 

Leonard

ADD COMMENT

Login before adding your answer.

Traffic: 726 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6