I get an error running Celaref's "contrasteachgrouptothe_rest()" (details below). Traceback() suggests the problem is related to how data is passed to MAST. This is my first time using celaref, so I guess my inputs are not correctly configured. But I can't see the issue. Can anyone offer any insights?
Thanks!
Hamid
--Details:
cellInfoTbl <- data.frame(cellId=names(clusterIDs), Sample="blank" , Cluster=clusterIDs, group="group")
datasetse <- loadsefromtables(countsmatrix = normalizedPBMC, + cellinfotable = cellInfoTbl, + groupcol_name = "group")
datasetse <- trimsmallgroupsandlowexpressiongenes(datasetse, + minlibsize = 100, # set low to keep all cells! + mingroupmembership = 5,
+ mindetectedbyminsamples = 5)detablequery <- contrasteachgrouptotherest(datasetse, "PBMC3K", num_cores = 4)
Error in checkArrayNames(exprsArray, cData, fData) :
exprsArray must be numeric
traceback()
7: stop("exprsArray must be numeric")
6: checkArrayNames(exprsArray, cData, fData)
5: MAST::FromMatrix(loggedcounts, cData = coldataforMAST, fData = rowData(dataset_se))
4: FUN(X[[i]], ...)
3: lapply(X = X, FUN = FUN, ...)
2: parallel::mclapply(groups2test, FUN = contrastthegrouptotherest, datasetse = datasetse, mc.cores = numcores)
1: contrasteachgrouptotherest(datasetse, "PBMC3K", num_cores = 4) >
str(clusterIDs)
Factor w/ 10 levels "0","1","2","3",..: 6 4 2 5 7 2 6 6 6 5 ...
- attr(*, "names")= chr [1:2638] "AAACATACAACCAC" "AAACATTGAGCTAC" "AAACATTGATCAGC" "AAACCGTGCTTCCG" ...
class(normalizedPBMC)
[1] "dgCMatrix" attr(,"package") [1] "Matrix" >
assay(dataset_se)[1:4,1:4]
4 x 4 sparse Matrix of class "dgCMatrix"
AAACATACAACCAC AAACATTGAGCTAC AAACATTGATCAGC AAACCGTGCTTCCG
NOC2L . . . .
HES4 . . . .
ISG15 . . 0.6931472 2.197225
TNFRSF18 . 0.6931472 . .
>
sessionInfo() R version 3.5.2 (2018-12-20) Platform: x86_64-apple-darwin15.6.0 (64-bit) Running under: macOS Mojave 10.14.4
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] enUS.UTF-8/enUS.UTF-8/enUS.UTF-8/C/enUS.UTF-8/en_US.UTF-8
attached base packages: [1] compiler parallel stats4 stats graphics grDevices utils datasets methods base
other attached packages:
[1] celarefData1.0.0 ExperimentHub1.8.0 AnnotationHub2.14.2 celaref1.0.1 SummarizedExperiment1.12.0 DelayedArray0.8.0
[7] BiocParallel1.16.5 matrixStats0.54.0 GenomicRanges1.34.0 GenomeInfoDb1.18.1 lsa0.73.1 SnowballC0.6.0
[13] Hmisc4.2-0 Formula1.2-3 survival2.43-3 lattice0.20-38 CellMix1.6.2 GSEABase1.44.0
[19] graph1.60.0 annotate1.60.0 XML3.98-1.16 stringr1.4.0 csSAM1.2.4 NMF0.21.0
[25] cluster2.0.7-1 rngtools1.3.1 pkgmaker0.27 registry0.5 pathview1.22.3 org.Hs.eg.db3.7.0
[31] AnnotationDbi1.44.0 IRanges2.16.0 S4Vectors0.20.1 Biobase2.42.0 BiocGenerics0.28.0 gage2.32.1
[37] clusterProfiler3.10.1 ggplot23.1.0 dplyr0.8.0.1 Seurat3.0.0
loaded via a namespace (and not attached):
[1] R.methodsS31.7.1 tidyr0.8.2 acepack1.4.1 bit640.9-7 knitr1.21 irlba2.3.3
[7] R.utils2.7.0 data.table1.12.0 rpart4.1-13 KEGGREST1.22.0 RCurl1.95-4.11 doParallel1.0.14
[13] metap1.0 preprocessCore1.44.0 callr3.1.1 cowplot0.9.4 usethis1.4.0 RSQLite2.1.1
[19] RANN2.6.1 europepmc0.3 future1.11.0 bit1.1-14 enrichplot1.2.0 httpuv1.4.5.1
[25] xml21.2.0 assertthat0.2.0 viridis0.5.1 xfun0.4 hms0.4.2 promises1.0.1
[31] evaluate0.13 BiocInstaller1.32.1 fansi0.4.0 progress1.2.0 caTools1.17.1.1 Rgraphviz2.26.0
[37] igraph1.2.2 DBI1.0.0 htmlwidgets1.3 purrr0.3.0 backports1.1.3 gbRd0.4-11
[43] gridBase0.4-7 SingleCellExperiment1.4.1 remotes2.0.2 ROCR1.0-7 abind1.4-5 withr2.1.2
[49] ggforce0.2.2 triebeard0.3.0 checkmate1.9.1 sctransform0.2.0 prettyunits1.0.2 DOSE3.8.2
[55] ape5.3 lazyeval0.2.1 crayon1.3.4 genefilter1.64.0 pkgconfig2.0.2 labeling0.3
[61] tweenr1.0.1 nlme3.1-137 pkgload1.0.2 nnet7.3-12 devtools2.0.2 rlang0.3.1
[67] globals0.12.4 rsvd1.0.0 rprojroot1.3-2 polyclip1.10-0 lmtest0.9-36 Matrix1.2-15
[73] urltools1.7.3 zoo1.8-4 base64enc0.1-3 beeswarm0.2.3 ggridges0.5.1 processx3.2.1
[79] png0.1-7 viridisLite0.3.0 bitops1.0-6 R.oo1.22.0 KernSmooth2.23-15 Biostrings2.50.2
[85] blob1.1.1 qvalue2.14.1 gridGraphics0.3-0 scales1.0.0 lpSolve5.6.13 memoise1.1.0
[91] magrittr1.5 plyr1.8.4 ica1.0-2 gplots3.0.1 bibtex0.4.2 gdata2.18.0
[97] zlibbioc1.28.0 lsei1.2-0 RColorBrewer1.1-2 KEGGgraph1.42.0 fitdistrplus1.0-14 cli1.0.1
[103] XVector0.22.0 listenv0.7.0 pbapply1.3-4 ps1.3.0 htmlTable1.13.1 MASS7.3-51.1
[109] limSolve1.5.5.3 tidyselect0.2.5 MAST1.8.2 stringi1.3.1 yaml2.2.0 GOSemSim2.8.0
[115] latticeExtra0.6-28 ggrepel0.8.0 grid3.5.2 fastmatch1.1-0 tools3.5.2 future.apply1.1.0
[121] rstudioapi0.9.0 foreach1.4.4 foreign0.8-71 gridExtra2.3 farver1.1.0 Rtsne0.15
[127] ggraph1.0.2 digest0.6.18 rvcheck0.1.3 BiocManager1.30.4 shiny1.2.0 quadprog1.5-5
[133] Rcpp1.0.0 SDMTools1.1-221 later0.7.5 httr1.4.0 npsurv0.4-0 Rdpack0.10-1
[139] colorspace1.4-0 fs1.2.6 reticulate1.10 splines3.5.2 ggplotify0.0.3 plotly4.8.0
[145] sessioninfo1.1.1 xtable1.8-3 jsonlite1.6 corpcor1.6.9 UpSetR1.3.3 R62.4.0
[151] mime0.6 pillar1.3.1 htmltools0.3.6 glue1.3.0 interactiveDisplayBase1.20.0 codetools0.2-16
[157] fgsea1.8.0 utf81.1.4 pkgbuild1.0.2 tsne0.1-3 tibble2.0.1 curl3.3
[163] gtools3.8.1 GO.db3.7.0 rmarkdown1.11 desc1.2.0 munsell0.5.0 GenomeInfoDbData1.2.0
[169] DO.db2.9 iterators1.0.10 reshape21.4.3 gtable0.2.0
see below, should be answer.