In this vignette it is shown the analysis performed on a single cell RNA seq dataset from a time-course of iPS reprogramming. The data are from day 10, day 12, day 14, day 16 and day 18 from Schiebinger et al. 2019. Below the transition matrices between time points t and t+1 are loaded. Each matrix is the output of \emph{get_transition_matrix} and was built starting from \emph{compute_all_transport_maps} function from pyhton package . See ?\emph{get_transition_matrix} for more info.
load(system.file("extdata", "cluster_label_example.Rda", package = "WOTPLY"))
#load(system.file("extdata", "time_label_example.Rda", package = "WOTPLY"))
load(system.file("extdata", "example_day_10_12.Rda", package = "WOTPLY"))
load(system.file("extdata", "example_day_12_14.Rda", package = "WOTPLY"))
load(system.file("extdata", "example_day_14_16.Rda", package = "WOTPLY"))
load(system.file("extdata", "example_day_16_18.Rda", package = "WOTPLY"))list_transition_matrices_example <- list(example_day_10_12,example_day_12_14,example_day_14_16,example_day_16_18)
selected_stages <- c("IPS","Trophoblast","Epithelial","STROMAL","NEURAL")
legend_time_example <- c("day_10","day_12","day_14","day_16","day_18")
customize_color <- WOTPLY:::gg_color_hue(length(levels(factor(cluster_label_example))))The output of WOTPLY function shows the connections between \emph{selected_stages} from the latest time point and the clusters from previous time points. The number of columns is equal to the numbers of time points. In each column, the cluster of the corresponding time point is shown as network node. The weight of the links between clusters at time points t and t+1 reflect the weight of the transition probabilities from \emph{list_transition_matrices}. It is possible to change the maximum number of links to select between clusters at time t and clusters at time t+1. Links are sorted according to the weight and then only the \emph{top_link} are kept. If (default), all the links are kept.
#png("/Users/gabriele.lubatti/Desktop/Phd/Embryo_Organoids/wot_publsihed_data/input_wot/WOTPLY_1.png")
WOTPLY(list_transition_matrices_example, selected_stages, cluster_label_example, legend_time_example, customize_color, top_link = NULL)
#> Registered S3 method overwritten by 'GGally':
#> method from
#> +.gg ggplot2
#png("/Users/gabriele.lubatti/Desktop/Phd/Embryo_Organoids/wot_publsihed_data/input_wot/WOTPLY_2.png")
WOTPLY(list_transition_matrices_example, selected_stages, cluster_label_example, legend_time_example,customize_color, top_link = 3)
#> Warning in WOTPLY(list_transition_matrices_example, selected_stages,
#> cluster_label_example, : From time day_10 to day_12 is only possible to select
#> less than 3 connections
utils::sessionInfo()
#> R version 4.0.2 (2020-06-22)
#> Platform: x86_64-apple-darwin17.0 (64-bit)
#> Running under: macOS Mojave 10.14.6
#>
#> Matrix products: default
#> BLAS: /Library/Frameworks/R.framework/Versions/4.0/Resources/lib/libRblas.dylib
#> LAPACK: /Library/Frameworks/R.framework/Versions/4.0/Resources/lib/libRlapack.dylib
#>
#> locale:
#> [1] C/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] WOTPLY_0.1.0
#>
#> loaded via a namespace (and not attached):
#> [1] Rcpp_1.0.7 sna_2.6 highr_0.9
#> [4] plyr_1.8.6 pillar_1.6.4 bslib_0.3.1
#> [7] compiler_4.0.2 RColorBrewer_1.1-2 jquerylib_0.1.4
#> [10] tools_4.0.2 digest_0.6.29 jsonlite_1.7.2
#> [13] evaluate_0.14 lifecycle_1.0.1 tibble_3.1.6
#> [16] gtable_0.3.0 lattice_0.20-41 pkgconfig_2.0.3
#> [19] rlang_0.4.12 DBI_1.1.2 GGally_2.1.2
#> [22] yaml_2.2.1 xfun_0.29 fastmap_1.1.0
#> [25] coda_0.19-4 dplyr_1.0.7 stringr_1.4.0
#> [28] knitr_1.37 generics_0.1.1 sass_0.4.0
#> [31] vctrs_0.3.8 tidyselect_1.1.1 grid_4.0.2
#> [34] reshape_0.8.8 glue_1.6.0 R6_2.5.1
#> [37] fansi_0.5.0 rmarkdown_2.11 farver_2.1.0
#> [40] purrr_0.3.4 ggplot2_3.3.5 magrittr_2.0.1
#> [43] scales_1.1.1 ellipsis_0.3.2 htmltools_0.5.2
#> [46] assertthat_0.2.1 colorspace_2.0-2 utf8_1.2.2
#> [49] stringi_1.7.6 network_1.17.1 munsell_0.5.0
#> [52] statnet.common_4.5.0 crayon_1.4.2