Open herong2006 opened 1 year ago
GSE50244.bulk.eset = readRDS('GSE50244bulkeset.rds') GSE50244.bulk.eset bulk.mtx = exprs(GSE50244.bulk.eset) bulk.meta = exprs(GSE50244.bulk.eset)
EMTAB.sce = readRDS('EMTABsce_healthy.rds') EMTAB.sce
XinT2D.sce = readRDS('XinT2Dsce.rds') XinT2D.sce
bulk.mtx = exprs(GSE50244.bulk.eset)
Est.prop.GSE50244 = music_prop(bulk.mtx = bulk.mtx, sc.eset = EMTAB.sce, clusters = 'cellType', samples = 'sampleID', select.ct = c('alpha', 'beta', 'delta', 'gamma', 'acinar', 'ductal'), verbose = F) run with error:
Est.prop.GSE50244 = music_prop(bulk.mtx = bulk.mtx, sc.eset = EMTAB.sce, clusters = 'cellType', samples = 'sampleID', select.ct = c('alpha', 'beta', 'delta', 'gamma', 'acinar', 'ductal'), verbose = F) _Error in musicbasis(sc.sce, non.zero = TRUE, markers = sc.markers, clusters = clusters, : argument "sc.sce" is missing, with no default
Est.prop.GSE50244 = music_prop(bulk.mtx = bulk.mtx, sc.eset = EMTAB.sce, clusters = 'cellType',
I have also encountered this problem as well and created a separate issue #111. Would appreciate any suggestions if you ever found a solution.
GSE50244.bulk.eset = readRDS('GSE50244bulkeset.rds') GSE50244.bulk.eset bulk.mtx = exprs(GSE50244.bulk.eset) bulk.meta = exprs(GSE50244.bulk.eset)
EMTAB.sce = readRDS('EMTABsce_healthy.rds') EMTAB.sce
XinT2D.sce = readRDS('XinT2Dsce.rds') XinT2D.sce
Bulk expression matrix
bulk.mtx = exprs(GSE50244.bulk.eset)
Estimate cell type proportions
Est.prop.GSE50244 = music_prop(bulk.mtx = bulk.mtx, sc.eset = EMTAB.sce, clusters = 'cellType', samples = 'sampleID', select.ct = c('alpha', 'beta', 'delta', 'gamma', 'acinar', 'ductal'), verbose = F) run with error: