Code Monkey home page Code Monkey logo

dyno's Introduction

Build Status codecov AppVeyor Build Status โ„น๏ธ Tutorials ย  ย  โ„น๏ธ Reference documentation

Inferring trajectories using dyno

The dyno package offers end-users a complete TI pipeline. It features:

For information on how to use dyno, check out the installation instructions, tutorials and documentation at dynverse.org

Latest changes

Check out news(package = "dyno") or NEWS.md for a full list of changes.

Recent changes in dyno 0.1.1 (06-04-2019)

  • BUGFIX: Fixed example dataset

  • MINOR CHANGE: Bumping version requirement of dynmethods to 1.0.1.

Recent changes in dyno 0.1.0 (29-03-2019)

  • Initial beta release of dyno

Dynverse dependencies

dyno's People

Contributors

helena-todd avatar kou avatar rcannood avatar zouter avatar

Stargazers

 avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar

Watchers

 avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar

dyno's Issues

questions not found

Originally posted by @ysu2015 in #9 (comment)

Hi @zouter,

when I run the commands in the tutorial, I got another error.

library(dyno)
library(tidyverse)
data("fibroblast_reprogramming_treutlein")
dim(fibroblast_reprogramming_treutlein$expression)
task <- wrap_expression(
  counts = fibroblast_reprogramming_treutlein$counts,
  expression = fibroblast_reprogramming_treutlein$expression
)
guidelines <- guidelines_shiny(task)

Error in setdiff(names(questions), given_question_ids) : 
  object 'questions' not found

Error with pCreode = TypeError: slice indices must be integers or None or have an __index__ method

Hi Dynverse Team,

First and foremost, you created sth formidable here - I am a huge fan of your platform!! Up to now I could resolve a couple of issues by myself, but this one is beyond my horizon:

model.tree.pcreode <- infer_trajectory(dataset, ti_pcreode(n_pca_components = 10, num_runs = 10L), seed = 7, verbose = TRUE)
docker.io/dynverse/ti_pcreode:v0.9.9.01
Executing 'pcreode' on '20190808_144935__data_wrapper__GQCmdkdhU4' With parameters: list(n_pca_components = 10, num_runs = 10L),
inputs: expression, and
priors : 
Input saved to C:\Users\FUCHSD~1\AppData\Local\Temp\Rtmpslaedt\file2eac2c0f391d/ti
Running method using babelwhale
Running "C:\PROGRA~1\Docker\Docker\RESOUR~1\bin\docker.exe" run --name 20190808_154412__container__DJeLusmOdl -e "TMPDIR=/tmp2" \
  --workdir /ti/workspace -v "/c/Users/FUCHSD~1/AppData/Local/Temp/Rtmpslaedt/file2eac2c0f391d/ti:/ti" -v \
  "/c/Users/FUCHSD~1/AppData/Local/Temp/Rtmpslaedt/file2eac1ad32718/tmp:/tmp2" "dynverse/ti_pcreode:v0.9.9.01" --dataset /ti/input.h5 \
  --output /ti/output.h5
Traceback (most recent call last):
  File "/code/run.py", line 34, in <module>
    pca_reduced_data = data_pca.pca_set_components(min(parameters["n_pca_components"],expression.shape[1]))
  File "/pCreode/pcreode/pcreode.py", line 83, in pca_set_components
    return( self.pca[:,:n_components])
TypeError: slice indices must be integers or None or have an __index__ method
Traceback (most recent call last):
  File "/code/run.py", line 34, in <module>
    pca_reduced_data = data_pca.pca_set_components(min(parameters["n_pca_components"],expression.shape[1]))
  File "/pCreode/pcreode/pcreode.py", line 83, in pca_set_components
    return( self.pca[:,:n_components])
TypeError: slice indices must be integers or None or have an __index__ method
Error: Error during trajectory inference, see output above <U+2191><U+2191><U+2191>

Here is my dataset:

> str(dataset)
List of 11
 $ id                  : chr "20190808_144935__data_wrapper__GQCmdkdhU4"
 $ cell_ids            : chr [1:1654] "BM_1_GTATCTTAGGGCTTGA.1" "BM_1_TCATTTGTCTTATCTG.1" "BM_1_GATCGATGTAACGACG.1" "BM_1_ATCCGAACAGTGACAG.1" ...
 $ cell_info           : NULL
 $ counts              :Formal class 'dgCMatrix' [package "Matrix"] with 6 slots
  .. ..@ i       : int [1:2322752] 12 16 104 114 120 122 130 152 159 171 ...
  .. ..@ p       : int [1:10442] 0 42 45 427 434 601 602 606 653 981 ...
  .. ..@ Dim     : int [1:2] 1654 10441
  .. ..@ Dimnames:List of 2
  .. .. ..$ : chr [1:1654] "BM_1_GTATCTTAGGGCTTGA.1" "BM_1_TCATTTGTCTTATCTG.1" "BM_1_GATCGATGTAACGACG.1" "BM_1_ATCCGAACAGTGACAG.1" ...
  .. .. ..$ : chr [1:10441] "A2M" "A3GALT2" "AAAS" "AACS" ...
  .. ..@ x       : num [1:2322752] 1 1 1 1 1 2 1 2 1 1 ...
  .. ..@ factors : list()
 $ expression          :Formal class 'dgCMatrix' [package "Matrix"] with 6 slots
  .. ..@ i       : int [1:2322752] 12 16 104 114 120 122 130 152 159 171 ...
  .. ..@ p       : int [1:10442] 0 42 45 427 434 601 602 606 653 981 ...
  .. ..@ Dim     : int [1:2] 1654 10441
  .. ..@ Dimnames:List of 2
  .. .. ..$ : chr [1:1654] "BM_1_GTATCTTAGGGCTTGA.1" "BM_1_TCATTTGTCTTATCTG.1" "BM_1_GATCGATGTAACGACG.1" "BM_1_ATCCGAACAGTGACAG.1" ...
  .. .. ..$ : chr [1:10441] "A2M" "A3GALT2" "AAAS" "AACS" ...
  .. ..@ x       : num [1:2322752] 1 1 1 1 1 2 1 2 1 1 ...
  .. ..@ factors : list()
 $ expression_projected: NULL
 $ feature_info        :Classes โ€˜tbl_dfโ€™, โ€˜tblโ€™ and 'data.frame':	10441 obs. of  1 variable:
  ..$ feature_id: chr [1:10441] "A2M" "A3GALT2" "AAAS" "AACS" ...
 $ feature_ids         : chr [1:10441] "A2M" "A3GALT2" "AAAS" "AACS" ...
 $ prior_information   :List of 5
  ..$ start_id : chr [1:100] "BM_1_GTATCTTAGGGCTTGA.1" "BM_1_TCATTTGTCTTATCTG.1" "BM_1_GATCGATGTAACGACG.1" "BM_1_ATCCGAACAGTGACAG.1" ...
  ..$ end_id   : chr [1:500] "PB_TACGGATGTAGAGCTG.1" "PB_ACCCACTTCCGCATCT.1" "BM_1_GACCAATAGCTACCTA.1" "BM_1_GACCAATGTGCAGACA.1" ...
  ..$ groups_id:'data.frame':	1654 obs. of  2 variables:
  .. ..$ cell_id : chr [1:1654] "BM_1_GTATCTTAGGGCTTGA.1" "BM_1_TCATTTGTCTTATCTG.1" "BM_1_GATCGATGTAACGACG.1" "BM_1_ATCCGAACAGTGACAG.1" ...
  .. ..$ group_id: chr [1:1654] "HSC" "HSC" "HSC" "HSC" ...
  ..$ start_n  : num 1
  ..$ end_n    : num 5
 $ group_ids           : chr [1:19] "HSC" "MPP" "BCP" "BC" ...
 $ grouping            : Named chr [1:1654] "HSC" "HSC" "HSC" "HSC" ...
  ..- attr(*, "names")= chr [1:1654] "BM_1_GTATCTTAGGGCTTGA.1" "BM_1_TCATTTGTCTTATCTG.1" "BM_1_GATCGATGTAACGACG.1" "BM_1_ATCCGAACAGTGACAG.1" ... - attr(*, "class")= chr [1:5] "dynwrap::with_grouping" "dynwrap::with_prior" "dynwrap::with_expression" "dynwrap::data_wrapper" ...
> sessionInfo()
R version 3.6.1 (2019-07-05)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows 10 x64 (build 17763)

Matrix products: default

Random number generation:
 RNG:     Mersenne-Twister 
 Normal:  Inversion 
 Sample:  Rounding 
 
locale:
[1] LC_COLLATE=English_United States.1252  LC_CTYPE=English_United States.1252    LC_MONETARY=English_United States.1252
[4] LC_NUMERIC=C                           LC_TIME=English_United States.1252    

attached base packages:
[1] parallel  stats4    stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
 [1] dyno_0.1.1                  dynwrap_1.1.4.9000          dynplot_1.0.2.9000          dynmethods_1.0.5           
 [5] dynguidelines_1.0.0         dynfeature_1.0.0.9000       scran_1.12.1                scater_1.12.2              
 [9] futile.matrix_1.2.7         Matrix_1.2-18               data.table_1.12.2           gplots_3.0.3               
[13] Seurat_3.0.2                forcats_0.4.0               stringr_1.4.0               dplyr_0.8.3                
[17] purrr_0.3.2                 readr_1.3.1                 tidyr_0.8.3                 tibble_2.1.3               
[21] ggplot2_3.2.0               tidyverse_1.2.1             SingleCellExperiment_1.6.0  SummarizedExperiment_1.14.1
[25] DelayedArray_0.10.0         BiocParallel_1.18.0         matrixStats_0.54.0          Biobase_2.44.0             
[29] GenomicRanges_1.36.0        GenomeInfoDb_1.20.0         IRanges_2.18.1              S4Vectors_0.22.0           
[33] BiocGenerics_0.30.0        

loaded via a namespace (and not attached):
  [1] reticulate_1.13          R.utils_2.9.0            tidyselect_0.2.5         htmlwidgets_1.3          grid_3.6.1              
  [6] ranger_0.11.2            Rtsne_0.15               munsell_0.5.0            codetools_0.2-16         ica_1.0-2               
 [11] statmod_1.4.32           future_1.14.0            withr_2.1.2              colorspace_1.4-1         rstudioapi_0.10         
 [16] ROCR_1.0-7               gbRd_0.4-11              listenv_0.7.0            labeling_0.3             Rdpack_0.11-0           
 [21] GenomeInfoDbData_1.2.1   polyclip_1.10-0          bit64_0.9-8              farver_1.1.0             vctrs_0.2.0             
 [26] generics_0.0.2           lambda.r_1.2.3           R6_2.4.0                 GA_3.2                   ggbeeswarm_0.6.0        
 [31] rsvd_1.0.2               locfit_1.5-9.1           hdf5r_1.2.0              bitops_1.0-6             assertthat_0.2.1        
 [36] promises_1.0.1           SDMTools_1.1-221.1       scales_1.0.0             ggraph_1.0.2             beeswarm_0.2.3          
 [41] gtable_0.3.0             babelwhale_1.0.0.9000    npsurv_0.4-0             globals_0.12.4           processx_3.4.1          
 [46] tidygraph_1.1.2          rlang_0.4.0              zeallot_0.1.0            scatterplot3d_0.3-41     splines_3.6.1           
 [51] lazyeval_0.2.2           broom_0.5.2              yaml_2.2.0               reshape2_1.4.3           modelr_0.1.4            
 [56] backports_1.1.4          httpuv_1.5.1             tools_3.6.1              RColorBrewer_1.1-2       dynamicTreeCut_1.63-1   
 [61] ggridges_0.5.1           Rcpp_1.0.2               plyr_1.8.4               lambda.tools_1.0.9       zlibbioc_1.30.0         
 [66] RCurl_1.95-4.12          ps_1.3.0                 pbapply_1.4-1            viridis_0.5.1            cowplot_1.0.0           
 [71] dynparam_1.0.0.9000      zoo_1.8-6                haven_2.1.1              ggrepel_0.8.1            cluster_2.1.0           
 [76] magrittr_1.5             futile.options_1.0.1     carrier_0.1.0            lmtest_0.9-37            RANN_2.6.1              
 [81] fitdistrplus_1.0-14      patchwork_0.0.1          hms_0.5.0                lsei_1.2-0               mime_0.7.1              
 [86] xtable_1.8-4             RMTstat_0.3              readxl_1.3.1             gridExtra_2.3            testthat_2.2.1          
 [91] compiler_3.6.1           KernSmooth_2.23-15       crayon_1.3.4             rje_1.9                  R.oo_1.22.0             
 [96] htmltools_0.3.6          proxyC_0.1.5             later_0.8.0              RcppParallel_4.4.3       lubridate_1.7.4         
[101] diffusionMap_1.1-0.1     tweenr_1.0.1             formatR_1.7.1            MASS_7.3-51.4            cli_1.1.0               
[106] R.methodsS3_1.7.1        gdata_2.18.0             metap_1.1                igraph_1.2.4.1           pkgconfig_2.0.2         
[111] plotly_4.9.0             foreach_1.5.1            xml2_1.2.1               vipor_0.4.5              dqrng_0.2.1             
[116] dynutils_1.0.4.9000      XVector_0.24.0           bibtex_0.4.2             rvest_0.3.4              dyndimred_1.0.2.9000    
[121] digest_0.6.20            sctransform_0.2.0        tsne_0.1-3               cellranger_1.1.0         edgeR_3.26.6            
[126] DelayedMatrixStats_1.6.0 shiny_1.3.2              gtools_3.8.1             nlme_3.1-141             jsonlite_1.6            
[131] BiocNeighbors_1.2.0      futile.logger_1.4.3      viridisLite_0.3.0        limma_3.40.6             pillar_1.4.2            
[136] lattice_0.20-38          httr_1.4.1               survival_2.44-1.1        glue_1.3.1               remotes_2.1.0           
[141] iterators_1.0.12         png_0.1-7                bit_1.1-14               ggforce_0.2.2            stringi_1.4.3           
[146] BiocSingular_1.0.0       caTools_1.17.1.2         irlba_2.3.3              future.apply_1.3.0       ape_5.3

Thanks in advance for any suggestions!

Plotting Issue/Bug for dynplot::plot_dimred with large data sets

Dear Dynverse Team,

When using a 10X Genomics dataset with ~18k cells x ~10k genes (FULL) I came across a strange behavior with plot_dimred: To optimize conditions I am using a test dataset of 1.6k cells x ~10k genes, which gives me a sweet plot
dimred.test.PAGA_tree.ggplot2.pdf
with this code:
plot_dimred(model.tree.paga_tree, color_cells = "grouping", alpha_cells = 1, size_cells = 1) + scale_color_manual(values = color.clusters.NOmp) + theme_classic()

However, using the FULL dataset yields the default plotting
dimred.FULL.PAGA_tree.NOggplot2.pdf
despite the same code:
plot_dimred(PT.paga_tree, color_cells = "grouping", alpha_cells = 1, size_cells = 1) + scale_color_manual(values = color.clusters.NOmp) + theme_classic()

It seems the large data dimension inactivates or overrides ggplot2 coloring. Indeed, while the test dataset code will have the usual ggplot2 color scale note
Scale for 'colour' is already present. Adding another scale for 'colour', which will replace the existing scale.
which does not appear with the FULL dataset.
Strange enough the theme works, but adding i.e. geom_rug again only works for test but not FULL.

Is there a way to fix this? I would like to use the plot_dimred trajectories for presentation and publishing, truely appreciate if you can help me out here.
Many Thanks in advance,
Stephan

Error using infer_trajectory

This part goes well.

dynwrap::test_docker_installation(detailed = TRUE)
โœ” Docker is installed
โœ” Docker daemon is running
โœ” Docker is at correct version (>1.0): 1.39
โœ” Docker is in linux mode
โœ” Docker can pull images
โœ” Docker can run image
โœ” Docker can mount temporary volumes
โœ” Docker test successful -----------------------------------------------------------------

