Author: Alan O’Callaghan (alan.b.ocallaghan@gmail.com)

1 Introduction

Due to the size of the objects, this file is seperated to three fils. You can view this series in the following links:

This is file 3 in the series.

# Let's load the packages
library(heatmaply)
library(heatmaplyExamples)

2 Visualization of voom-transformed data (median-centered data, PAM50 genes only)

Following normalization, gene expression patterns appear roughly similar. This indicates that relative expression levels have not been altered unduly. Furthermore, slightly increased concordance with the pre-assigned cluster labels is observed in normalized data. Samples appear to cluster based Sample-sample correlation appears to show less concordance with row annotations than clustering based on gene expression. However, the use of different linkage criteria or distance measures may alter the observed clusters.

center_voom_mat <- voomed_pam50_expression - 
    apply(voomed_pam50_expression, 1, median)

voom_max <- max(abs(center_voom_mat))
voom_limits <- c(-voom_max, voom_max)


heatmaply(t(center_voom_mat), 
    row_side_colors=tcga_brca_clinical,
    fontsize_col = 7.5,
    showticklabels = c(TRUE, FALSE),
    col = cool_warm(50),
    limits = voom_limits,
    main = 'Normalised, centred log2 CPM, PAM50 genes',
    plot_method = 'plotly')
heatmaply_cor(cor(center_voom_mat), 
    row_side_colors = tcga_brca_clinical,
    showticklabels = c(FALSE, FALSE),
    main = 'Sample-sample correlation based on centred, normalised PAM50 gene expression',
    plot_method = 'plotly')

3 Discussion

It may be useful when examining expression heatmaps to identify particularly high or low measures for a single gene in a group of patients, or a gene which shows unusually high or low variance. The mouse-over text available in the heatmaply package allows visual assessment of measures of interest and quick identification of samples or genes with unusual gene expression patterns. Similarly, visualizing correlation heatmaps with heatmaply allows the user to rapidly identify samples with unusually high or low pairwise correlation.

4 References

5 sessionInfo

sessionInfo()
#> R version 3.4.1 (2017-06-30)
#> Platform: x86_64-w64-mingw32/x64 (64-bit)
#> Running under: Windows 7 x64 (build 7601) Service Pack 1
#> 
#> Matrix products: default
#> 
#> locale:
#> [1] LC_COLLATE=Hebrew_Israel.1255  LC_CTYPE=Hebrew_Israel.1255   
#> [3] LC_MONETARY=Hebrew_Israel.1255 LC_NUMERIC=C                  
#> [5] LC_TIME=Hebrew_Israel.1255    
#> 
#> attached base packages:
#> [1] stats     graphics  grDevices datasets  utils     methods   base     
#> 
#> other attached packages:
#>  [1] bindrcpp_0.2            knitr_1.16             
#>  [3] heatmaplyExamples_0.2.0 heatmaply_0.11.0       
#>  [5] viridis_0.4.0           viridisLite_0.2.0      
#>  [7] plotly_4.7.1            ggplot2_2.2.1.9000     
#>  [9] installr_0.19.0         stringr_1.2.0          
#> 
#> loaded via a namespace (and not attached):
#>  [1] httr_1.2.1         tidyr_0.6.3        jsonlite_1.5      
#>  [4] foreach_1.4.3      gtools_3.5.0       shiny_1.0.3       
#>  [7] assertthat_0.2.0   stats4_3.4.1       yaml_2.1.14       
#> [10] robustbase_0.92-7  backports_1.1.0    lattice_0.20-35   
#> [13] glue_1.1.1         digest_0.6.12      RColorBrewer_1.1-2
#> [16] colorspace_1.3-2   httpuv_1.3.5       htmltools_0.3.6   
#> [19] plyr_1.8.4         pkgconfig_2.0.1    devtools_1.13.2   
#> [22] xtable_1.8-2       purrr_0.2.2.2      mvtnorm_1.0-6     
#> [25] scales_0.5.0       gdata_2.18.0       whisker_0.3-2     
#> [28] tibble_1.3.4       withr_2.0.0        nnet_7.3-12       
#> [31] lazyeval_0.2.0     mime_0.5           magrittr_1.5      
#> [34] mclust_5.3         memoise_1.1.0      evaluate_0.10.1   
#> [37] MASS_7.3-47        gplots_3.0.1       class_7.3-14      
#> [40] tools_3.4.1        registry_0.3       data.table_1.10.4 
#> [43] trimcluster_0.1-2  kernlab_0.9-25     munsell_0.4.3     
#> [46] cluster_2.0.6      fpc_2.1-10         compiler_3.4.1    
#> [49] caTools_1.17.1     rlang_0.1.2        grid_3.4.1        
#> [52] iterators_1.0.8    htmlwidgets_0.9    crosstalk_1.0.0   
#> [55] bitops_1.0-6       rmarkdown_1.6      gtable_0.2.0      
#> [58] codetools_0.2-15   flexmix_2.3-14     TSP_1.1-5         
#> [61] R6_2.2.2           seriation_1.2-2    gridExtra_2.2.1   
#> [64] prabclus_2.2-6     dplyr_0.7.3.9000   bindr_0.1         
#> [67] rprojroot_1.2      KernSmooth_2.23-15 dendextend_1.6.0  
#> [70] modeltools_0.2-21  stringi_1.1.5      Rcpp_0.12.12      
#> [73] gclus_1.3.1        DEoptimR_1.0-8     diptest_0.75-7