我正在使用R软件包“ anrichment”进行细胞富集分析。但是,当我使用richmentAnalysis()函数时发生错误。
下面是错误消息。
Error in base::mapply(..., SIMPLIFY = SIMPLIFY) :
zero-length inputs cannot be mixed with those of non-zero length
请告诉我如何解决此错误? 追溯消息中的“长度”可能是导致此错误的原因,但我不知道该“长度”的来源。
9.
.mymapply(rep, 1:n, lengths)
回溯消息是这个。
10.
base::mapply(..., SIMPLIFY = SIMPLIFY)
9.
.mymapply(rep, 1:n, lengths)
8.
do.call("c", .mymapply(rep, 1:n, lengths))
7.
.indexedFlattenedList(lapply(groupList, getElement, "alternateNames"))
6.
impliedGroups(collection$groups, queryGroups = NULL, includeSelf = TRUE,
get = "parents")
5.
.matchDataSets(collection, tags = tags, matchComponents = matchComponents,
searchType = searchType, invert = invertSearch, exactMatch = exactMatch,
fixed = fixed, firstDate = firstDate, lastDate = lastDate,
dateFormat = dateFormat, invert.date = invertDateSearch)
4.
.geneLists(collection, evidence = evidence, tags = tags, matchComponents = matchComponents,
searchType = searchType, invertSearch = invertSearch, exactMatch = exactMatch,
fixed = fixed, firstDate = firstDate, lastDate = lastDate,
dateFormat = dateFormat, invertDateSearch = invertDateSearch)
3.
allGeneSetGenes(refCollection, evidence = useEvidence)
2.
enrichmentAnalysis(classLabels = Group, identifiers = entrez,
refCollection = collection, useBackground = "given", threshold = 1,
thresholdType = "Bonferroni", getOverlapEntrez = FALSE, getOverlapSymbols = TRUE,
maxReportedOverlapGenes = 10000, removeDuplicatesInDifferentClasses = FALSE, ...
1.
我的代码在下面。
library(anRichment)
analysisresult = enrichmentAnalysis(
classLabels = Group,
identifiers = entrez,
refCollection = collection,
useBackground = "given",threshold = 1,
thresholdType = "Bonferroni",
getOverlapEntrez = FALSE,
getOverlapSymbols= TRUE,
maxReportedOverlapGenes = 10000,
removeDuplicatesInDifferentClasses=FALSE,
entrySeparator = ",",
ignoreLabels = "grey",
combineEnrichmentTables = FALSE)