dataset <- wrap_expression(expression = example_dataset$expression,
                            counts = example_dataset$counts)
dataset <- add_prior_information(dataset,start_id = "Cell1")

answers <- dynguidelines::answer_questions(
  multiple_disconnected = TRUE, 
  expect_topology = NULL, 
  expected_topology = NULL, 
  n_cells = 1000, 
  n_features = 3, 
  memory = "30GB", 
  prior_information = "start_id", 
  docker = TRUE
)
guidelines <- dynguidelines::guidelines(answers = answers) 

This runs into error.

model <- infer_trajectory(dataset,"paga")

/usr/local/lib/python3.7/site-packages/scanpy/preprocessing/_deprecated/highly_variable_genes.py:168: RuntimeWarning: invalid value encountered in true_divide
  / disp_mad_bin[df['mean_bin'].values].values
Traceback (most recent call last):
  File "/code/run.py", line 60, in <module>
    sc.pp.recipe_zheng17(adata, n_top_genes=n_top_genes)
  File "/usr/local/lib/python3.7/site-packages/scanpy/preprocessing/_recipes.py", line 108, in recipe_zheng17
    adata.X, flavor='cell_ranger', n_top_genes=n_top_genes, log=False)
  File "/usr/local/lib/python3.7/site-packages/scanpy/preprocessing/_deprecated/highly_variable_genes.py", line 175, in filter_genes_dispersion
    disp_cut_off = dispersion_norm[n_top_genes-1]
IndexError: index 2 is out of bounds for axis 0 with size 0

I think infer_trajectory() may not be using docker. Is there an argument to specify that it should use docker and not a local python source. Also, Where is the workflow all the options in answers/guidelines actually used?

R version 3.6.0

Session Info

>> devtools::session_info()
โ”€ Session info โ”€โ”€โ”€โ”€โ”€โ”€โ”€โ”€โ”€โ”€โ”€โ”€โ”€โ”€
 setting  value                       
 version  R version 3.6.0 (2019-04-26)
 os       Ubuntu 18.04.3 LTS          
 system   x86_64, linux-gnu           
 ui       RStudio                     
 language en_GB:en                    
 collate  en_GB.UTF-8                 
 ctype    en_GB.UTF-8                 
 tz       Europe/Stockholm            
 date     2019-09-18                  

โ”€ Packages โ”€โ”€โ”€โ”€โ”€โ”€โ”€โ”€โ”€โ”€โ”€โ”€โ”€โ”€โ”€โ”€
 package       * version    date       lib source                                 
 acepack         1.4.1      2016-10-29 [1] CRAN (R 3.6.0)                         
 akima           0.6-2      2016-12-20 [1] CRAN (R 3.6.0)                         
 assertthat      0.2.1      2019-03-21 [1] CRAN (R 3.6.0)                         
 babelwhale      1.0.1      2019-09-17 [1] Github (dynverse/babelwhale@a043848)   
 backports       1.1.4      2019-04-10 [1] CRAN (R 3.6.0)                         
 base64enc       0.1-3      2015-07-28 [1] CRAN (R 3.6.0)                         
 bit             1.1-14     2018-05-29 [1] CRAN (R 3.6.0)                         
 bit64           0.9-7      2017-05-08 [1] CRAN (R 3.6.0)                         
 broom           0.5.2      2019-04-07 [1] CRAN (R 3.6.0)                         
 callr           3.3.1      2019-07-18 [1] CRAN (R 3.6.0)                         
 carrier         0.1.0      2018-10-16 [1] CRAN (R 3.6.0)                         
 cellranger      1.1.0      2016-07-27 [1] CRAN (R 3.6.0)                         
 checkmate       1.9.4      2019-07-04 [1] CRAN (R 3.6.0)                         
 cli             1.1.0      2019-03-19 [1] CRAN (R 3.6.0)                         
 cluster         2.0.8      2019-04-05 [2] CRAN (R 3.6.0)                         
 codetools       0.2-16     2018-12-24 [2] CRAN (R 3.6.0)                         
 colorspace      1.4-1      2019-03-18 [1] CRAN (R 3.6.0)                         
 cowplot         1.0.0      2019-07-11 [1] CRAN (R 3.6.0)                         
 crayon          1.3.4      2017-09-16 [1] CRAN (R 3.6.0)                         
 data.table      1.12.2     2019-04-07 [1] CRAN (R 3.6.0)                         
 desc            1.2.0      2018-05-01 [1] CRAN (R 3.6.0)                         
 devtools        2.1.0      2019-07-06 [1] CRAN (R 3.6.0)                         
 digest          0.6.20     2019-07-04 [1] CRAN (R 3.6.0)                         
 dplyr         * 0.8.3      2019-07-04 [1] CRAN (R 3.6.0)                         
 dyndimred       1.0.2.9000 2019-09-17 [1] Github (dynverse/dyndimred@5680317)    
 dynfeature    * 1.0.0.9000 2019-09-17 [1] Github (dynverse/dynfeature@93d8964)   
 dynguidelines * 1.0.0      2019-09-17 [1] Github (dynverse/dynguidelines@5181f75)
 dynmethods    * 1.0.5      2019-09-17 [1] Github (dynverse/dynmethods@53c6d1d)   
 dyno          * 0.1.1      2019-09-17 [1] Github (dynverse/dyno@d605ec9)         
 dynparam        1.0.0      2019-09-17 [1] Github (dynverse/dynparam@0979184)     
 dynplot       * 1.0.2.9000 2019-09-17 [1] Github (dynverse/dynplot@7c239aa)      
 dynutils        1.0.4.9000 2019-09-17 [1] Github (dynverse/dynutils@c160231)     
 dynwrap       * 1.1.4      2019-09-17 [1] Github (dynverse/dynwrap@019c3e6)      
 fansi           0.4.0      2018-10-05 [1] CRAN (R 3.6.0)                         
 farver          1.1.0      2018-11-20 [1] CRAN (R 3.6.0)                         
 forcats       * 0.4.0      2019-02-17 [1] CRAN (R 3.6.0)                         
 foreach         1.4.7      2019-07-27 [1] CRAN (R 3.6.0)                         
 foreign         0.8-71     2018-07-20 [2] CRAN (R 3.6.0)                         
 Formula         1.2-3      2018-05-03 [1] CRAN (R 3.6.0)                         
 fs              1.3.1      2019-05-06 [1] CRAN (R 3.6.0)                         
 GA              3.2        2019-01-10 [1] CRAN (R 3.6.0)                         
 generics        0.0.2      2018-11-29 [1] CRAN (R 3.6.0)                         
 ggforce         0.3.1      2019-08-20 [1] CRAN (R 3.6.0)                         
 ggplot2       * 3.2.1      2019-08-10 [1] CRAN (R 3.6.0)                         
 ggraph          2.0.0      2019-09-02 [1] CRAN (R 3.6.0)                         
 ggrepel         0.8.1      2019-05-07 [1] CRAN (R 3.6.0)                         
 glue            1.3.1      2019-03-12 [1] CRAN (R 3.6.0)                         
 graphlayouts    0.5.0      2019-08-20 [1] CRAN (R 3.6.0)                         
 gridExtra       2.3        2017-09-09 [1] CRAN (R 3.6.0)                         
 gtable          0.3.0      2019-03-25 [1] CRAN (R 3.6.0)                         
 haven           2.1.1      2019-07-04 [1] CRAN (R 3.6.0)                         
 hdf5r           1.2.0      2019-04-20 [1] CRAN (R 3.6.0)                         
 Hmisc           4.2-0      2019-01-26 [1] CRAN (R 3.6.0)                         
 hms             0.5.1      2019-08-23 [1] CRAN (R 3.6.0)                         
 htmlTable       1.13.1     2019-01-07 [1] CRAN (R 3.6.0)                         
 htmltools       0.3.6      2017-04-28 [1] CRAN (R 3.6.0)                         
 htmlwidgets     1.3        2018-09-30 [1] CRAN (R 3.6.0)                         
 httpuv          1.5.2      2019-09-11 [1] CRAN (R 3.6.0)                         
 httr            1.4.1      2019-08-05 [1] CRAN (R 3.6.0)                         
 igraph          1.2.4.1    2019-04-22 [1] CRAN (R 3.6.0)                         
 irlba           2.3.3      2019-02-05 [1] CRAN (R 3.6.0)                         
 iterators       1.0.12     2019-07-26 [1] CRAN (R 3.6.0)                         
 jsonlite        1.6        2018-12-07 [1] CRAN (R 3.6.0)                         
 knitr           1.24       2019-08-08 [1] CRAN (R 3.6.0)                         
 later           0.8.0      2019-02-11 [1] CRAN (R 3.6.0)                         
 lattice         0.20-38    2018-11-04 [2] CRAN (R 3.6.0)                         
 latticeExtra    0.6-28     2016-02-09 [1] CRAN (R 3.6.0)                         
 lazyeval        0.2.2      2019-03-15 [1] CRAN (R 3.6.0)                         
 lifecycle       0.1.0      2019-08-01 [1] CRAN (R 3.6.0)                         
 lmds            0.1.0      2019-09-17 [1] Github (dynverse/lmds@4c9229a)         
 lubridate       1.7.4      2018-04-11 [1] CRAN (R 3.6.0)                         
 magrittr        1.5        2014-11-22 [1] CRAN (R 3.6.0)                         
 MASS            7.3-51.4   2019-03-31 [2] CRAN (R 3.6.0)                         
 Matrix          1.2-17     2019-03-22 [2] CRAN (R 3.6.0)                         
 memoise         1.1.0      2017-04-21 [1] CRAN (R 3.6.0)                         
 mime            0.7        2019-06-11 [1] CRAN (R 3.6.0)                         
 modelr          0.1.5      2019-08-08 [1] CRAN (R 3.6.0)                         
 munsell         0.5.0      2018-06-12 [1] CRAN (R 3.6.0)                         
 nlme            3.1-139    2019-04-09 [2] CRAN (R 3.6.0)                         
 nnet            7.3-12     2016-02-02 [2] CRAN (R 3.6.0)                         
 patchwork       0.0.1      2019-09-17 [1] Github (thomasp85/patchwork@36b4918)   
 pillar          1.4.2      2019-06-29 [1] CRAN (R 3.6.0)                         
 pkgbuild        1.0.5      2019-08-26 [1] CRAN (R 3.6.0)                         
 pkgconfig       2.0.2      2018-08-16 [1] CRAN (R 3.6.0)                         
 pkgload         1.0.2      2018-10-29 [1] CRAN (R 3.6.0)                         
 plyr            1.8.4      2016-06-08 [1] CRAN (R 3.6.0)                         
 polyclip        1.10-0     2019-03-14 [1] CRAN (R 3.6.0)                         
 prettyunits     1.0.2      2015-07-13 [1] CRAN (R 3.6.0)                         
 processx        3.4.1      2019-07-18 [1] CRAN (R 3.6.0)                         
 promises        1.0.1      2018-04-13 [1] CRAN (R 3.6.0)                         
 proxyC          0.1.5      2019-07-21 [1] CRAN (R 3.6.0)                         
 ps              1.3.0      2018-12-21 [1] CRAN (R 3.6.0)                         
 purrr         * 0.3.2      2019-03-15 [1] CRAN (R 3.6.0)                         
 R6              2.4.0      2019-02-14 [1] CRAN (R 3.6.0)                         
 ranger          0.11.2     2019-03-07 [1] CRAN (R 3.6.0)                         
 RColorBrewer    1.1-2      2014-12-07 [1] CRAN (R 3.6.0)                         
 Rcpp            1.0.2      2019-07-25 [1] CRAN (R 3.6.0)                         
 RcppParallel    4.4.3      2019-05-22 [1] CRAN (R 3.6.0)                         
 readr         * 1.3.1      2018-12-21 [1] CRAN (R 3.6.0)                         
 readxl          1.3.1      2019-03-13 [1] CRAN (R 3.6.0)                         
 remotes         2.1.0      2019-06-24 [1] CRAN (R 3.6.0)                         
 reshape2        1.4.3      2017-12-11 [1] CRAN (R 3.6.0)                         
 rje             1.10.10    2019-08-28 [1] CRAN (R 3.6.0)                         
 rlang           0.4.0      2019-06-25 [1] CRAN (R 3.6.0)                         
 rpart           4.1-15     2019-04-12 [2] CRAN (R 3.6.0)                         
 rprojroot       1.3-2      2018-01-03 [1] CRAN (R 3.6.0)                         
 rstudioapi      0.10       2019-03-19 [1] CRAN (R 3.6.0)                         
 rvest           0.3.4      2019-05-15 [1] CRAN (R 3.6.0)                         
 scales          1.0.0      2018-08-09 [1] CRAN (R 3.6.0)                         
 sessioninfo     1.1.1      2018-11-05 [1] CRAN (R 3.6.0)                         
 shiny         * 1.3.2      2019-04-22 [1] CRAN (R 3.6.0)                         
 shinyjs         1.0        2018-01-08 [1] CRAN (R 3.6.0)                         
 shinyWidgets    0.4.9      2019-09-10 [1] CRAN (R 3.6.0)                         
 sp              1.3-1      2018-06-05 [1] CRAN (R 3.6.0)                         
 stringi         1.4.3      2019-03-12 [1] CRAN (R 3.6.0)                         
 stringr       * 1.4.0      2019-02-10 [1] CRAN (R 3.6.0)                         
 survival        2.44-1.1   2019-04-01 [2] CRAN (R 3.6.0)                         
 testthat        2.2.1      2019-07-25 [1] CRAN (R 3.6.0)                         
 tibble        * 2.1.3      2019-06-06 [1] CRAN (R 3.6.0)                         
 tidygraph       1.1.2      2019-02-18 [1] CRAN (R 3.6.0)                         
 tidyr         * 1.0.0      2019-09-11 [1] CRAN (R 3.6.0)                         
 tidyselect      0.2.5      2018-10-11 [1] CRAN (R 3.6.0)                         
 tidyverse     * 1.2.1      2017-11-14 [1] CRAN (R 3.6.0)                         
 tweenr          1.0.1      2018-12-14 [1] CRAN (R 3.6.0)                         
 usethis         1.5.1      2019-07-04 [1] CRAN (R 3.6.0)                         
 utf8            1.1.4      2018-05-24 [1] CRAN (R 3.6.0)                         
 vctrs           0.2.0      2019-07-05 [1] CRAN (R 3.6.0)                         
 viridis         0.5.1      2018-03-29 [1] CRAN (R 3.6.0)                         
 viridisLite     0.3.0      2018-02-01 [1] CRAN (R 3.6.0)                         
 withr           2.1.2      2018-03-15 [1] CRAN (R 3.6.0)                         
 xfun            0.9        2019-08-21 [1] CRAN (R 3.6.0)                         
 xml2            1.2.2      2019-08-09 [1] CRAN (R 3.6.0)                         
 xtable          1.8-4      2019-04-21 [1] CRAN (R 3.6.0)                         
 yaml            2.2.0      2018-07-25 [1] CRAN (R 3.6.0)                         
 zeallot         0.1.0      2018-01-28 [1] CRAN (R 3.6.0) 

"Error during trajectory inference" => out of memory

Hello,

I am currently getting this error:

model <- infer_trajectory(dataset, methods_selected[1])
Error traceback:
1:
Attaching package: โ€˜dplyrโ€™

The following objects are masked from โ€˜package:statsโ€™:

filter, lag

The following objects are masked from โ€˜package:baseโ€™:

intersect, setdiff, setequal, union

Attaching package: โ€˜purrrโ€™

The following object is masked from โ€˜package:jsonliteโ€™:

flatten

Loading required package: dynutils
Killed

Error: Error during trajectory inference

Any help is appreciated!

