-
Notifications
You must be signed in to change notification settings - Fork 22
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
long vectors not supported yet #32
Comments
The error is actually from R itself: https://github.com/wch/r-source/blob/trunk/src/include/Rinlinedfuns.h This used to be a more common error in R before version...3 maybe? There's possibly something we could do to fix this. We'll investigate. For context: Best, Evan |
Thanks for the issue! Are the three samples from the same individual (so that they have the same germline SNP profile)? If so, there's no need to merge the bams manually; You can supply multiple BAMs and barcode files to Best, |
Dear Teng, As for the 'long vectors error', it is strange since it is running with R 4.0.3, here the sessionInfo() for further details:
Matrix products: default locale: attached base packages: other attached packages: loaded via a namespace (and not attached): |
The error occurred because there were more than one individual's genotypes in the allele data. Only data from the same individual should be provided to Thanks, |
Ok, thanks for the info.
I thought something like that (cbind/rbind) would be the solution. As for
showing the posteriors of different samples combined on an integrated
object, that was my first working hypothesis, but I was not sure I could
then interpret the results well, as the posteriors for normal vs tumor are
intra-sample. Indeed, my point of analyzing them together was because in
one sample I had very few normal cells and I wondered whether running
phasing/numbat altogether could aid, but now I see it is not possible for
the individual germline in normal cells is key. I think I will then
visualize the single-dataset calls (for samples coming from different
individuals) on an integrated object.
Thanks again for a great package and help!
Jose
…--------------------------------------
Jose M. Garcia Manteiga PhD
Computational Biologist
Center for Translational Genomics and BioInformatics
Dibit2-Basilica, 4A3
San Raffaele Scientific Institute
Via Olgettina 58, 20132 Milano (MI), Italy
Tel: +39-02-2643-9211
e-mail: ***@***.***
Il giorno mar 17 mag 2022 alle ore 05:52 Teng Gao ***@***.***>
ha scritto:
Hi @josegarciamanteiga <https://github.com/josegarciamanteiga>,
The error occurred because there were more than one individual's genotypes
in the allele data. Only data from the same individual should be provided
to pileup_and_phase.R and run_numbat. If you have two samples from the
same individual, you can concatenate the gene count matrices (e.g. cbind)
and allele dataframes (e.g. rbind) as input to run_numbat. If the third
sample belongs to a separate individual, I would run it separately. If you
want to plot the single-cell posteriors in an integrated expression
embedding from different samples/individuals, you can combine the posterior
dataframes (e.g. nb$joint_post, nb$clone_post) after reading in the
results for each individual separately. For more info on the output, please
see this tutorial
<https://kharchenkolab.github.io/numbat/articles/visualization.html#single-cell-cnv-calls>
.
Thanks,
Teng
—
Reply to this email directly, view it on GitHub
<#32 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/AA2UOMIPYRKOCJCLJGYDQVDVKMJXFANCNFSM5V3YLCSQ>
.
You are receiving this because you were mentioned.Message ID:
***@***.***>
|
Hi,
Thanks for the package. Spectacular results with single cell RNASeq in tumors. I'd like to publish the identification of CAFs as normal cells in my tumors using it and it works smoothly in single datasets from 10X but I tried to mix three samples and got this error:
..../....
Retesting CNVs..
Retesting CNVs..
Retesting CNVs..
Retesting CNVs..
Retesting CNVs..
Finishing..
Finishing..
Finishing..
Finishing..
Finishing..
Error in vec_slice(x_out, x_slicer) :
long vectors not supported yet: ../../src/include/Rinlinedfuns.h:535
Error: Tibble columns must have compatible sizes.
3
.1
.2
.0
.ℹ Only values of size one are recycled.
Backtrace:
█
%>%
(...)%>%
(...)Warning message:
In mclapply(bulks %>% split(.$sample), mc.cores = ncores, function(bulk) { :
scheduled core 2 encountered error in user code, all values of the job will be affected
Execution halted
I used pileup_and_phase.R without problems on a bam merged from the cellranger bams where I substituted the "-1" at the end of the barcodes to avoid collisions after using cellranger aggr to generate barcodes.
The error is thrown by run_numbat run with 64GB and 12 cores.
Thanks for the help
Jose
The text was updated successfully, but these errors were encountered: