Wedge-lab / dpclust

Dirichlet Process based methods for subclonal reconstruction of tumours
GNU Affero General Public License v3.0
28 stars 16 forks source link

Error Estimating density between pairs of samples #20

Closed mdurante1 closed 1 year ago

mdurante1 commented 1 year ago

Hi,

I recently encountered an error during the estimating density step between pairs of samples step (see error output below) when running one of our multi-sample cases. Any insight into what could be occurring would be greatly appreciated. Of note, there’s good overlap of mutations between the 2 samples.

IMG_9093

I’m happy to provide any additional information requested. Thank you!

[1] "Loading data..."
[1] "Removed 0 with missing WTCount"
[1] "Removed 0 with missing mutCount"
[1] "Removed 610 with missing totalCopyNumber"
[1] "Removed 0 with missing copyNumberAdjustment"
[1] "Removed 476 with missing kappa"
[1] "Removed 0 with no coverage"
[1] "Removed 0 on not supported genomic regions"
[1] ""
[1] "Removed 610 mutations with missing data"
[1] ""
[1] "Running DPClust..."
[1] "100  /  2000"
[1] "200  /  2000"
[1] "300  /  2000"
[1] "400  /  2000"
[1] "500  /  2000"
[1] "600  /  2000"
[1] "700  /  2000"
[1] "800  /  2000"
[1] "900  /  2000"
[1] "1000  /  2000"
[1] "1100  /  2000"
[1] "1200  /  2000"
[1] "1300  /  2000"
[1] "1400  /  2000"
[1] "1500  /  2000"
[1] "1600  /  2000"
[1] "1700  /  2000"
[1] "1800  /  2000"
[1] "1900  /  2000"
[1] "2000  /  2000"
[1] "Estimating density between pairs of samples..."
[1] "Samples P2510T and P3298T"
Error in if (is.na(indices)) { : the condition has length > 1
Calls: RunDP ... DirichletProcessClustering -> Gibbs.subclone.density.est
Execution halted
mdurante1 commented 1 year ago

Issue was solved by going back to original install on R 4.1.3. This error was only re-produced with R 4.2.3.