Note: I am using methods_selected[1] because when I use first(methods_selected) I get:
Error in (function (classes, fdef, mtable) :
unable to find an inherited method for function โ€˜firstโ€™ for signature โ€˜"character"โ€™

I am not sure if this is related.

Docker toolbox install windows with virtualbox: mounting temporary volumes

Hey there,

So I am very novel (new) with rstudio and github, but I have been studying TI and have become very interested in the area, leading me to the Dynverse. However, I am having trouble with the installation, not being able to enable the final docker install check (i.e. mounting temp. Volumes: using docker toolbox with virtualbox on windows 10). Tried everything; any help would be deeply appreciated. Thank you.

Error in eval_tidy(enquo(var), var_env) : object 'method_id' not found

Hi,

I met the error message when I run the tutorial code. "Error in eval_tidy(enquo(var), var_env) : object 'method_id' not found". Could you try to help me out? Thanks a lot.

the details show in below.

guidelines <- guidelines_shiny(dataset)
Loading required package: shiny

Listening on http://0.0.0.0:3209
Code to reproduce the guidelines, copy it over to your script!

Reproduces the guidelines as created in the shiny app

answers <- dynguidelines::answer_questions(
multiple_disconnected = NULL,
expect_topology = NULL,
expected_topology = NULL,
n_cells = 392,
n_features = 2000,
prior_information = NULL
)
guidelines <- dynguidelines::guidelines(answers = answers)

methods <- guidelines$methods %>% filter(selected) %>% pull(method_id) %>% first()
Error in eval_tidy(enquo(var), var_env) : object 'method_id' not found

error : $ operator is invalid for atomic vectors

I've been trying to run Dyno so as to proceed to lineage inference from a public dataset.

I've done the preprocessing to reshape de data to the required input format (counts and expression matrices with cells as rows and genes as columns). I've assigned cells IDs and genes as rownames and colnames.

When trying to run PAGA, the following error pops out: $ operator is invalid for atomic vectors

Here's my piece of code so far. Any insights on why this is happening? Thanks in advance!

#Create a dyno dataset with wrap_expression

library(tidyverse)
library(dyno)

dataset <- wrap_expression(counts = counts, expression = exp)

#Run dynguidelines

guidelines <- guidelines_shiny(dataset)

answers <- dynguidelines::answer_questions(
multiple_disconnected = NULL,
expect_topology = NULL,
expected_topology = NULL,
n_cells = 19743,
n_features = 20921,
memory = "8GB",
prior_information = c("start_id", "end_id", "end_n", "start_n", "groups_n", "features_id"),
docker = TRUE
)

#Adding prior information to the dataset (PAGA needs it to run)

dataset <- add_prior_information(dataset, start_id = "arc1_CAAATCACAAGC")

#Reproduces the guidelines as created in the shiny app

guidelines <- dynguidelines::guidelines(answers = answers)

methods_selected <- guidelines$methods_selected

#Run the selected method

model <- infer_trajectory(dataset, first(methods_selected))

As an output from the last line, I get:

model <- infer_trajectory(dataset, first(methods_selected))
Error traceback:
1: Killed

Error: $ operator is invalid for atomic vectors

Error: Error during trajectory inference any(duplicated(c(cell_ids, milestone_ids))) isn't false.

I'm trying to run projected_paga on my dataset and get this error:

Error: Error during trajectory inference
any(duplicated(c(cell_ids, milestone_ids))) isn't false.

I can't figure out where milestone_ids comes from.

wrapped_data <- wrap_expression( counts = param_ion_counts, expression = param_asinh_expression )

any(duplicated(wrapped_data$cell_ids))

evaluates to FALSE

I can't share the dataset publicly, but if you can't reproduce this error I can probably email the dataset to you.

Best,
Geoff

Correct usage of prior information with TI methods

I have inferred a trajectory with dyno using slingshot between 3 groups of cells, however prior knowledge suggests that the trajectory should go from groups A -> B -> C where as slingshot has returned A -> B and A -> C (i.e a bifurcation). According to the dyno guidelines, slingshot is able to make use of prior information (e.g. start/end states) however, whatever priors I provide slingshot just gives me the same trajectory. Is this expected behaviour or am I doing something wrong?

Below is an example of the commands I would use to run slingshot with a hypothetical dataset containing 3 cells per group:

dat <- add_prior_information(
  dataset = dat
  start_id = c("A1, "A2", "A3"),
  end_id = c("C1", "C2", "C3"),
  groups_id = data.frame(cell_id = c("A1, "A2", "A3", "B1", "B2", "B3", "C1", "C2", "C3"),
                         group_id = c("A", "A", "A", "B", "B", "B", "C", "C", "C"),
  groups_network = data.frame(from = c("A", "B"), to = c("B", "C"))
  start_n = 1,
  end_n = 1
)

mod <- infer_trajectory(dat, "slingshot", give_priors = c("start_id", "end_id", "groups_id", "groups_network", "start_n", "end_n"), seed = 1701)

I am using the following package versions:

  • dynwrap_1.1.2
  • dynfeature_1.0.0
  • dynplot_1.0.1
  • dynmethods_1.0.1
  • dyno_0.1.1
  • dynguidelines_1.0.0

Error: Column `is_start` can't be converted from integer to logical

Hi,

I keep getting an error when trying to plot the trajectory on a dim reduction using plot_dimred() function.

below is an example of analysis flow:

raw.data.spleen <- read.table("output/tables/raw_dataframe_Spleen_only.txt")
normed.data.spleen <- read.table("output/tables/filtered_dataframe_Spleen_only.txt")

Dyno.spleen.data <- wrap_expression(
counts = raw.data.spleen,
expression = normed.data.spleen
)

answers.spleen <- dynguidelines::answer_questions(
multiple_disconnected = NULL,
expect_topology = NULL,
expected_topology = NULL,
n_cells = 3447,
n_features = 12236,
memory = "8GB",
prior_information = NULL,
docker = TRUE
)

slingshot.spleen <- infer_trajectory(Dyno.spleen.data,
method = "dynverse/ti_slingshot:v0.9.9.01",
parameters = NULL,
give_priors = NULL,
seed = set.seed(42),
verbose = TRUE)

plot_dimred(slingshot.spleen)

Coloring by milestone
Using milestone_percentages from trajectory
Error: Column is_start can't be converted from integer to logical

traceback()
13: stop(list(message = "Column is_start can't be converted from integer to logical",
call = NULL, cppstack = NULL))
12: bind_rows_(x, .id)
11: bind_rows(divergence_regions, extra_divergences)
10: calculate_geodesic_distances_(cell_ids = trajectory$cell_ids,
milestone_ids = trajectory$milestone_ids, milestone_network = trajectory$milestone_network,
milestone_percentages = trajectory$milestone_percentages,
divergence_regions = trajectory$divergence_regions, waypoint_cells = waypoint_cells,
waypoint_milestone_percentages = waypoint_milestone_percentages)
9: calculate_geodesic_distances(trajectory, waypoint_milestone_percentages = waypoint_milestone_percentages)
8: dynwrap::select_waypoints(trajectory)
7: is.data.frame(x)
6: colnames(waypoints$geodesic_distances)
5: eval_bare(get_expr(quo), get_env(quo))
4: quasi_label(enquo(expected))
3: testthat::expect_setequal(cell_positions$cell_id, colnames(waypoints$geodesic_distances))
2: project_waypoints(trajectory = trajectory, cell_positions = cell_positions,
waypoints = waypoints, trajectory_projection_sd = trajectory_projection_sd,
color_trajectory = color_trajectory)
1: plot_dimred(slingshot.liver)

IndexError running paga or paga_tree

Hi there,

I am trying to use your package on flow cytometry data using an examplary dataset of ~30.000 cells and 5 dimensions with a pre-set start_id. The fluorescence data is wrapped into counts and expression of a dataset.

Unfortunately when running

model_paga_tree <- infer_trajectory(dataset, ti_paga(), verbose = TRUE)

I get the following error

Executing 'paga' on '20190627_145358__data_wrapper__8U7MT2djO2'
With parameters: list(n_neighbors = 15L, n_comps = 50L, n_dcs = 15L, resolution = 1L,     embedding_type = "fa", connectivity_cutoff = 0.05),
inputs: counts, and
priors : start_id
Input saved to C:\Users\...\AppData\Local\Temp\RtmpagjEd1\file28381f963609/ti
Running method using babelwhale
Running "C:\PROGRA~1\DOCKER~1\docker.exe" run -e "TMPDIR=/tmp2" --workdir /ti/workspace -v "/c/Users/.../AppData/Local/Temp/RtmpagjEd1/file28381f963609/ti:/ti" -v \
  "/c/Users/.../AppData/Local/Temp/RtmpagjEd1/file2838ad45788/tmp:/tmp2" "dynverse/ti_paga:v0.9.9.04" --dataset /ti/input.h5 --output /ti/output.h5 --use_priors all
/usr/local/lib/python3.7/site-packages/scanpy/preprocessing/_deprecated/highly_variable_genes.py:168: RuntimeWarning: invalid value encountered in true_divide
  / disp_mad_bin[df['mean_bin'].values].values
Traceback (most recent call last):
  File "/code/run.py", line 55, in <module>
    sc.pp.recipe_zheng17(adata, n_top_genes=n_top_genes)
  File "/usr/local/lib/python3.7/site-packages/scanpy/preprocessing/_recipes.py", line 108, in recipe_zheng17
    adata.X, flavor='cell_ranger', n_top_genes=n_top_genes, log=False)
  File "/usr/local/lib/python3.7/site-packages/scanpy/preprocessing/_deprecated/highly_variable_genes.py", line 175, in filter_genes_dispersion
    disp_cut_off = dispersion_norm[n_top_genes-1]
IndexError: index 4 is out of bounds for axis 0 with size 0
Error: Error during trajectory inference 
/usr/local/lib/python3.7/site-packages/scanpy/preprocessing/_deprecated/highly_variable_genes.py:168: RuntimeWarning: invalid value encountered in true_divide
  / disp_mad_bin[df['mean_bin'].values].values
Traceback (most recent call last):
  File "/code/run.py", line 55, in <module>
    sc.pp.recipe_zheng17(adata, n_top_genes=n_top_genes)
  File "/usr/local/lib/python3.7/site-packages/scanpy/preprocessing/_recipes.py", line 108, in recipe_zheng17
    adata.X, flavor='cell_ranger', n_top_genes=n_top_genes, log=False)
  File "/usr/local/lib/python3.7/site-packages/scanpy/preprocessing/_deprecated/highly_variable_genes.py", line 175, in filter_genes_dispersion
    disp_cut_off = dispersion_norm[n_top_genes-1]
IndexError: index 4 is out of bounds for axis 0 with size 0

Exactly the same is happening with paga_tree. Can you tell me what's going wrong?

Thanks and best,
Laura

Question: SCE to H5AD

Hi,
this is just a question, not really an issue -- I was just wondering whether you could point me to a function that would enable me to convert my SingleCellExperiment objects into scanpy-ready H5AD objects?

Cheers,

Friederike

When and where to add biological data?

Given that several of the methods available through dyno require root cells, or can make use of time point data, it is unclear at what step you add these to the model. Should these always be added after running infer_trajectories? Is there a key to know which methods require what data a priori?

Issue with singularity when running infer_trajectory

Hi,

thanks for developing this amazing package and for the incredible amount of work you put into benchmarking all of these tools, a great effort!

I was trying to run the tutorial to test dyno's functionalities. I am using dyno_0.9.9 and tidyverse_1.2.1.

I first tried on my local machine:

R version 3.5.2 (2018-12-20)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu 18.04.2 LTS

I installed singularity and the checks from dynwrap::test_singularity_installation(detailed = TRUE) are all positive. I'm using singularity version 3.0.3.When I try to run any method using infer_trajectory I get the following error:

methods_selected <- c("slingshot" "paga_tree" "scorpius" "angle" )

model <- infer_trajectory(dataset, methods_selected[1], verbose = TRUE, debug = TRUE)

hi hi Error: FATAL: failed to retrieved path for /home/florian/Dyno_trajectories/.singularity/dynverse/ti_slingshot_tag-v0.9.9.simg: lstat /home/florian/Dyno_trajectories/.singularity/dynverse/ti_slingshot_tag-v0.9.9.simg: no such file or directory ERROR : Child exit with status 255

I have tried the same procedure on an HPC cluster that is runnig R 3.5.0 where singularity 3.1 is installed as a module that we can load. I get the exact same error that the simg file is not found.

Could you please help me troubleshoot this issue?

Thanks,

Florian

Running slingshot on low dimensional data

Hi there,
I am just trying to use your package on flow cytometry data. I read in the User Guide that "Currently, alternative input data such as ATAC-Seq or cytometry data are not yet supported, although it is possible to simply include this data as expression and counts."

I am trying an examplary dataset of ~30.000 cells and 5 dimensions with a pre-set start_id, because I know how the cells develop concerning these 5 dimensions. I wrapped this data into counts and expression of a dataset.

Unfortunately when running

model_slingshot <- infer_trajectory(dataset, ti_slingshot(), verbose = TRUE)

I get the following error

Executing 'slingshot' on '20190626_154646__data_wrapper__e6QaaLwGP3'
With parameters: list(shrink = 1L, reweight = TRUE, reassign = TRUE, thresh = 0.001,     maxit = 10L, stretch = 2L, smoother = "smooth.spline", shrink.method = "cosine"),
inputs: expression, and
priors : 
Input saved to C:\Users\...\AppData\Local\Temp\RtmpKGtDxe\file29f863416e7c/ti
Running method using babelwhale
Running "C:\PROGRA~1\DOCKER~1\docker.exe" run -e "TMPDIR=/tmp2" --workdir /ti/workspace -v "/c/Users/.../AppData/Local/Temp/RtmpKGtDxe/file29f863416e7c/ti:/ti" -v \
  "/c/Users/.../AppData/Local/Temp/RtmpKGtDxe/file29f843c423f4/tmp:/tmp2" "dynverse/ti_slingshot:v0.9.9.01" --dataset /ti/input.h5 --output /ti/output.h5
Loading required package: princurve
Error in (function (A, nv = 5, nu = nv, maxit = 1000, work = nv + 7, reorth = TRUE,  : 
  max(nu, nv) must be strictly less than min(nrow(A), ncol(A))
Calls: <Anonymous> -> do.call -> <Anonymous>
Execution halted
Error: Error during trajectory inference 
Loading required package: princurve
Error in (function (A, nv = 5, nu = nv, maxit = 1000, work = nv + 7, reorth = TRUE,  : 
  max(nu, nv) must be strictly less than min(nrow(A), ncol(A))
Calls: <Anonymous> -> do.call -> <Anonymous>
Execution halted

I think that the following line in ti_slingshot causes the problem since n=20 is fixed here and my dataset does have less than 20 dimensions :

pca <- irlba::prcomp_irlba(expression, n = 20)

Or am I missing something?

Thanks in advance,
Laura

dynwrap: cell_ids swapped with feature_ids

I created the dynwrap using the wrap_expression and the counts and expression matrix.

object_counts<- as(as.matrix(object@assays$RNA@counts), 'sparseMatrix')
object_expression<- as(as.matrix(object@assays$RNA@data), 'sparseMatrix')
object_dyn<- wrap_expression(counts = object_counts, 
                                              expression = object_expression)

However, I found that the cell_ids and feature_ids were swapped. All of the cell_ids were in the feature_ids column and vice-versa. I swapped them manually. However, when I created the model, everything went back again to how they were, all the labels were mixed up. Not sure if the model is now accurate or not, and how to fix this.
image

Error from running paga method

Hi dyno community,

I was trying to run slingshot and paga methods on my data. Since paga method required a prior, I set "start_id" based on some assumptions. Using the dataset, I could get slingshot trajectory without any issue, but for the paga method I got following issue and not so sure how to interpret this error report. Has anyone got the same error before and what would be a way to resolve this? Thanks!!

R[write to console]: Error in parse_args(parser, args = args) :

Error in getopt(spec = spec, opt = args) :
long flag "use_priors" is invalid
Calls: -> -> -> parse_args

Traceback (most recent call last):
File "/code/run.py", line 22, in
task = dynclipy.main()
File "/usr/local/lib/python3.7/site-packages/dynclipy/read.py", line 123, in main
""")
File "/usr/local/lib/python3.7/site-packages/rpy2/robjects/init.py", line 389, in call
res = self.eval(p)
File "/usr/local/lib/python3.7/site-packages/rpy2/robjects/functions.py", line 192, in call
.call(*args, **kwargs))
File "/usr/local/lib/python3.7/site-packages/rpy2/robjects/functions.py", line 121, in call
res = super(Function, self).call(*new_args, **new_kwargs)
File "/usr/local/lib/python3.7/site-packages/rpy2/rinterface_lib/conversion.py", line 28, in _
cdata = function(*args, **kwargs)
File "/usr/local/lib/python3.7/site-packages/rpy2/rinterface.py", line 773, in call
raise embedded.RRuntimeError(_rinterface._geterrmessage())
rpy2.rinterface_lib.embedded.RRuntimeError: Error in parse_args(parser, args = args) :
Error in getopt(spec = spec, opt = args) :
long flag "use_priors" is invalid
Calls: -> -> -> parse_args

Error: Error during trajectory inference
R[write to console]: Error in parse_args(parser, args = args) :
Error in getopt(spec = spec, opt = args) :
long flag "use_priors" is invalid
Calls: -> -> -> parse_args

Traceback (most recent call last):
File "/code/run.py", line 22, in
task = dynclipy.main()
File "/usr/local/lib/python3.7/site-packages/dynclipy/read.py", line 123, in main
""")
File "/usr/local/lib/python3.7/site-packages/rpy2/robjects/init.py", line 389, in call
res = self.eval(p)
File "/usr/local/lib/python3.7/site-packages/rpy2/robjects/functions.py", line 192, in call
.call(*args, **kwargs))
File "/usr/local/lib/python3.7/site-packages/rpy2/robjects/functions.py", line 121, in call
res = super(Function, self).call(*new_args, **new_kwargs)
File "/usr/local/lib/python3.7/site-packages/rpy2/rinterface_lib/conversion.py", line 28, in _
cdata = function(*args, **kw

debugging "Error: Error during trajectory inference"

I'm having trouble running dynwrap::infer_trajectory() on my data.
My installation seems to be okay, the methods work fine on the toy data.

However, on my data (and verbose=TRUE), I get a lot of package loading messages and:

1. dynwrap::infer_trajectory(dataset, method = ti_monocle_ddrtree(reduction_method = "DDRTree", 
 .     max_components = 5, norm_method = "log", filter_features_mean_expression = 0.1), 
 .     verbose = TRUE)
2. stop("Error during trajectory inference \n", crayon::bold(error), 
 .     call. = FALSE)

As my data was running as expected through Monocle in the past, I followed
https://github.com/dynverse/ti_monocle_ddrtree/blob/master/run.R
line by line, and everything runs up to line 79.

Do you have any tips on how I can debug this? (I tried several other ti_ methods as well, with similar results)

Thanks for this really cool package!

authentication error running singularity?

this worked for a me a couple weeks ago on our cluster but not any more:

$ singularity exec 'docker://dynverse/ti_slingshot:v1.0.3' echo hi
FATAL:   Unable to handle docker://dynverse/ti_slingshot:v1.0.3 uri: failed to get SHA of 
docker://dynverse/ti_slingshot:v1.0.3: error getting username and password: error reading 
JSON file "/run/user/68333/containers/auth.json": error unmarshaling JSON at "/run/user
/68333/containers/auth.json": unexpected end of JSON input

Error during trajectory inference: unused argument (env = processx_env)

I am following the dyno demo (https://github.com/dynverse/dyno#installation) and run into error at the following step. Prior to this step, because the "guidelines <- guidelines_shiny(dataset)" doesn't work, so I manually set the value of "methods_selected" to "scorpius".

model <- infer_trajectory(dataset, first(methods_selected))
Error traceback:
2: unused argument (env = processx_env)
1: processx::run
processx_command
processx_args
processx_env
verbose
verbose
TRUE
FALSE
Error: Error during trajectory inference
unused argument (env = processx_env)

trying to get this working on the cluster

I am submitting this job on our cluster using slurm with:

models <- dynwrap::infer_trajectories(
  dataset, 
  method = list(ti_slingshot(), ti_mst(), ti_slice(), ti_pcreode())
)

and this is returned to me. is it normal for the protocol string for singularity to say docker?

STDOUT
Running singularity exec 'docker://dynverse/ti_slingshot:v1.0.3' echo hi

STDERR
Error in processx::run("singularity", c("exec", paste0("docker://", container_id),  : 
  System command error
Calls: <Anonymous> ... create_ti_method_container -> <Anonymous> -> <Anonymous>
Execution halted
$ singularity --version
singularity version 3.3.0-1.el7

Missing required dependency [projected_paga]

Hi dyno team,

After managing to run dyno on our HPC, I am running into a weird dependency issue when trying to run "projected_paga" on the test dataset. The same code works fine when using a different method, for example "angle". My colleague is able to run the exact same analysis on the same HPC, leading me to believe it must be related to some python packages I have previously installed and he has not. However, I don't fully understand why dyno is trying to use my local python instead of the python inside the singularity container.

Please find the code and error below:

library(dyno)
library(tidyverse)

data("fibroblast_reprogramming_treutlein")

dataset <- wrap_expression(
counts = fibroblast_reprogramming_treutlein$counts,
expression = fibroblast_reprogramming_treutlein$expression
)

Executing 'projected_paga' on '20190409_063127__data_wrapper__ygGdphzfxh'
With parameters: list(n_neighbors = 15L, n_comps = 50L, n_dcs = 15L, resolution = 1L,     embedding_type = "fa", connectivity_cutoff = 0.05),
inputs: counts, and
priors : 
Input saved to /tmp/RtmpuyJ9Xp/file2a3f56e8a7d33/ti
Running /opt/software/singularity-3.1/bin/singularity run --pwd /ti/workspace \
  -B \
  '/tmp/RtmpuyJ9Xp/file2a3f56e8a7d33/ti:/ti,/tmp/RtmpuyJ9Xp/file2a3f514c36879/tmp:/tmp2' \
  'docker://dynverse/ti_projected_paga:v0.9.9.01' --dataset /ti/input.h5 \
  --output /ti/output.h5
WARNING: Could not set container working directory /ti/workspace: chdir /ti/workspace: no such file or directory
Traceback (most recent call last):
  File "/code/run.py", line 3, in <module>
    import dynclipy
  File "/usr/local/lib/python3.7/site-packages/dynclipy/__init__.py", line 3, in <module>
    from rpy2.robjects import pandas2ri
  File "/usr/local/lib/python3.7/site-packages/rpy2/robjects/pandas2ri.py", line 13, in <module>
    from pandas.core.frame import DataFrame as PandasDataFrame
  File "/home/wueflo00/.local/lib/python3.7/site-packages/pandas/__init__.py", line 19, in <module>
    "Missing required dependencies {0}".format(missing_dependencies))
ImportError: Missing required dependencies ['numpy']
Error: Error during trajectory inference 
WARNING: Could not set container working directory /ti/workspace: chdir /ti/workspace: no such file or directory
Traceback (most recent call last):
  File "/code/run.py", line 3, in <module>
    import dynclipy
  File "/usr/local/lib/python3.7/site-packages/dynclipy/__init__.py", line 3, in <module>
    from rpy2.robjects import pandas2ri
  File "/usr/local/lib/python3.7/site-packages/rpy2/robjects/pandas2ri.py", line 13, in <module>
    from pandas.core.frame import DataFrame as PandasDataFrame
  File "/home/wueflo00/.local/lib/python3.7/site-packages/pandas/__init__.py", line 19, in <module>
    "Missing required dependencies {0}".format(missing_dependencies))
ImportError: Missing required dependencies ['numpy']

Error: "cell_id" %in% colnames(end_state_probabilities) isn't true.

Hi,

When I use the sample data and try the grandprix method, I meet this error:

library(dyno)
library(tidyverse)
data("fibroblast_reprogramming_treutlein")
dataset <- wrap_expression(
counts = fibroblast_reprogramming_treutlein$counts,
expression = fibroblast_reprogramming_treutlein$expression
) %>%
add_prior_information(
end_n = 1
)
model <- infer_trajectory(dataset = dataset, method = ti_grandprix(), verbose = TRUE)
Executing 'grandprix' on '20181104_034735__data_wrapper__dwq1Aondcn'
With parameters: list()
And inputs: expression, end_n
Input saved to /tmp/Rtmp1ghIhj/file10af56aa4118/ti/input:
end_n.json
expression.csv
params.json
Running /usr/bin/docker run -e 'TMPDIR=/tmp2' --workdir /ti/workspace -v
'/tmp/Rtmp1ghIhj/file10af56aa4118/ti:/ti' -v
'/tmp/Rtmp1ghIhj/file10af8de69b4/tmp:/tmp2' dynverse/ti_grandprix
2018-11-04 03:47:44.038317: I tensorflow/core/platform/cpu_feature_guard.cc:141] Your CPU supports instructions that this TensorFlow binary was not compiled to use: AVX2 FMA
output saved in /tmp/Rtmp1ghIhj/file10af56aa4118/ti/output:
cell_ids.csv
end_state_probabilities.csv
pseudotime.csv
timings.json
Error traceback:
2: "cell_id" %in% colnames(end_state_probabilities) isn't true.
1:
Error: Error during trajectory inference
"cell_id" %in% colnames(end_state_probabilities) isn't true.

And I have another two questions:
1, Could you please explain me more about the parameter end_n?
2, I want to add timecourse information to the grandprix method, but I don't know what the timecourse table look like. Could you please generate a sample timecourse table? Thanks.

Yifang

error running slingshot

Hi,

I've tried my hands at slingshot via dynwrap and that didn't pan out.

> library(dynwrap)
> assayNames(sce)
[1] "counts"    "logcounts"
> dataset <- wrap_expression(counts = counts(sce),
                             +   expression =logcounts(sce)
                             + )
> model <- infer_trajectory(dataset,"slingshot")
Loading required namespace: dynmethods
Using default tag: latest
latest: Pulling from dynverse/ti_slingshot
bc9ab73e5b14: Pulling fs layer
846568129ac2: Pulling fs layer
85acc287e507: Pulling fs layer
1815dfbed2cb: Pulling fs layer
f340461deaa1: Pulling fs layer
1f59c4d5538f: Pulling fs layer
c902428c9782: Pulling fs layer
dfec0a49fb3b: Pulling fs layer
138cd0e8b92c: Pulling fs layer
94d2b05985da: Pulling fs layer
af230cdeb663: Pulling fs layer
b9ae738f1559: Pulling fs layer
96270503d443: Pulling fs layer
5ce5e3b13e39: Pulling fs layer
c3d16d17d155: Pulling fs layer
95406573f667: Pulling fs layer
e642cfd8a979: Pulling fs layer
1815dfbed2cb: Waiting
f340461deaa1: Waiting
1f59c4d5538f: Waiting
c902428c9782: Waiting
dfec0a49fb3b: Waiting
138cd0e8b92c: Waiting
94d2b05985da: Waiting
af230cdeb663: Waiting
b9ae738f1559: Waiting
96270503d443: Waiting
5ce5e3b13e39: Waiting
c3d16d17d155: Waiting
95406573f667: Waiting
e642cfd8a979: Waiting
bc9ab73e5b14: Retrying in 5 seconds
bc9ab73e5b14: Retrying in 4 seconds
bc9ab73e5b14: Retrying in 3 seconds
bc9ab73e5b14: Retrying in 2 seconds
bc9ab73e5b14: Retrying in 1 second
bc9ab73e5b14: Verifying Checksum
bc9ab73e5b14: Download complete
1815dfbed2cb: Verifying Checksum
1815dfbed2cb: Download complete
f340461deaa1: Verifying Checksum
f340461deaa1: Download complete
1f59c4d5538f: Verifying Checksum
1f59c4d5538f: Download complete
bc9ab73e5b14: Pull complete
c902428c9782: Verifying Checksum
c902428c9782: Download complete
85acc287e507: Verifying Checksum
85acc287e507: Download complete
138cd0e8b92c: Retrying in 5 seconds
138cd0e8b92c: Retrying in 4 seconds
138cd0e8b92c: Retrying in 3 seconds
138cd0e8b92c: Retrying in 2 seconds
138cd0e8b92c: Retrying in 1 second
138cd0e8b92c: Verifying Checksum
138cd0e8b92c: Download complete
94d2b05985da: Download complete
846568129ac2: Download complete
846568129ac2: Pull complete
b9ae738f1559: Verifying Checksum
b9ae738f1559: Download complete
96270503d443: Verifying Checksum
96270503d443: Download complete
85acc287e507: Pull complete
1815dfbed2cb: Pull complete
f340461deaa1: Pull complete
1f59c4d5538f: Pull complete
c902428c9782: Pull complete
5ce5e3b13e39: Verifying Checksum
5ce5e3b13e39: Download complete
c3d16d17d155: Verifying Checksum
c3d16d17d155: Download complete
af230cdeb663: Retrying in 4 seconds
af230cdeb663: Retrying in 3 seconds
af230cdeb663: Retrying in 2 seconds
af230cdeb663: Retrying in 1 second
af230cdeb663: Verifying Checksum
af230cdeb663: Download complete
e642cfd8a979: Verifying Checksum
e642cfd8a979: Download complete
95406573f667: Verifying Checksum
95406573f667: Download complete
dfec0a49fb3b: Verifying Checksum
dfec0a49fb3b: Download complete
dfec0a49fb3b: Pull complete
138cd0e8b92c: Pull complete
94d2b05985da: Pull complete
af230cdeb663: Pull complete
b9ae738f1559: Pull complete
96270503d443: Pull complete
5ce5e3b13e39: Pull complete
c3d16d17d155: Pull complete
95406573f667: Pull complete
e642cfd8a979: Pull complete
Digest: sha256:9408fbc09ce683f658344fcda7c2590928f4ff577098123d883c7039f324e5c5
Status: Downloaded newer image for dynverse/ti_slingshot:latest
Error traceback:
  1: 
  Attaching package: โ€˜dplyrโ€™

The following objects are masked from โ€˜package:statsโ€™:
  
  filter, lag

The following objects are masked from โ€˜package:baseโ€™:
  
  intersect, setdiff, setequal, union


Attaching package: โ€˜purrrโ€™

The following object is masked from โ€˜package:jsonliteโ€™:
  
  flatten

Warning messages:
  1: In rgl.init(initValue, onlyNULL) : RGL: unable to open X11 display
2: 'rgl_init' failed, running with rgl.useNULL = TRUE 
Killed

Error: $ operator is invalid for atomic vectors
In addition: Warning message:
  In babelwhale::get_label(container_id, "version") : 
  
  Error: No such image: dynverse/ti_slingshot

Error in "model <- infer_trajectory(dataset, methods_selected, verbose = TRUE)"

Hi,

I have a problem to run the code

model <- infer_trajectory(dataset, methods_selected, verbose = TRUE)

In June I tried your package with some single cell data and followed your tutorial for that and it worked totally fine. A few days ago I tried to run exactly the same code with the same dataset, but it is not working anymore. So as there was an update for docker, I tried to set it back to the version I used before, but that didn't help.
I run the following code:

library(Seurat)
library(dyno)
library(Matrix)
library(dplyr)

seu.data <- Read10X(data.dir = path)
seu<-CreateSeuratObject(seu.data, min.cells = 2,min.features = 200)
seuNorm <- NormalizeData(object = seu, normalization.method = "LogNormalize", scale.factor = 10000)

dataset <-wrap_expression(id         = "MyData",
                          expression = t(seuNorm@assays[[1]]@data),
                          counts     = t(seuNorm@assays[[1]]@counts) )

guidelines <- guidelines_shiny(dataset)
methods_selected <- guidelines$methods_selected[1] 

dataset <- add_prior_information(
  dataset,
  start_id = "AAACCTGAGAGGTTAT" 
)

model <- infer_trajectory(dataset, methods_selected,  verbose = TRUE)

I chose the method PAGA for a quick test, but R computes for hours. The output during that time is:

v0.9.9.04: Pulling from dynverse/ti_paga_tree
ab1fc7e4bf91: Pulling fs layer
66575050199c: Pulling fs layer
4d6ccd149fc3: Pulling fs layer
d588c72fe3f7: Pulling fs layer
c4efdd5f451e: Pulling fs layer
834a41c4efa6: Pulling fs layer
82fa28fdbeb3: Pulling fs layer
d9dcdde68429: Pulling fs layer
0512db24e6b5: Pulling fs layer
b8a1c4f7f22e: Pulling fs layer
035dd037a0ee: Pulling fs layer
fea30c11d7f1: Pulling fs layer
4fe9981d22c9: Pulling fs layer
696bc3423178: Pulling fs layer
0a65a0f5dcd1: Pulling fs layer
d0bb45d42c8e: Pulling fs layer
1a9616399b58: Pulling fs layer
32cd32510543: Pulling fs layer
6508acd2eecd: Pulling fs layer
479ce1a1b674: Pulling fs layer
6b6f799b42df: Pulling fs layer
868e7e7207d7: Pulling fs layer
aa9bdfcbcdf9: Pulling fs layer
fea30c11d7f1: Waiting
4fe9981d22c9: Waiting
696bc3423178: Waiting
0a65a0f5dcd1: Waiting
d0bb45d42c8e: Waiting
1a9616399b58: Waiting
32cd32510543: Waiting
6508acd2eecd: Waiting
479ce1a1b674: Waiting
6b6f799b42df: Waiting
868e7e7207d7: Waiting
aa9bdfcbcdf9: Waiting
82fa28fdbeb3: Waiting
d588c72fe3f7: Waiting
c4efdd5f451e: Waiting
834a41c4efa6: Waiting
d9dcdde68429: Waiting
0512db24e6b5: Waiting
b8a1c4f7f22e: Waiting
035dd037a0ee: Waiting
ab1fc7e4bf91: Verifying Checksum
ab1fc7e4bf91: Download complete
d588c72fe3f7: Verifying Checksum
d588c72fe3f7: Download complete
c4efdd5f451e: Verifying Checksum
c4efdd5f451e: Download complete
ab1fc7e4bf91: Pull complete
834a41c4efa6: Download complete
82fa28fdbeb3: Verifying Checksum
82fa28fdbeb3: Download complete
4d6ccd149fc3: Verifying Checksum
4d6ccd149fc3: Download complete
0512db24e6b5: Download complete
b8a1c4f7f22e: Verifying Checksum
b8a1c4f7f22e: Download complete
66575050199c: Verifying Checksum
66575050199c: Download complete
fea30c11d7f1: Verifying Checksum
fea30c11d7f1: Download complete
4fe9981d22c9: Verifying Checksum
4fe9981d22c9: Download complete
696bc3423178: Verifying Checksum
696bc3423178: Download complete
0a65a0f5dcd1: Download complete
66575050199c: Pull complete
035dd037a0ee: Verifying Checksum
035dd037a0ee: Download complete
d9dcdde68429: Verifying Checksum
d9dcdde68429: Download complete
1a9616399b58: Verifying Checksum
1a9616399b58: Download complete
d0bb45d42c8e: Verifying Checksum
d0bb45d42c8e: Download complete
32cd32510543: Verifying Checksum
32cd32510543: Download complete
6508acd2eecd: Verifying Checksum
6508acd2eecd: Download complete
6b6f799b42df: Verifying Checksum
6b6f799b42df: Download complete
4d6ccd149fc3: Pull complete
d588c72fe3f7: Pull complete
868e7e7207d7: Verifying Checksum
868e7e7207d7: Download complete
c4efdd5f451e: Pull complete
aa9bdfcbcdf9: Verifying Checksum
aa9bdfcbcdf9: Download complete
834a41c4efa6: Pull complete
82fa28fdbeb3: Pull complete
479ce1a1b674: Verifying Checksum
479ce1a1b674: Download complete
d9dcdde68429: Pull complete
0512db24e6b5: Pull complete
b8a1c4f7f22e: Pull complete
035dd037a0ee: Pull complete
fea30c11d7f1: Pull complete
4fe9981d22c9: Pull complete
696bc3423178: Pull complete
0a65a0f5dcd1: Pull complete
d0bb45d42c8e: Pull complete
1a9616399b58: Pull complete
32cd32510543: Pull complete
6508acd2eecd: Pull complete
479ce1a1b674: Pull complete
6b6f799b42df: Pull complete
868e7e7207d7: Pull complete
aa9bdfcbcdf9: Pull complete
Digest: sha256:18b8a4745e2f1d1e655a021f5949441b46dcaf32d8edbcb149650b5137ebc250
Status: Downloaded newer image for dynverse/ti_paga_tree:v0.9.9.04
Executing 'paga_tree' on 'MyData'
With parameters: list(n_neighbors = 15L, n_comps = 50L, n_dcs = 15L, resolution = 1L,     embedding_type = "fa"),
inputs: counts, and
priors : start_id
Loading required namespace: hdf5r
Input saved to D:\Users\sandra.wiedenmann\AppData\Local\Temp\Rtmp6hdAbw\file3f782c6a2b54/ti
Running method using babelwhale
Running "C:\PROGRA~1\Docker\Docker\RESOUR~1\bin\docker.exe" run --name 20190812_125834__container__jby8fcpVMu -e \
  "TMPDIR=/tmp2" --workdir /ti/workspace -v \
  "/d/Users/sandra.wiedenmann/AppData/Local/Temp/Rtmp6hdAbw/file3f782c6a2b54/ti:/ti" -v \
  "/d/Users/sandra.wiedenmann/AppData/Local/Temp/Rtmp6hdAbw/file3f78189f64f5/tmp:/tmp2" \
  "dynverse/ti_paga_tree:v0.9.9.04" --dataset /ti/input.h5 --output /ti/output.h5 --use_priors all
|

after a few hours I stopped it.
when I let it once run over night it showed the Error:

Error: Error during trajectory inference Container was killed, possibly because it ran out of memory (error code 137)

I also tried other packages and a colleague of mine has exactly the same problem with another single cell data set.

Best
Sandra

the session info:

R version 3.6.0 (2019-04-26)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows >= 8 x64 (build 9200)

Matrix products: default

locale:
[1] LC_COLLATE=German_Germany.1252  LC_CTYPE=German_Germany.1252    LC_MONETARY=German_Germany.1252
[4] LC_NUMERIC=C                    LC_TIME=German_Germany.1252    

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
 [1] shiny_1.3.2           dplyr_0.8.3           Matrix_1.2-17         dyno_0.1.1            dynwrap_1.1.3        
 [6] dynplot_1.0.1         dynmethods_1.0.2      dynguidelines_1.0.0   dynfeature_1.0.0.9000 Seurat_3.0.1         

loaded via a namespace (and not attached):
  [1] backports_1.1.4     Hmisc_4.2-0         dyndimred_1.0.1     babelwhale_1.0.0    plyr_1.8.4         
  [6] igraph_1.2.4.1      lazyeval_0.2.2      sp_1.3-1            proxyC_0.1.5        splines_3.6.0      
 [11] listenv_0.7.0       ggplot2_3.2.0       digest_0.6.20       foreach_1.4.7       htmltools_0.3.6    
 [16] viridis_0.5.1       gdata_2.18.0        magrittr_1.5        checkmate_1.9.4     carrier_0.1.0      
 [21] cluster_2.0.8       ROCR_1.0-7          remotes_2.1.0       globals_0.12.4      readr_1.3.1        
 [26] RcppParallel_4.4.3  R.utils_2.9.0       dynutils_1.0.3      pdist_1.2           colorspace_1.4-1   
 [31] ggrepel_0.8.1       xfun_0.8            crayon_1.3.4        jsonlite_1.6        zeallot_0.1.0      
 [36] iterators_1.0.12    survival_2.44-1.1   zoo_1.8-6           ape_5.3             glue_1.3.1         
 [41] polyclip_1.10-0     gtable_0.3.0        future.apply_1.3.0  dynparam_1.0.0      scales_1.0.0       
 [46] bibtex_0.4.2        Rcpp_1.0.2          metap_1.1           viridisLite_0.3.0   xtable_1.8-4       
 [51] htmlTable_1.13.1    reticulate_1.13     bit_1.1-14          foreign_0.8-71      rsvd_1.0.2         
 [56] akima_0.6-2         SDMTools_1.1-221.1  Formula_1.2-3       tsne_0.1-3          htmlwidgets_1.3    
 [61] httr_1.4.1          FNN_1.1.3           gplots_3.0.1.1      RColorBrewer_1.1-2  acepack_1.4.1      
 [66] ica_1.0-2           pkgconfig_2.0.2     R.methodsS3_1.7.1   farver_1.1.0        nnet_7.3-12        
 [71] tidyselect_0.2.5    rlang_0.4.0         reshape2_1.4.3      later_0.8.0         munsell_0.5.0      
 [76] tools_3.6.0         cli_1.1.0           ranger_0.11.2       ggridges_0.5.1      stringr_1.4.0      
 [81] yaml_2.2.0          npsurv_0.4-0        bit64_0.9-7         processx_3.4.1      knitr_1.24         
 [86] fitdistrplus_1.0-14 tidygraph_1.1.2     caTools_1.17.1.2    purrr_0.3.2         RANN_2.6.1         
 [91] ggraph_1.0.2        pbapply_1.4-1       future_1.14.0       nlme_3.1-139        mime_0.7           
 [96] GA_3.2              R.oo_1.22.0         hdf5r_1.2.0         compiler_3.6.0      rstudioapi_0.10    
[101] plotly_4.9.0        png_0.1-7           testthat_2.2.1      lsei_1.2-0          tibble_2.1.3       
[106] tweenr_1.0.1        stringi_1.4.3       ps_1.3.0            lattice_0.20-38     shinyjs_1.0        
[111] vctrs_0.2.0         pillar_1.4.2        Rdpack_0.11-0       lmtest_0.9-37       data.table_1.12.2  
[116] cowplot_1.0.0       bitops_1.0-6        irlba_2.3.3         gbRd_0.4-11         patchwork_0.0.1    
[121] httpuv_1.5.1        R6_2.4.0            latticeExtra_0.6-28 promises_1.0.1      KernSmooth_2.23-15 
[126] gridExtra_2.3       codetools_0.2-16    MASS_7.3-51.4       gtools_3.8.1        assertthat_0.2.1   
[131] rje_1.9             shinyWidgets_0.4.8  sctransform_0.2.0   parallel_3.6.0      hms_0.5.0          
[136] grid_3.6.0          rpart_4.1-15        tidyr_0.8.3         Rtsne_0.15          ggforce_0.2.2      
[141] base64enc_0.1-3    

Selecting expression value for infer_trajectory

Dear authors,

Thanks for the great package.

I was wondering if there is a way to select for the expression values while using infer_trajectory Please find an example below;

Creating the dataset (sce referes to SingleCellExperiment object)

RawCount <- t(assay(sce, "logcounts"))
NormExprs <-  t(assay(sce, "norm_exprs"))

dataset <- wrap_expression(
  counts = RawCount,
  expression = NormExprs)

After selecting the guidelines with the shinyapp, proceed to infer trajectory

model <- infer_trajectory(dataset, "slingshot", verbose=TRUE")

Which returns the following.

Executing 'slingshot' on '20190226_084323__data_wrapper__SETXwUVAWV' With parameters: list(shrink = 1, reweight = TRUE, reassign = TRUE, thresh = 0.001, maxit = 10L, stretch = 2, smoother = "smooth.spline", shrink.method = "cosine") And inputs: counts, start_id

I was wondering, how could we selct for a different expression input for slingshot? Specifcally if I want to use the normalised expression that has been log-transformed. Therefore I would like to skip the default log transformation.

Many thanks again for the great package

Best
Zaki

wrap_expression inheritance

Hello!

Thank you for putting together this wonderful package. I ran through the test dataset "fibroblast_reprogramming_treutlein" without any issues.
However, when I attempt to use my own datasets - I keep getting this bug when I try to "wrap" the dataset:
cell_ids inherits from NULL not character.
I tried to input the correct data into the dataset object created and run a trajectory method, but I receive the same error: cell_ids inherits from NULL not character.

I have run a check on dataset$cell_ids to check the inheritance and it seems to pass:
inherits(dataset$cell_ids, 'character')
[1] TRUE

Any help or guidance would be very much appreciated!
Thanks!

Problem with GUI

Hi,

Whenever I try running the GUI in guidelines_shiny() I get this error over and over in the console:

image

It seems to mostly show up when I'm pressing a button, and the numbers vary based on which button I'm pressing. When I finally try to "close and use", the app never closes.

I've been able to get the tutorial to work using the command line guidelines() function, but I'm finding the documentation for answer_questions() quite opaque so I have no idea how to set options manually...

Thanks!
Geoff

Running Dyno on HPC Cluster

Hi all,

Thanks for an amazing package. Forgive me as I'm new to using containers, so I'm having some trouble getting the singularity containers to run on our scientific computing cluster. I have pulled all the TI method images to my local machine and then uploaded them to the cluster, but I can't get singularity to recognize that directory as the cache directory to run out of. When I use the babelwhale call as suggested (like below)

##Specify where to place singularity containers
singularity_config <- dynwrap::set_default_config(dynwrap::create_singularity_config(cache_dir = "/hpc/users/cohenp05/Comp1_Container/"))
babelwhale::set_default_config(singularity_config)

and point the singularity cache_dir to my directory, I get this error:

FATAL:   Unable to handle docker://alpine:3.7 uri: failed to get SHA of docker://alpine:3.7: pinging docker registry returned: Get https://registry-1.docker.io/v2/: dial tcp 52.22.67.152:443: connect: network is unreachable

Here is my sessionInfo:

R version 3.5.3 (2019-03-11)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: CentOS Linux 7 (Core)

Matrix products: default
BLAS/LAPACK: /hpc/packages/minerva-centos7/intel/parallel_studio_xe_2019/compilers_and_libraries_2019.0.117/linux/mkl/lib/intel64_lin/libmkl_gf_lp64.so

locale:
[1] C

attached base packages:
 [1] splines   stats4    parallel  stats     graphics  grDevices utils    
 [8] datasets  methods   base     

other attached packages:
 [1] rlist_0.4.6.1       forcats_0.4.0       stringr_1.4.0      
 [4] purrr_0.3.2         readr_1.3.1         tidyr_0.8.3        
 [7] tibble_2.1.1        tidyverse_1.2.1     dyno_0.1.1         
[10] dynwrap_1.1.3       dynplot_1.0.1       dynmethods_1.0.4   
[13] dynguidelines_1.0.0 dynfeature_1.0.0    RColorBrewer_1.1-2 
[16] data.table_1.12.2   monocle_2.10.1      DDRTree_0.1.5      
[19] irlba_2.3.3         VGAM_1.1-1          Biobase_2.42.0     
[22] BiocGenerics_0.28.0 Matrix_1.2-15       gridExtra_2.3      
[25] Seurat_3.0.2        dplyr_0.8.1         plyr_1.8.4         
[28] cowplot_0.9.4       ggplot2_3.1.1      

loaded via a namespace (and not attached):
  [1] reticulate_1.12       R.utils_2.8.0         tidyselect_0.2.5     
  [4] htmlwidgets_1.3       grid_3.5.3            combinat_0.0-8       
  [7] ranger_0.11.2         docopt_0.6.1          Rtsne_0.15           
 [10] munsell_0.5.0         codetools_0.2-16      ica_1.0-2            
 [13] future_1.13.0         withr_2.1.2           colorspace_1.4-1     
 [16] fastICA_1.2-1         knitr_1.23            rstudioapi_0.10      
 [19] ROCR_1.0-7            gbRd_0.4-11           listenv_0.7.0        
 [22] Rdpack_0.11-0         slam_0.1-45           polyclip_1.10-0      
 [25] farver_1.1.0          pheatmap_1.0.12       generics_0.0.2       
 [28] xfun_0.7              R6_2.4.0              GA_3.2               
 [31] rsvd_1.0.0            pdist_1.2             bitops_1.0-6         
 [34] assertthat_0.2.1      promises_1.0.1        SDMTools_1.1-221.1   
 [37] scales_1.0.0          ggraph_1.0.2          nnet_7.3-12          
 [40] gtable_0.3.0          babelwhale_0.0.0.9000 npsurv_0.4-0         
 [43] globals_0.12.4        processx_3.3.1        tidygraph_1.1.2      
 [46] rlang_0.3.4           lazyeval_0.2.2        acepack_1.4.1        
 [49] broom_0.5.2           checkmate_1.9.3       yaml_2.2.0           
 [52] reshape2_1.4.3        modelr_0.1.4          backports_1.1.4      
 [55] httpuv_1.5.1          Hmisc_4.2-0           tools_3.5.3          
 [58] gplots_3.0.1.1        ggridges_0.5.1        Rcpp_1.0.1           
 [61] base64enc_0.1-3       densityClust_0.3      ps_1.3.0             
 [64] rpart_4.1-15          pbapply_1.4-0         viridis_0.5.1        
 [67] dynparam_1.0.0        zoo_1.8-5             haven_2.1.0          
 [70] ggrepel_0.8.1         cluster_2.0.7-1       magrittr_1.5         
 [73] carrier_0.1.0         lmtest_0.9-37         RANN_2.6.1           
 [76] fitdistrplus_1.0-14   matrixStats_0.54.0    hms_0.4.2            
 [79] patchwork_0.0.1       lsei_1.2-0            mime_0.6             
 [82] xtable_1.8-4          readxl_1.3.1          sparsesvd_0.1-4      
 [85] HSMMSingleCell_1.2.0  testthat_2.1.1        compiler_3.5.3       
 [88] KernSmooth_2.23-15    crayon_1.3.4          rje_1.9              
 [91] R.oo_1.22.0           htmltools_0.3.6       proxyC_0.1.4         
 [94] later_0.8.0           Formula_1.2-3         RcppParallel_4.4.2   
 [97] lubridate_1.7.4       tweenr_1.0.1          MASS_7.3-51.1        
[100] cli_1.1.0             R.methodsS3_1.7.1     gdata_2.18.0         
[103] metap_1.1             igraph_1.2.4.1        pkgconfig_2.0.2      
[106] foreign_0.8-71        plotly_4.9.0          xml2_1.2.0           
[109] foreach_1.4.4         dynutils_1.0.3        rvest_0.3.4          
[112] bibtex_0.4.2          digest_0.6.19         dyndimred_1.0.1      
[115] sctransform_0.2.0     tsne_0.1-3            cellranger_1.1.0     
[118] htmlTable_1.13.1      shiny_1.3.2           gtools_3.8.1         
[121] nlme_3.1-137          jsonlite_1.6          viridisLite_0.3.0    
[124] limma_3.38.3          pillar_1.4.0          lattice_0.20-38      
[127] httr_1.4.0            survival_2.44-1.1     glue_1.3.1           
[130] remotes_2.0.4         qlcMatrix_0.9.7       FNN_1.1.3            
[133] png_0.1-7             iterators_1.0.10      ggforce_0.2.2        
[136] stringi_1.4.3         latticeExtra_0.6-28   caTools_1.17.1.2     
[139] future.apply_1.2.0    ape_5.3

I've also tried running the containers in a bash script by exporting my dynverse object and a dimensional reduction object to h5 files and running singularity run directly on the image like so:

ml singularity
module load R/3.5.3
SINGULARITY_CACHEDIR="/hpc/users/cohenp05/Comp1_Container/"
export SINGULARITY_CACHEDIR

singularity cache list

singularity run /hpc/users/cohenp05/Comp1_Container/ti_comp1_latest.sif --dataset=/hpc/users/cohenp05/Dyno_Objects/Monocyte_Dyno.h5 --output=/hpc/users/cohenp05/Dyno_Objects/Monocyte_Comp1.h5 --dimred=/hpc/users/cohenp05/Dyno_Objects/Monocytes_UMAP.h5

And then I get this error

Error: --dataset should contain a pathname of a .loom or .h5 file. Add a '-h' flag for help.
Execution halted

Please help! I'm out of ideas.

Error: 'manhattan_distance' is not an exported object from 'namespace:dynutils'

Hi,
thank you for this amazing package and study! I would like to use your package for my analysis, so I'm testing it.
I get this error when I try to make the plot with the following commands:

plot_dimred(
model,
expression_source = dataset$expression,
color_cells = "feature",
feature_oi = "Vim",
color_density = "grouping",
grouping = fibroblast_reprogramming_treutlein$grouping,
label_milestones = FALSE
)

Any help for this? Thanks!

installation of dyno fails since dynwrap is called twice

Hi the installation of the dyno package fails at the point where dynwrap is called again - @devel
dynfeature (NA -> 88274ff88...) [GitHub]
dynguidel... (NA -> d177eb719...) [GitHub]
dynmethods (NA -> fbac20a1a...) [GitHub]
dynplot (NA -> 87717ae06...) [GitHub]
dynwrap (NA -> dee7f681b...) [GitHub]
BH (NA -> 1.66.0-1 ) [CRAN]
bindr (NA -> 0.1.1 ) [CRAN]
bindrcpp (NA -> 0.2.2 ) [CRAN]
dplyr (NA -> 0.7.8 ) [CRAN]
fansi (NA -> 0.4.0 ) [CRAN]
pillar (NA -> 1.3.1 ) [CRAN]
pkgconfig (NA -> 2.0.2 ) [CRAN]
plogr (NA -> 0.2.0 ) [CRAN]
plyr (NA -> 1.8.4 ) [CRAN]
praise (NA -> 1.0.0 ) [CRAN]
purrr (NA -> 0.2.5 ) [CRAN]
ranger (NA -> 0.10.1 ) [CRAN]
RcppEigen (NA -> 0.3.3.5.0 ) [CRAN]
reshape2 (NA -> 1.4.3 ) [CRAN]
stringi (NA -> 1.2.4 ) [CRAN]
stringr (NA -> 1.3.1 ) [CRAN]
testthat (NA -> 2.0.1 ) [CRAN]
tibble (NA -> 1.4.2 ) [CRAN]
tidyr (NA -> 0.8.2 ) [CRAN]
tidyselect (NA -> 0.2.5 ) [CRAN]
utf8 (NA -> 1.1.4 ) [CRAN]
dynutils (NA -> 051f64ef5...) [GitHub]
dynwrap (NA -> dee7f681b...) [GitHub]
dyneval (NA -> 72948e12a...) [GitHub]
dynwrap (dee7f681b... -> 887a5b927...) [GitHub]

...
Downloading GitHub repo dynverse/dyneval@master
Installing package into โ€˜/usr/local/lib/R/site-libraryโ€™
(as โ€˜libโ€™ is unspecified)
* installing *source* package โ€˜dynevalโ€™ ...
** R
** data
*** moving datasets to lazyload DB
** inst
** preparing package for lazy loading
** help
*** installing help indices
*** copying figures
** building package indices
** testing if installed package can be loaded
* DONE (dyneval)
Downloading GitHub repo dynverse/dynwrap@devel
Skipping dynwrap, it is already being installed
Error in vapply(remotes, install_remote, ..., FUN.VALUE = character(1)) :
  values must be type 'character',
 but FUN(X[[2]]) result is type 'logical'
Calls: <Anonymous> ... update -> update.package_deps -> install_remotes -> vapply
Execution halted
...

static port for shiny app

Hi, I am experimenting with a docker-base setup of training sessions. Given that your package has been shown in our latest course, it tried to get the script running in a Jupyter notebook environment.
Installation in the docker goes fine but I have the idea that the port the shiny app is opened on e.g. 127.0.0.1:3880 varies each time. Would it be possible to use a static port?
Thanks for your feedback.

BTW: I am happy to explain more of the background via email if needed. It is still test phase and technology scouting.

Problem with singularity/docker?

Hi there,

I installed the package as recommended and I opted for singularity. When I run dynwrap::test_singularity_installation(detailed = TRUE), I get

โœ” Singularity is installed
โœ” Singularity is at correct version (>=3.0): 3.3.0-1 is installed
โœ” Singularity can pull and run a container from Dockerhub
โœ” Singularity can mount temporary volumes
โœ” Singularity test successful ------------------------------------------------------------
[1] TRUE

However, when I run model <- infer_trajectory(dataset, ti_comp1()) as in the user guide I get now the following error

Got permission denied while trying to connect to the Docker daemon socket at unix:///var/run/docker.sock: Get http://%2Fvar%2Frun%2Fdocker.sock/v1.40/version: dial unix /var/run/docker.sock: connect: permission denied
Error in processx::run("docker", c("images", paste0("--format=", format),  : 
  System command error

Any idea of what's happening?

installation issue for course

Hi,
When trying to install dyno on our servers for the single cell training I get this issue:
Cheers, Alex

library("devtools")
devtools::install_github("dynverse/dyno")
Skipping install of 'dyno' from a github remote, the SHA1 (707c424) has not changed since last install.
Use force = TRUE to force installation
devtools::install_github("dynverse/dyno", force=TRUE)
Downloading GitHub repo dynverse/dyno@master
checking for file \u2018/tmp/RtmpOAdEih/remotes17e265c31f69/dynverse-dyno-707c424/DESCRIPTION\u2019 ..\u2714 checking for file \u2018/tmp/RtmpOAdEih/remotes17e265c31f69/dynverse-dyno-707c424/DESCRIPTION\u2019
\u2500 preparing \u2018dyno\u2019:
\u2714 checking DESCRIPTION meta-information ...
\u2500 checking for LF line-endings in source and make files and shell scripts
\u2500 checking for empty or unneeded directories
\u2500 looking to see if a \u2018data/datalist\u2019 file should be added
\u2500 building \u2018dyno_0.1.0.tar.gz\u2019 (1.2s)

Installing package into \u2018/home/VIBTraining1/R/x86_64-pc-linux-gnu-library/3.5\u2019
(as \u2018lib\u2019 is unspecified)

  • installing source package \u2018dyno\u2019 ...
    ** R
    ** data
    *** moving datasets to lazyload DB
    ** byte-compile and prepare package for lazy loading
    Error: package or namespace load failed for \u2018dynfeature\u2019:
    object \u2018safe_tempdir\u2019 is not exported by 'namespace:dynutils'
    Error : package \u2018dynfeature\u2019 could not be loaded
    ERROR: lazy loading failed for package \u2018dyno\u2019
  • removing \u2018/home/VIBTraining1/R/x86_64-pc-linux-gnu-library/3.5/dyno\u2019
    Error in i.p(...) :
    (converted from warning) installation of package \u2018/tmp/RtmpOAdEih/file17e2110c1204/dyno_0.1.0.tar.gz\u2019 had non-zero exit status

sessionInfo()
R version 3.5.2 (2018-12-20)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu 18.04.1 LTS

Matrix products: default
BLAS: /usr/lib/x86_64-linux-gnu/blas/libblas.so.3.7.1
LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.7.1

locale:
[1] LC_CTYPE=C.UTF-8 LC_NUMERIC=C LC_TIME=C.UTF-8
[4] LC_COLLATE=C.UTF-8 LC_MONETARY=C.UTF-8 LC_MESSAGES=C.UTF-8
[7] LC_PAPER=C.UTF-8 LC_NAME=C LC_ADDRESS=C
[10] LC_TELEPHONE=C LC_MEASUREMENT=C.UTF-8 LC_IDENTIFICATION=C

attached base packages:
[1] stats graphics grDevices utils datasets methods base

loaded via a namespace (and not attached):
[1] Rcpp_1.0.0 plyr_1.8.4 pillar_1.3.1 compiler_3.5.2
[5] bindr_0.1.1 prettyunits_1.0.2 remotes_2.0.2 tools_3.5.2
[9] testthat_2.0.1 digest_0.6.18 pkgbuild_1.0.2 pkgload_1.0.2
[13] gtable_0.2.0 memoise_1.1.0 tibble_2.0.1 pkgconfig_2.0.2
[17] rlang_0.3.1 cli_1.0.1 rstudioapi_0.9.0 patchwork_0.0.1
[21] yaml_2.2.0 bindrcpp_0.2.2 withr_2.1.2 dplyr_0.7.8
[25] stringr_1.3.1 hms_0.4.2 fs_1.2.6 desc_1.2.0
[29] dynutils_1.0.0 devtools_2.0.1 grid_3.5.2 rprojroot_1.3-2
[33] tidyselect_0.2.5 glue_1.3.0 R6_2.3.0 processx_3.2.1
[37] sessioninfo_1.1.1 ggplot2_3.1.0 readr_1.3.1 tidyr_0.8.2
[41] purrr_0.3.0 callr_3.1.1 magrittr_1.5 scales_1.0.0
[45] usethis_1.4.0 backports_1.1.3 ps_1.3.0 dyndimred_0.1.0
[49] assertthat_0.2.0 colorspace_1.4-0 stringi_1.2.4 lazyeval_0.2.1
[53] munsell_0.5.0 crayon_1.3.4

Error Message: "Warning: Error in : Column `start_n` not found 128: <Anonymous>"

Hello,
I tried running the following code from your tutorial:

library(dyno)
library(tidyverse)
data("fibroblast_reprogramming_treutlein")
dim(fibroblast_reprogramming_treutlein$expression)
task <- wrap_expression(
counts = fibroblast_reprogramming_treutlein$counts,
expression = fibroblast_reprogramming_treutlein$expression)
guidelines <- guidelines_shiny(task)

The shiny app will launch, but when I try clicking the "Use Methods" button, I get the following error message: "Warning: Error in : Column start_n not found 128: " and then the window closes.

Do you have any suggestions?

Thank you,
Sarah

Error in plot_dimred function

I keep getting this error and am unable to resolve it.

I have tried plotting different trajectories and other arguments.

I have also tried both the most up to date version and previous versions of dynutils and the dyno package.

dynplot::plot_dimred(slingshot.total, color_cells = "pseudotime")

Error: 'manhattan_distance' is not an exported object from 'namespace:dynutils'

sessionInfo()

R version 3.5.3 (2019-03-11)
Platform: x86_64-apple-darwin15.6.0 (64-bit)
Running under: macOS Mojave 10.14.3

Matrix products: default
BLAS: >/System/Library/Frameworks/Accelerate.framework/Versions/A/Frameworks/vecLib.framework/Versions/A/libBLAS.dylib
LAPACK: /Library/Frameworks/R.framework/Versions/3.5/Resources/lib/libRlapack.dylib

locale:
[1] en_AU.UTF-8/en_AU.UTF-8/en_AU.UTF-8/C/en_AU.UTF-8/en_AU.UTF-8

attached base packages:
[1] stats graphics grDevices utils datasets methods base

other attached packages:
[1] forcats_0.4.0 stringr_1.4.0 dplyr_0.8.0.1 purrr_0.3.2 readr_1.3.1 tidyr_0.8.3
[7] tibble_2.1.1 tidyverse_1.2.1 dynutils_1.0.3 dyno_0.1.1 dynwrap_1.0.1 >dynplot_1.0.0.9000
[13] dynmethods_1.0.1 dynguidelines_1.0.0 dynfeature_1.0.0 Seurat_2.3.4 Matrix_1.2-17 >cowplot_0.9.4
[19] ggplot2_3.1.1

loaded via a namespace (and not attached):
[1] utf8_1.1.4 reticulate_1.12 R.utils_2.8.0 tidyselect_0.2.5 htmlwidgets_1.3 >grid_3.5.3
[7] trimcluster_0.1-2.1 ranger_0.11.2 Rtsne_0.15 devtools_2.0.2 dyntoy_0.9.9 >munsell_0.5.0
[13] codetools_0.2-16 ica_1.0-2 withr_2.1.2 colorspace_1.4-1 knitr_1.22 >rstudioapi_0.10
[19] stats4_3.5.3 ROCR_1.0-7 robustbase_0.93-4 dtw_1.20-1 gbRd_0.4-11 >labeling_0.3
[25] Rdpack_0.11-0 lars_1.2 polyclip_1.10-0 bit64_0.9-7 farver_1.1.0 >rprojroot_1.3-2
[31] generics_0.0.2 xfun_0.6 diptest_0.75-7 R6_2.4.0 GA_3.2 pdist_1.2
[37] hdf5r_1.2.0 flexmix_2.3-15 bitops_1.0-6 assertthat_0.2.1 promises_1.0.1 SDMTools_1.1-221.1
[43] scales_1.0.0 ggraph_1.0.2 nnet_7.3-12 gtable_0.3.0 babelwhale_0.0.0.9000 npsurv_0.4-0
[49] processx_3.3.0 tidygraph_1.1.2 rlang_0.3.4 akima_0.6-2 splines_3.5.3 lazyeval_0.2.2
[55] acepack_1.4.1 hexbin_1.27.2 broom_0.5.2 checkmate_1.9.3 modelr_0.1.4 BiocManager_1.30.4
[61] yaml_2.2.0 reshape2_1.4.3 backports_1.1.4 httpuv_1.5.1 Hmisc_4.2-0 tools_3.5.3
[67] usethis_1.5.0 gplots_3.0.1.1 RColorBrewer_1.1-2 proxy_0.4-23 sessioninfo_1.1.1 ggridges_0.5.1
[73] Rcpp_1.0.1 plyr_1.8.4 base64enc_0.1-3 ps_1.3.0 prettyunits_1.0.2 rpart_4.1-15
[79] pbapply_1.4-0 viridis_0.5.1 dynparam_1.0.0 zoo_1.8-5 haven_2.1.0 ggrepel_0.8.0
[85] cluster_2.0.9 fs_1.3.0 magrittr_1.5 RSpectra_0.14-0 data.table_1.12.2 carrier_0.1.0
[91] lmtest_0.9-37 RANN_2.6.1 mvtnorm_1.0-10 fitdistrplus_1.0-14 pkgload_1.0.2 hms_0.4.2
[97] patchwork_0.0.1 lsei_1.2-0 mime_0.6 xtable_1.8-4 readxl_1.3.1 mclust_5.4.3
[103] gridExtra_2.3 testthat_2.1.1 compiler_3.5.3 KernSmooth_2.23-15 crayon_1.3.4 rje_1.9
[109] R.oo_1.22.0 htmltools_0.3.6 proxyC_0.1.3 segmented_0.5-3.0 later_0.8.0 Formula_1.2-3
[115] snow_0.4-3 RcppParallel_4.4.2 lubridate_1.7.4 tweenr_1.0.1 MASS_7.3-51.4 fpc_2.1-11.1
[121] cli_1.1.0 R.methodsS3_1.7.1 gdata_2.18.0 parallel_3.5.3 metap_1.1 igraph_1.2.4.1
[127] pkgconfig_2.0.2 dyneval_0.9.9 sp_1.3-1 foreign_0.8-71 xml2_1.2.0 foreach_1.4.4
[133] vipor_0.4.5 rvest_0.3.3 bibtex_0.4.2 callr_3.2.0 digest_0.6.18 dyndimred_1.0.1
[139] RcppAnnoy_0.0.11 tsne_0.1-3 cellranger_1.1.0 htmlTable_1.13.1 uwot_0.1.3 curl_3.3
[145] kernlab_0.9-27 shiny_1.3.2 gtools_3.8.1 modeltools_0.2-22 nlme_3.1-139 jsonlite_1.6
[151] fansi_0.4.0 desc_1.2.0 viridisLite_0.3.0 pillar_1.3.1 lattice_0.20-38 shades_1.3.1
[157] httr_1.4.0 DEoptimR_1.0-8 pkgbuild_1.0.3 survival_2.44-1.1 glue_1.3.1 remotes_2.0.4
[163] FNN_1.1.3 png_0.1-7 prabclus_2.2-7 iterators_1.0.10 bit_1.1-14 ggforce_0.2.2
[169] class_7.3-15 stringi_1.4.3 mixtools_1.1.0 doSNOW_1.0.16 latticeExtra_0.6-28 caTools_1.17.1.2
[175] memoise_1.1.0 irlba_2.3.3 ape_5.3

Error: Error during trajectory inference - Error in normalise(counts) : trying to get slot "x" from an object of a basic class ("matrix") with no slots

Hi!
I ran into this issue when trying to infer a trajectory in one of your gold/datasets: gold/stimulated-dendritic-cells-LPS_shalek. I'm loading the dataset.rds which I assume contains everything that I need to run infer_trajectory() and other downstream analysis.

The error occurs when I use the infer_trajectory() using 'count' methods such as "monocle_ddrtree" or "paga_tree", as opposed to 'expression' methods such as "slingshot" and "comp1" (for which it works).

The error should be reproducible with the following code.

dataset_file <- c("real/gold/stimulated-dendritic-cells-LPS_shalek/dataset.rds")

## Import data
dataset <- read_rds(derived_file(dataset_file, "01-datasets"))

#c("slingshot", "monocle_ddrtree", "comp1", "paga_tree")
## Select Method
methods_selected <- "monocle_ddrtree" ## "paga_tree"
#methods_selected <- ti_monocle_ddrtree() ##  ti_paga_tree()

## Infer trajectory
trajectory <- infer_trajectory(dataset, methods_selected, verbose = TRUE)

which returns: (for "monocle_ddrtree")

Executing 'monocle_ddrtree' on 'real/gold/stimulated-dendritic-cells-LPS_shalek'
With parameters: list(reduction_method = "ICA", max_components = 2L, norm_method = "log",     auto_param_selection = TRUE, filter_features = TRUE, filter_features_mean_expression = 0.1),
inputs: counts, and
priors : 
Input saved to /tmp/RtmpCyjjsc/file88774270f8b5/ti
Running /usr/bin/docker run -e 'TMPDIR=/tmp2' --workdir /ti/workspace -v '/tmp/RtmpCyjjsc/file88774270f8b5/ti:/ti' -v \
  '/tmp/RtmpCyjjsc/file88772160a383/tmp:/tmp2' 'dynverse/ti_monocle_ddrtree:v0.9.9.01' --dataset /ti/input.h5 --output /ti/output.h5
Error in normalise(counts) : 
  trying to get slot "x" from an object of a basic class ("matrix") with no slots
Calls: <Anonymous> -> parse_dataset -> normalise
Execution halted
Error: Error during trajectory inference 
Error in normalise(counts) : 
  trying to get slot "x" from an object of a basic class ("matrix") with no slots
Calls: <Anonymous> -> parse_dataset -> normalise
Execution halted

Similarly for "paga_tree":

Executing 'paga_tree' on 'real/gold/stimulated-dendritic-cells-LPS_shalek'
With parameters: list(n_neighbors = 15L, n_comps = 50L, n_dcs = 15L, resolution = 1L,     embedding_type = "fa"),
inputs: counts, and
priors : start_id
Input saved to /tmp/RtmpCyjjsc/file887759e3f4e5/ti
Running /usr/bin/docker run -e 'TMPDIR=/tmp2' --workdir /ti/workspace -v '/tmp/RtmpCyjjsc/file887759e3f4e5/ti:/ti' -v \
  '/tmp/RtmpCyjjsc/file88773e675320/tmp:/tmp2' 'dynverse/ti_paga_tree:v0.9.9.01' --dataset /ti/input.h5 --output /ti/output.h5
R[write to console]: Error in normalise(counts) : 
  trying to get slot "x" from an object of a basic class ("matrix") with no slots
Calls: <Anonymous> ... <Anonymous> -> <Anonymous> -> parse_dataset -> normalise

Traceback (most recent call last):
  File "/code/run.py", line 22, in <module>
    task = dynclipy.main()
  File "/usr/local/lib/python3.7/site-packages/dynclipy/read.py", line 123, in main
    """)
  File "/usr/local/lib/python3.7/site-packages/rpy2/robjects/__init__.py", line 389, in __call__
    res = self.eval(p)
  File "/usr/local/lib/python3.7/site-packages/rpy2/robjects/functions.py", line 192, in __call__
    .__call__(*args, **kwargs))
  File "/usr/local/lib/python3.7/site-packages/rpy2/robjects/functions.py", line 121, in __call__
    res = super(Function, self).__call__(*new_args, **new_kwargs)
  File "/usr/local/lib/python3.7/site-packages/rpy2/rinterface_lib/conversion.py", line 28, in _
    cdata = function(*args, **kwargs)
  File "/usr/local/lib/python3.7/site-packages/rpy2/rinterface.py", line 773, in __call__
    raise embedded.RRuntimeError(_rinterface._geterrmessage())
rpy2.rinterface_lib.embedded.RRuntimeError: Error in normalise(counts) : 
  trying to get slot "x" from an object of a basic class ("matrix") with no slots
Calls: <Anonymous> ... <Anonymous> -> <Anonymous> -> parse_dataset -> normalise

Error: Error during trajectory inference 
R[write to console]: Error in normalise(counts) : 
  trying to get slot "x" from an object of a basic class ("matrix") with no slots
Calls: <Anonymous> ... <Anonymous> -> <Anonymous> -> parse_dataset -> normalise

Traceback (most recent call last):
  File "/code/run.py", line 22, in <module>
    task = dynclipy.main()
  File "/usr/local/lib/python3.7/site-packages/dynclipy/read.py", line 123, in main
    """)
  File "/usr/local/lib/python3.7/site-packages/rpy2/robjects/__init__.py", line 389, in __call__
    res = self.eval(p)
  File "/usr/local/lib/python3.7/site-packages/rpy2/robjects/functions.py", line 192, in __call__
    .__call__(*args, **kwargs))
  File "/usr/local/lib/python3.7/site-packages/rpy2/robjects/functions.py", line 121, in __call__
    res = super(Function, self).__call__(*new_args, **new_kwargs)
  File "/usr/local/lib/python3.7/site-packages/rpy2/rinterface_lib/conversion.py", line 28, in _
    cdata = function(*args, *

The error seems to be similar:

Error in normalise(counts) : 
  trying to get slot "x" from an object of a basic class ("matrix") with no slots
Calls: <Anonymous> ... <Anonymous> -> <Anonymous> -> parse_dataset -> normalise

I guess that it has to do with the fact that both methods are using the 'count' data instead of the 'expression' data, but I haven't found a way to solve this. I'm using your data wrong?

Thanks in advance!!

Error: 'random_seed' is not an exported object from 'namespace:dynutils' when running `infer_trajectory`

Whenever I ran infer_trajectory the same error message showed up
Error: 'random_seed' is not an exported object from 'namespace:dynutils'

Just wondering how can I fix it. I have tried dynutils_1.0.3 and dynutils_1.0.4.

sessionInfo()
R version 3.6.0 (2019-04-26)
Platform: x86_64-apple-darwin15.6.0 (64-bit)
Running under: macOS Mojave 10.14

Matrix products: default
BLAS: /System/Library/Frameworks/Accelerate.framework/Versions/A/Frameworks/vecLib.framework/Versions/A/libBLAS.dylib
LAPACK: /Library/Frameworks/R.framework/Versions/3.6/Resources/lib/libRlapack.dylib

locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8

attached base packages:
[1] stats graphics grDevices utils datasets methods base

other attached packages:
[1] forcats_0.4.0 stringr_1.4.0 dplyr_0.8.3 purrr_0.3.2 readr_1.3.1 tidyr_0.8.3 tibble_2.1.3
[8] ggplot2_3.2.1 tidyverse_1.2.1 shiny_1.3.2 dynutils_1.0.3 dyno_0.1.1 dynwrap_1.1.4.9000 dynplot_1.0.2.9000
[15] dynmethods_1.0.5 dynguidelines_1.0.0 dynfeature_1.0.0.9000

loaded via a namespace (and not attached):
[1] colorspace_1.4-1 rprojroot_1.3-2 dynparam_1.0.0.9000 fs_1.3.1 rstudioapi_0.10 rje_1.9 farver_1.1.0
[8] remotes_2.1.0 ggrepel_0.8.1 lubridate_1.7.4 xml2_1.2.2 ranger_0.11.2 codetools_0.2-16 polyclip_1.10-0
[15] pkgload_1.0.2 zeallot_0.1.0 jsonlite_1.6 broom_0.5.2 ggforce_0.3.1 compiler_3.6.0 httr_1.4.1
[22] backports_1.1.4 assertthat_0.2.1 Matrix_1.2-17 lazyeval_0.2.2 cli_1.1.0 later_0.8.0 tweenr_1.0.1
[29] htmltools_0.3.6 prettyunits_1.0.2 tools_3.6.0 igraph_1.2.4.1 gtable_0.3.0 glue_1.3.1 reshape2_1.4.3
[36] Rcpp_1.0.2 GA_3.2 cellranger_1.1.0 vctrs_0.2.0 nlme_3.1-141 iterators_1.0.12 ggraph_1.0.2
[43] ps_1.3.0 rvest_0.3.4 testthat_2.2.1 mime_0.7 akima_0.6-2 irlba_2.3.3 devtools_2.1.0
[50] MASS_7.3-51.4 scales_1.0.0 tidygraph_1.1.2 babelwhale_1.0.0.9000 hms_0.5.0 promises_1.0.1 RColorBrewer_1.1-2
[57] yaml_2.2.0 curl_4.0 memoise_1.1.0 gridExtra_2.3 dyndimred_1.0.2.9000 stringi_1.4.3 desc_1.2.0
[64] foreach_1.4.7 pkgbuild_1.0.4 rlang_0.4.0 pkgconfig_2.0.2 lattice_0.20-38 patchwork_0.0.1 cowplot_1.0.0
[71] processx_3.4.1 tidyselect_0.2.5 plyr_1.8.4 magrittr_1.5 R6_2.4.0 generics_0.0.2 pillar_1.4.2
[78] haven_2.1.1 carrier_0.1.0 withr_2.1.2 proxyC_0.1.5 sp_1.3-1 modelr_0.1.5 crayon_1.3.4
[85] shinyWidgets_0.4.8 viridis_0.5.1 usethis_1.5.1 grid_3.6.0 readxl_1.3.1 callr_3.3.1 digest_0.6.20
[92] xtable_1.8-4 httpuv_1.5.1 RcppParallel_4.4.3 munsell_0.5.0 viridisLite_0.3.0 sessioninfo_1.1.1 shinyjs_1.0

How to prepare the input dataset?

Hi Zouter,
Thank you so much for your quick response and help. The tutorial runs smoothly in my hand finally.
I am trying to run some real datasets which download from GEO. Do you have some pieces of instruction on how to prepare the input dataset from digital expression reads, such as Drop-seq output? That would be very help. thanks a lot.

Y

HDF5 Errors when attempting to run trajectory inference

Fantastic effort on this package, I'm really looking forward to being able to use TI methods from a common interface in R. I've followed the installation instructions and successfully installed the packages and passed the docker installation test, but I'm having difficulty running the code with the example data or my own data.

For example, running the ti_comp1() on the example_dataset, I get an error on related to HDF5.

> library(dyno)
> model <- infer_trajectory(example_dataset, ti_comp1(), verbose = T)
Executing 'comp1' on 'example'
With parameters: list(dimred = "pca", ndim = 2L, component = 1L),
inputs: expression, and
priors : 
Input saved to C:\Users\reismab\AppData\Local\Temp\RtmpOK0ZZ0\file19089056d27/ti
Running "C:\PROGRA~1\Docker\Docker\RESOUR~1\bin\docker.exe" run -e "TMPDIR=/tmp2" --workdir /ti/workspace -v \
  "/c/Users/reismab/AppData/Local/Temp/RtmpOK0ZZ0/file19089056d27/ti:/ti" -v \
  "/c/Users/reismab/AppData/Local/Temp/RtmpOK0ZZ0/file19083e865acc/tmp:/tmp2" "dynverse/ti_comp1:v0.9.9" --dataset \
  /ti/input.h5 --output /ti/output.h5
Loading required package: dynutils
Error: Error during trajectory inference 
HDF5-API Errors:
    error #000: C:\pkg\hdf5-1.8.14\src\H5Tnative.c in H5Tget_native_type(): line 119: cannot retrieve native type
        class: HDF5
        major: Invalid arguments to routine
        minor: Inappropriate type

    error #001: C:\pkg\hdf5-1.8.14\src\H5Tnative.c in H5T_get_native_type(): line 400: cannot get member value
        class: HDF5
        major: Invalid arguments to routine
        minor: Inappropriate type

    error #002: C:\pkg\hdf5-1.8.14\src\H5T.c in H5T_convert(): line 4816: data type conversion failed
        class: HDF5
        major: Attribute
        minor: Unable to encode value

    error #003: C:\pkg\hdf5-1.8.14\src\H5Tconv.c in H5T__conv_i_i(): line 3639: can't find property list for ID
        class: HDF5
        major: Object atom
        minor: Unable to find atom information (already closed?)

    error #004: C:\pkg\hdf5-1.8.14\src\H5Pint.c in H5P_object_verify(): line 3381: property list is no

I get similar errors when trying to run ti_PAGA() or ti_wanderlust() with the example datasets as well as the fibroblast_reprogramming_treutlein dataset. I've also tried restarting R, docker, and windows without success. I'm not sure if this is a bug, or [perhaps more likely] an issue with my installation of the packages/docker. Unfortunately, I don't know how to begin to debug this but I've included my session info below.

> sessionInfo()
R version 3.5.3 (2019-03-11)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows >= 8 x64 (build 9200)

Matrix products: default

locale:
[1] LC_COLLATE=English_United States.1252  LC_CTYPE=English_United States.1252    LC_MONETARY=English_United States.1252
[4] LC_NUMERIC=C                           LC_TIME=English_United States.1252    

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] shiny_1.2.0         dyno_0.9.9          dynwrap_1.0.0       dynplot_1.0.0       dynmethods_1.0.0    dynguidelines_1.0.0
[7] dynfeature_1.0.0   

loaded via a namespace (and not attached):
  [1] colorspace_1.4-1      rprojroot_1.3-2       dynparam_1.0.0        htmlTable_1.13.1      base64enc_0.1-3      
  [6] fs_1.2.7              rstudioapi_0.10       rje_1.9               farver_1.1.0          remotes_2.0.2        
 [11] dynutils_1.0.2        bit64_0.9-7           ggrepel_0.8.0         ranger_0.11.2         codetools_0.2-16     
 [16] splines_3.5.3         knitr_1.22            polyclip_1.10-0       pkgload_1.0.2         Formula_1.2-3        
 [21] jsonlite_1.6          cluster_2.0.7-1       ggforce_0.2.1         readr_1.3.1           compiler_3.5.3       
 [26] backports_1.1.3       assertthat_0.2.1      Matrix_1.2-15         lazyeval_0.2.2        cli_1.1.0            
 [31] later_0.8.0           tweenr_1.0.1          acepack_1.4.1         htmltools_0.3.6       prettyunits_1.0.2    
 [36] tools_3.5.3           igraph_1.2.4          gtable_0.3.0          glue_1.3.1            reshape2_1.4.3       
 [41] dplyr_0.8.0.1         Rcpp_1.0.1            GA_3.2                iterators_1.0.10      ggraph_1.0.2         
 [46] xfun_0.6              stringr_1.4.0         ps_1.3.0              testthat_2.0.1        akima_0.6-2          
 [51] mime_0.6              devtools_2.0.1        MASS_7.3-51.1         scales_1.0.0          tidygraph_1.1.2      
 [56] babelwhale_0.0.0.9000 hms_0.4.2             promises_1.0.1        RColorBrewer_1.1-2    yaml_2.2.0           
 [61] curl_3.3              memoise_1.1.0         gridExtra_2.3         ggplot2_3.1.0         dyndimred_1.0.0      
 [66] rpart_4.1-13          latticeExtra_0.6-28   stringi_1.4.3         desc_1.2.0            foreach_1.4.4        
 [71] checkmate_1.9.1       shades_1.3.1          pkgbuild_1.0.3        rlang_0.3.3           pkgconfig_2.0.2      
 [76] lattice_0.20-38       purrr_0.3.2           patchwork_0.0.1       htmlwidgets_1.3       bit_1.1-14           
 [81] cowplot_0.9.4         pdist_1.2             processx_3.3.0        tidyselect_0.2.5      plyr_1.8.4           
 [86] magrittr_1.5          R6_2.4.0              Hmisc_4.2-0           pillar_1.3.1          foreign_0.8-71       
 [91] carrier_0.1.0         withr_2.1.2           sp_1.3-1              survival_2.43-3       nnet_7.3-12          
 [96] tibble_2.1.1          crayon_1.3.4          hdf5r_1.1.1           shinyWidgets_0.4.8    viridis_0.5.1        
[101] usethis_1.4.0         grid_3.5.3            data.table_1.12.0     callr_3.2.0           digest_0.6.18        
[106] xtable_1.8-3          tidyr_0.8.3           httpuv_1.5.0          munsell_0.5.0         viridisLite_0.3.0    
[111] shinyjs_1.0           sessioninfo_1.1.1  
> dynwrap::test_docker_installation(detailed = TRUE)
<U+2714> Docker is installed
<U+2714> Docker daemon is running
<U+2714> Docker is at correct version (>1.0): 1.39
<U+2714> Docker is in linux mode
<U+2714> Docker can pull images
<U+2714> Docker can run image
<U+2714> Docker can mount temporary volumes
<U+2714> Docker test successful -----------------------------------------------------------------
[1] TRUE

dyno installation fails due to API limit rate

Hi,
Since one week or so, I experience issues when installing dyno via an ansible script on a google cloud server but also in an interactive R session on the same server.
This has not been the case in the summer when I installed your tool for the summer school.
Any ideas?
Thanks a lot.
Cheers, Alex

R version 3.4.4 (2018-03-15) -- "Someone to Lean On"
Copyright (C) 2018 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> devtools::install_github('dynverse/dyno')
Error: HTTP error 403.
  API rate limit exceeded for 35.195.YY.XX. (But here's the good news: Authenticated requests get a higher rate lim
it. Check out the documentation for more details.)
  Rate limit remaining: 0/60
  Rate limit reset at: 2018-12-16 13:41:48 UTC
  To increase your GitHub API rate limit
  - Use `usethis::browse_github_pat()` to create a Personal Access Token.
  - Use `usethis::edit_r_environ()` and add the token as `GITHUB_PAT`.
> traceback()
8: stop(github_error(res))
7: github_DESCRIPTION(username = remote$username, repo = remote$repo, 
       subdir = remote$subdir, host = remote$host, ref = remote$ref, 
       pat = remote$auth_token %||% github_pat(), use_curl = use_curl)
6: remote_package_name.github_remote(remote)
5: remote_package_name(remote)
4: FUN(X[[i]], ...)
3: vapply(remotes, install_remote, ..., FUN.VALUE = character(1))
2: install_remotes(remotes, auth_token = auth_token, host = host, 
       dependencies = dependencies, upgrade = upgrade, force = force, 
       quiet = quiet, build = build, build_opts = build_opts, repos = repos, 
       type = type, ...)
1: devtools::install_github("dynverse/dyno")

error when run the command in tutorial

Hi zouter,
when I run the commands in the tutorial, I got another error.
library(dyno)
library(tidyverse)
data("fibroblast_reprogramming_treutlein")
dim(fibroblast_reprogramming_treutlein$expression)
task <- wrap_expression(
counts = fibroblast_reprogramming_treutlein$counts,
expression = fibroblast_reprogramming_treutlein$expression)
guidelines <- guidelines_shiny(task)

Error in setdiff(names(questions), given_question_ids) :
object 'questions' not found

Usage of add_prior_information

Hi,

I am trying to use dyno to run PAGA, which requires the user to define one or multiple starting cells. I tried using add_prior_information to add a vector of cell names that should be used as root cells, and the command appeared to accept that vector as valid input. However, when I try to run PAGA I keep getting this error as if the list of starting cells was not actually added to the dataset object. Any pointers on how to properly add metadata to the dyno objects? Thanks!

Running:

#Annotate starting cells
start_cells = as.vector([email protected])
#Add prior information to dataset
add_prior_information(dataset = dataset_object, start_id = start_cells)
#Run model
model_PAGA <- infer_trajectory(dataset = dataset_object, method = methods_selected[2], verbose = T, give_priors = "start_id")

Gives the error:
Error in .method_extract_args(dataset, method$inputs, give_priors) :
Prior information start_id is missing from dataset 20190219_155152__data_wrapper__9eeKucfMOk but is required by the method.
-> If known, you can add this prior information using add_prior_information(dataset, start_id = ).
-> Otherwise, this method cannot be used.

Cannot open connection

when executing guidelines <- guidelines_shiny(task) I get the error:
Error in file(con, "r") :
cannot open the connection to 'https://www.googletagmanager.com/gtag/js?id = UA-578149-3'
In addition: Warning message:
In file(con, "r") :

This website successfully opens when I type it into a browser so I'm at a loss for why it cannot open in R. Is this an error other users have encountered?

infer_trajectory: Error in readRDS(file) : unknown input format

Dear dyno-team,

while dyno worked for me previously on macosx + docker,
i am having trouble getting the package to work on a linux server using singularity

containers are not being pulled, although singularity check is fine:


> dynwrap::test_singularity_installation(detailed = TRUE)
โœ” Singularity is installed
โœ” Singularity is at correct version (>=3.0): 3.1.0 is installed
INFO:    Converting OCI blobs to SIF format
INFO:    Starting build...
Getting image source signatures
Copying blob sha256:5d20c808ce198565ff70b3ed23a991dd49afac45dece63474b27ce6ed036adc6
 2.01 MiB / 2.01 MiB  0s
Copying config sha256:54c9afe02933973b548088d74fca8dc95f0c3ae0d1bc765efbf2673c31a8e627
 585 B / 585 B  0s
Writing manifest to image destination
Storing signatures
INFO:    Creating SIF file...
INFO:    Build complete: 
INFO:    Image cached as SIF at 
โœ” Singularity can pull and run a container from Dockerhub
โœ” Singularity can mount temporary volumes
โœ” Singularity test successful ------------------------------------------------------------
[1] TRUE

another singularity test

system("singularity exec 'docker://dynverse/ti_slingshot:v0.9.9' echo hi")
INFO:    Converting OCI blobs to SIF format
INFO:    Starting build...
Getting image source signatures
Copying blob sha256:ab1fc7e4bf9195e554669fafa47f69befe22053d7100f5f7002cb9254a36f37c
...

INFO:    Creating SIF file...
INFO:    Build complete: 
...
INFO:    Image cached as SIF at 

nevertheless i get an error when trying to use trajectory inference:

dataset <- wrap_expression(
  counts = example_dataset$counts,
  expression = example_dataset$expression
)

model <- infer_trajectory(dataset,c("slingshot"))
Error in readRDS(file) : unknown input format

model <- infer_trajectory(dataset,"slingshot", verbose = TRUE, debug = TRUE)
Error in readRDS(file) : unknown input format

model <- infer_trajectory(dataset, ti_slingshot(), verbose = TRUE)
Error in readRDS(file) : unknown input format

create_ti_method_container()
Error in readRDS(file) : unknown input format

ti_slingshot()
Error in readRDS(file) : unknown input format
R version 3.5.1 (2018-07-02)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: CentOS Linux 7 (Core)

other attached packages:
[1] dyno_0.1.1            dynwrap_1.1.4.9000    dynplot_1.0.2.9000   
[4] dynmethods_1.0.5      dynguidelines_1.0.0   dynfeature_1.0.0.9000

loaded via a namespace (and not attached):
  [1] fs_1.3.1             usethis_1.5.1        dynutils_1.0.4.9000 
  [4] devtools_2.2.0       RColorBrewer_1.1-2   rprojroot_1.3-2     
  [7] tools_3.5.1          backports_1.1.4      utf8_1.1.4          
 [10] R6_2.4.0             irlba_2.3.3          DT_0.8              
 [13] lazyeval_0.2.2       colorspace_1.4-1     sp_1.3-1            
 [16] withr_2.1.2          tidyselect_0.2.5     carrier_0.1.0       
 [19] gridExtra_2.3        prettyunits_1.0.2    processx_3.4.1      
 [22] proxyC_0.1.5         curl_4.0             compiler_3.5.1      
 [25] cli_1.1.0            desc_1.2.0           scales_1.0.0        
 [28] readr_1.3.1          callr_3.3.1          stringr_1.4.0       
 [31] digest_0.6.20        pkgconfig_2.0.2      htmltools_0.3.6     
 [34] akima_0.6-2          sessioninfo_1.1.1    GA_3.2              
 [37] rje_1.10.10          htmlwidgets_1.3      rlang_0.4.0         
 [40] shiny_1.3.2          farver_1.1.0         dplyr_0.8.3         
 [43] magrittr_1.5         patchwork_0.0.1      Matrix_1.2-14       
 [46] fansi_0.4.0          Rcpp_1.0.2           dyndimred_1.0.2.9000
 [49] munsell_0.5.0        viridis_0.5.1        stringi_1.4.3       
 [52] yaml_2.2.0           ggraph_2.0.0         MASS_7.3-50         
 [55] pkgbuild_1.0.5       plyr_1.8.4           grid_3.5.1          
 [58] promises_1.0.1       ggrepel_0.8.1        crayon_1.3.4        
 [61] lattice_0.20-35      cowplot_1.0.0        graphlayouts_0.5.0  
 [64] hms_0.5.1            knitr_1.24           zeallot_0.1.0       
 [67] ps_1.3.0             pillar_1.4.2         igraph_1.2.4.1      
 [70] ranger_0.11.2        babelwhale_1.0.1     codetools_0.2-15    
 [73] reshape2_1.4.3       pkgload_1.0.2        glue_1.3.1          
 [76] remotes_2.1.0        RcppParallel_4.4.3   foreach_1.4.7       
 [79] vctrs_0.2.0          dynparam_1.0.0.9000  tweenr_1.0.1        
 [82] httpuv_1.5.1         testthat_2.2.1       gtable_0.3.0        
 [85] purrr_0.3.2          polyclip_1.10-0      tidyr_0.8.3         
 [88] assertthat_0.2.1     ggplot2_3.2.1        xfun_0.9            
 [91] ggforce_0.3.1        lmds_0.1.0           mime_0.7            
 [94] xtable_1.8-4         tidygraph_1.1.2      later_0.8.0         
 [97] viridisLite_0.3.0    tibble_2.1.3         iterators_1.0.12    
[100] memoise_1.1.0        ellipsis_0.2.0.1    

thankful for any suggestions!

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.