Question: DEP package: error in plot_heatmap
0
gravatar for robert.arsie
14 days ago by
robert.arsie0 wrote:

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:

library(DEP)
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.data.frame(df) is not TRUE

Can I do something to solve this issue?

dep • 69 views
ADD COMMENTlink modified 9 days ago by arne.smits60 • written 14 days ago by robert.arsie0
1

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

ADD REPLYlink written 13 days ago by arne.smits60
> 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

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] 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
ADD REPLYlink written 14 days ago by robert.arsie0
Answer: DEP package: error in plot_heatmap
0
gravatar for arne.smits
9 days ago by
arne.smits60
arne.smits60 wrote:

Dear robert.arsie,

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

# Update Bioconductor packages:
BiocManager::install()

Best, Arne

ADD COMMENTlink written 9 days ago by arne.smits60
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 302 users visited in the last hour