Entering edit mode
andreas.pittroff
•
0
@andreaspittroff-12269
Last seen 7.7 years ago
Hello all,
I am currently trying to merge two expressionsets using the merge-function of InSilicoMerging.
Here is the code I was using (error translates to "error in data frame [...] arguments imply differing number of rows: 181, 107":
mergeset = merge(esets) Fehler in data.frame(<S4 object of class "ExpressionSet">, <S4 object of class "ExpressionSet">, : Argumente implizieren unterschiedliche Anzahl Zeilen: 181, 107
Traceback: > traceback() 10: stop(gettextf("arguments imply differing number of rows: %s", paste(unique(nrows), collapse = ", ")), domain = NA) 9: data.frame(<S4 object of class "ExpressionSet">, <S4 object of class "ExpressionSet">, check.names = TRUE, stringsAsFactors = TRUE) 8: eval(expr, envir, enclos) 7: eval(as.call(c(expression(data.frame), x, check.names = !optional, stringsAsFactors = stringsAsFactors))) 6: as.data.frame.list(x) 5: as.data.frame(x) 4: merge(as.data.frame(x), as.data.frame(y), ...) 3: merge.default(esets) 2: merge(esets) 1: merge(esets)
The two expressionsets look like this:
> eset1b ExpressionSet (storageMode: lockedEnvironment) assayData: 54675 features, 181 samples element names: exprs protocolData: none phenoData sampleNames: GSM1213669 GSM1213670 ... GSM1213849 (181 total) varLabels: title geo_accession ... data_row_count (35 total) varMetadata: labelDescription featureData featureNames: 1007_s_at 1053_at ... AFFX-TrpnX-M_at (54675 total) fvarLabels: ID GB_ACC ... Gene Ontology Molecular Function (16 total) fvarMetadata: Column Description labelDescription experimentData: use 'experimentData(object)' Annotation: GPL570 > eset3b ExpressionSet (storageMode: lockedEnvironment) assayData: 22283 features, 107 samples element names: exprs protocolData: none phenoData sampleNames: GSM254625 GSM254626 ... GSM254731 (107 total) varLabels: title geo_accession ... data_row_count (35 total) varMetadata: labelDescription featureData featureNames: 1007_s_at 1053_at ... AFFX-TrpnX-M_at (22283 total) fvarLabels: ID GB_ACC ... Gene Ontology Molecular Function (16 total) fvarMetadata: Column Description labelDescription experimentData: use 'experimentData(object)' Annotation: GPL96
Merge did definitely run before on my system using this line, when I did load the expressionsets using InsilicoDB. InSilicoDB is currently not working however, so I am using GEOquery to import the datasets I downloaded with InSilicoDB before (maybe there is a problem here?)
Regards
The function merge() from insilicodb is obscured by basic R merge, try to run inSilicoMerging::merge(esets, method="NONE")