[Bioc-devel] Duplicate issue in ELMER get.pair for HNSC data

ghafarpour gh@f@rpour @ending from ut@@c@ir
Tue Nov 20 13:57:52 CET 2018


Hi, 

When I try to analyze HNSC data by ELMER using this pipeline 

TCGA.pipe(disease = "HNSC",analysis = "all", genome = "hg19", cores = 1,
permu.size=10000, Pe=0.001) 

I got this error

Error in `.rowNamesDF<-`(x, value = value) : 
  duplicate 'row.names' are not allowed
In addition: Warning message:
non-unique values when setting 'row.names':
'cg00042263.ENSG00000146556', 'cg00045114.ENSG00000234741',
'cg00055603.ENSG00000162365', 'cg00081439.ENSG00000173876',
'cg00082912.ENSG00000011052', 'cg00082912.ENSG00000243678',
'cg00180909.ENSG00000166984', 'cg00216693.ENSG00000075886',
'cg00216693.ENSG00000173272', 'cg00253248.ENSG00000011052',
'cg00253248.ENSG00000243678', 'cg00263146.ENSG00000183666',
'cg00288376.ENSG00000197978', 'cg00288376.ENSG00000225273',
'cg00379619.ENSG00000105640', 'cg00379713.ENSG00000165478',
'cg00411945.ENSG00000172661', 'cg00432262.ENSG00000136270',
'cg00452252.ENSG00000243716', 'cg00459445.ENSG00000146574',
'cg00460916.ENSG00000262304', 'cg00463746.ENSG00000116678',
'cg00482162.ENSG00000015568', 'cg00522451.ENSG00000169629',
'cg00570580.ENSG00000183054', 'cg00610294.ENSG00000258755',
'cg00611227.ENSG00000146556', 'cg00634968.ENSG00000244067',
'cg00655552.ENSG00000151718', 'cg00664723.ENSG00000183666',
'cg00686132.ENSG00000166984', 'cg00750880.ENSG0000023491 [... truncated]


It is the result of running : 

pair <- get.pair(data = mae.hg19,
                 group.col = group.col,
                 group1 = group1,
                 mode = "unsupervised",
                 group2 = group2,
                 nearGenes = nearGenes,
                 diff.dir = diff.dir,
                 minSubgroupFrac = 0.4, # % of samples to use in to
create groups U/M
                 permu.dir = file.path(dir.out,"permu"),
                 permu.size = 100000, # Please set to 100000 to get
significant results
                 raw.pvalue = 0.05,   
                 Pe = 0.001, # Please set to 0.001 to get significant
results
                 filter.probes = TRUE, # See
preAssociationProbeFiltering function
                 filter.percentage = 0.05,
                 filter.portion = 0.3,
                 dir.out = dir.out,
                 cores = 1,
                 label = diff.dir) 

would you please help me to solve it 

https://support.bioconductor.org/p/115251/
	[[alternative HTML version deleted]]



More information about the Bioc-devel mailing list