DEP package: error in plot_heatmap
Entering edit mode
Rob • 0
Last seen 6 months ago

I am using the DEP package from Bioconductor for some proteomics analysis. I found a bug in the plotheatmap function: it works well when the argument plot=TRUE (default), but it does return an error when plot= FALSE. I need this information to retrieve the coordinates to get the Genenames from my different clusters. Example from the default of DEP package:

data <- UbiLength
data_unique <- make_unique(data, "Gene.names", "Protein.IDs", delim = ";")
columns <- grep("LFQ.", colnames(data_unique))
exp_design <- UbiLength_ExpDesign
se <- make_se(data_unique, columns, exp_design)
filt <- filter_missval(se, thr = 0)
imputed <- impute(filt, fun = "MinProb", q = 0.01)
diff <- test_diff(imputed, "control", "Ctrl")
dep <- add_rejections(diff, alpha = 0.05, lfc = 1)
plot_heatmap(dep, 'centered', kmeans = TRUE, k = 6, row_font_size = 3, plot=TRUE) #all good until here
plot_heatmap(dep, 'centered', kmeans = TRUE, k = 6, row_font_size = 3, plot=FALSE)

The last command will return the error:

  Error in rownames_to_column(., var = "protein") : is not TRUE

Can I do something to solve this issue?

DEP • 416 views
Entering edit mode

Dear robert.arsie,

Thank you for pointing out the bug. I am working on it and will let you know once a fix is available.

Best, Arne

Entering edit mode
> sessionInfo()
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

[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] imputeLCMD_2.0      impute_1.58.0       pcaMethods_1.76.0   Biobase_2.44.0      BiocGenerics_0.30.0 norm_1.0-9.5        tmvtnorm_1.4-10    
 [8] gmm_1.6-2           sandwich_2.5-1      Matrix_1.2-17       mvtnorm_1.0-10      DEP_1.6.0          

loaded via a namespace (and not attached):
 [1] ProtGenerics_1.16.0         bitops_1.0-6                matrixStats_0.54.0          doParallel_1.0.14           RColorBrewer_1.1-2         
 [6] GenomeInfoDb_1.20.0         MSnbase_2.10.0              tools_3.6.0                 R6_2.4.0                    DT_0.6                     
[11] affyio_1.54.0               lazyeval_0.2.2              colorspace_1.4-1            GetoptLong_0.1.7            tidyselect_0.2.5           
[16] compiler_3.6.0              preprocessCore_1.46.0       fdrtool_1.2.15              DelayedArray_0.10.0         scales_1.0.0               
[21] readr_1.3.1                 affy_1.62.0                 digest_0.6.19               XVector_0.24.0              pkgconfig_2.0.2            
[26] htmltools_0.3.6             limma_3.40.2                htmlwidgets_1.3             rlang_0.3.4                 GlobalOptions_0.1.0        
[31] rstudioapi_0.10             shiny_1.3.2                 shape_1.4.4                 zoo_1.8-5                   mzID_1.22.0                
[36] BiocParallel_1.17.18        dplyr_0.8.0.1               RCurl_1.95-4.12             magrittr_1.5                GenomeInfoDbData_1.2.1     
[41] MALDIquant_1.19.3           Rcpp_1.0.1                  munsell_0.5.0               S4Vectors_0.22.0            vsn_3.52.0                 
[46] MASS_7.3-51.4               SummarizedExperiment_1.14.0 zlibbioc_1.30.0             plyr_1.8.4                  grid_3.6.0                 
[51] promises_1.0.1              shinydashboard_0.7.1        crayon_1.3.4                lattice_0.20-38             circlize_0.4.6             
[56] hms_0.4.2                   mzR_2.18.0                  ComplexHeatmap_2.0.0        pillar_1.4.0                GenomicRanges_1.36.0       
[61] rjson_0.2.20                codetools_0.2-16            XML_3.98-1.19               glue_1.3.1                  BiocManager_1.30.4         
[66] png_0.1-7                   httpuv_1.5.1                foreach_1.4.4               tidyr_0.8.3                 gtable_0.3.0               
[71] purrr_0.3.2                 clue_0.3-57                 assertthat_0.2.1            ggplot2_3.1.1               mime_0.6                   
[76] xtable_1.8-4                later_0.8.0                 ncdf4_1.16.1                tibble_2.1.1                iterators_1.0.10           
[81] IRanges_2.18.0              cluster_2.0.9
Entering edit mode
arne.smits ▴ 70
Last seen 2.1 years ago

Dear robert.arsie,

A fix is available in the new version of DEP (version 1.6.1).

# Update Bioconductor packages:

Best, Arne


Login before adding your answer.

Traffic: 208 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6