robertaboukhalil / ginkgo

Cloud-based single-cell copy-number variation analysis tool
qb.cshl.edu/ginkgo
BSD 2-Clause "Simplified" License
47 stars 28 forks source link

problem in the R code #9

Closed RichardCorbett closed 7 years ago

RichardCorbett commented 7 years ago

Hi there,

I'm setting up a local install of ginko and running a set of just over 900 single cell libraries. After installing all the required libraries I had to make a small change to the code:

line 25 of analyze.sh is looking for a folder called "genomes" this isn't in my install so I took it out and got much farther.

Now I get to the part when lots of .jpegs are being created and I get an error that I have narrowed down to this location:

/opt/ginkgo/scripts/process.R /opt/ginkgo/genomes/hg19 /opt/ginkgo/uploads/CIC1 status.xml data 0 variable_100000_150_bwa ward.D2 euclidean 1 refDummy.bed_mapped 0 ploidyDummy.txt 1 0
Loading required package: amap
Loading required package: methods

Attaching package: ‘gplots’

The following object is masked from ‘package:stats’:

    lowess

Error in clip(tu[1], mean(temp) - (diff(reads$mids)/2), tu[3], tu[4]) : 
  invalid 'x2' argument
In addition: There were 50 or more warnings (use warnings() to see the first 50)
Execution halted

I get this error after getting though 102 libraries for which I get lots of pretty .jpeg images.

Looking in the status.xml file I am able to narrow this down to happening on just one of my libraries. That particular library doesn't looking very different from the others.

Do you have any suggestions I could try ?

tgarvin commented 7 years ago

Does that library have a very low number of mapped reads compared to the others?

Sent from my iPhone

On Aug 2, 2017, at 11:51 AM, Richard Corbett notifications@github.com<mailto:notifications@github.com> wrote:

Hi there,

I'm setting up a local install of ginko and running a set of just over 900 single cell libraries. After installing all the required libraries I had to make a small change to the code:

line 25 of analyze.sh is looking for a folder called "genomes" this isn't in my install so I took it out and got much farther.

Now I get to the part when lots of .jpegs are being created and I get an error that I have narrowed down to this location:

/opt/ginkgo/scripts/process.R /opt/ginkgo/genomes/hg19 /opt/ginkgo/uploads/CIC1 status.xml data 0 variable_100000_150_bwa ward.D2 euclidean 1 refDummy.bed_mapped 0 ploidyDummy.txt 1 0 Loading required package: amap Loading required package: methods

Attaching package: ‘gplots’

The following object is masked from ‘package:stats’:

lowess

Error in clip(tu[1], mean(temp) - (diff(reads$mids)/2), tu[3], tu[4]) : invalid 'x2' argument In addition: There were 50 or more warnings (use warnings() to see the first 50) Execution halted

I get this error after getting though 102 libraries for which I get lots of pretty .jpeg images.

Looking in the status.xml file I am able to narrow this down to happening on just one of my libraries. That particular library doesn't looking very different from the others.

Do you have any suggestions I could try ?

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHubhttps://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_robertaboukhalil_ginkgo_issues_9&d=DwMFaQ&c=mkpgQs82XaCKIwNV8b32dmVOmERqJe4bBOtF0CetP9Y&r=_00cSwveK_yUD35hWysVLQ&m=gBWOtwAcWWVfAs9dAorscAlMeQ9hVCWUEp7qX-WBE8A&s=xflLgPVKo2TsLOXaetsP5FtC-YFcPrKGIje3xVFiujA&e=, or mute the threadhttps://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_notifications_unsubscribe-2Dauth_AEHWMg9uU51OWK40DJd2GsETPzbMTHM1ks5sUMVFgaJpZM4Orgdx&d=DwMFaQ&c=mkpgQs82XaCKIwNV8b32dmVOmERqJe4bBOtF0CetP9Y&r=_00cSwveK_yUD35hWysVLQ&m=gBWOtwAcWWVfAs9dAorscAlMeQ9hVCWUEp7qX-WBE8A&s=VS5FNbda5G7hLOskwG4OV6gGqjdx5vsLZ_U3pU9zZ38&e=.

RichardCorbett commented 7 years ago

The BAM has roughly the same amount of reads as the others, but I just realized that the .bed.gz for that one is empty. I'm starting from scratch to see if this happens again.

tgarvin commented 7 years ago

Great let me know if you have any other questions.

Sent from my iPhone

On Aug 2, 2017, at 1:04 PM, Richard Corbett notifications@github.com<mailto:notifications@github.com> wrote:

The BAM has roughly the same amount of reads as the others, but I just realized that the .bed.gz for that one is empty. I'm starting from scratch to see if this happens again.

— You are receiving this because you commented. Reply to this email directly, view it on GitHubhttps://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_robertaboukhalil_ginkgo_issues_9-23issuecomment-2D319782614&d=DwMCaQ&c=mkpgQs82XaCKIwNV8b32dmVOmERqJe4bBOtF0CetP9Y&r=_00cSwveK_yUD35hWysVLQ&m=PDys-YcyMJpbN8pfInoEYxiL92QGYj_IJPCGAet3J4c&s=MGhZOKtL6ikZhOjqrkzZZwW9xe5O5HdacDok8fj5oZI&e=, or mute the threadhttps://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_notifications_unsubscribe-2Dauth_AEHWMuatvK9IgfUAHXecJHLj71JkKHP5ks5sUNZcgaJpZM4Orgdx&d=DwMCaQ&c=mkpgQs82XaCKIwNV8b32dmVOmERqJe4bBOtF0CetP9Y&r=_00cSwveK_yUD35hWysVLQ&m=PDys-YcyMJpbN8pfInoEYxiL92QGYj_IJPCGAet3J4c&s=xK9wOWrMImyvD7XFM4Gh0JjS9QGmpNkqCoOqaX6Qy9A&e=.

RichardCorbett commented 7 years ago

Restarting from scratch solved the problem. Thanks for the help.