Giter VIP home page Giter VIP logo

viber's Introduction

VIBER

R-CMD-check pkgdown Lifecycle: maturing

VIBER is a package that implements a variational Bayesian model to fit multi-variate Binomial mixtures. The statistical model is semi-parametric and fit by a variational mean-field approximation to the model posterior. The components are Binomial distributions which can model count data; these can be used to model sequencing counts in the context of cancer, for instance. The package implements methods to fit and visualize clustering results.

Citation

If you use VIBER, please cite:

  • G. Caravagna, T. Heide, M.J. Williams, L. Zapata, D. Nichol, K. Chkhaidze, W. Cross, G.D. Cresswell, B. Werner, A. Acar, L. Chesler, C.P. Barnes, G. Sanguinetti, T.A. Graham, A. Sottoriva. Subclonal reconstruction of tumors by using machine learning and population genetics. Nature Genetics 52, 898โ€“907 (2020).

Help and support

Installation

You can install the released version of VIBER from GitHub with:

# install.packages("devtools")
devtools::install_github("caravagnalab/VIBER")

Copyright and contacts

Giulio Caravagna. Cancer Data Science (CDS) Laboratory.

viber's People

Contributors

caravagn avatar l0ka avatar luca-dex avatar nicola-calonaci avatar

Stargazers

 avatar  avatar

Watchers

 avatar  avatar  avatar  avatar

Forkers

flywind2 l0ka

viber's Issues

"trace = TRUE" results in "variational_fit" error

Hi @caravagn, I noticed that when I set trace = TRUE in the variational_fit step it fails, resulting in this error:

[easypar] 10/10 computations returned errors and will be removed.
Error in sum(errs) : invalid 'type' (list) of argument

Specifically, the issue seems to be coming from here, the piece of code starting at row 286:

if(save_trace) {
      fit_trace = bind_rows(fit_trace,
                            tibble(
                              cluster.Binomial = labels,
                              step = i))
}

Replacing with cluster.Binomial = unlist(labels) seems to solve the problem but I'm not sure that's what you were looking for.
That's the code I used (with VIBER version: VIBER_0.1.0):

library(VIBER)
data(mvbmm_example)
variational_fit(x     = mvbmm_example$successes,
                y     = mvbmm_example$trials,
                trace = TRUE)

Thanks!

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    ๐Ÿ–– Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo D3

    Bring data to life with SVG, Canvas and HTML. ๐Ÿ“Š๐Ÿ“ˆ๐ŸŽ‰

Recommend Topics

  • javascript

    JavaScript (JS) is a lightweight interpreted programming language with first-class functions.

  • web

    Some thing interesting about web. New door for the world.

  • server

    A server is a program made to process requests and deliver data to clients.

  • Machine learning

    Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google โค๏ธ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.