statOmics/tradeSeq

Combining multiple gams together

Opened this issue · 0 comments

Hi,

I have multiple genotypes that I ran gams for and I want to combine all the gams together before clustering so I can have my clusters comparable with all the genotypes.

source("path/gam2.R")
wt_gam2 <- strain_fit(pt, "WT")
dmd3_gam2 <- strain_fit(pt, "dmd-3(-)")
mab3_gam2 <- strain_fit(pt, "mab-3(-)")
double_gam2 <- strain_fit(pt, "mab-3(-);dmd-3(-)")
wt_herm_gam2 <- strain_fit(pt, "WT_herm")

all <- rbind(wt_gam2,dmd3_gam2,mab3_gam2, double_gam2, wt_herm_gam2)
saveRDS(all, "/output/gam2s/all_exp.rds")

The error I keep getting is: Error in .rbind.SummarizedExperiment(args) : '...' objects must have the same colnames

I tried other ways of combining the gams (such as cbind, etc) but nothing works.. Any suggestions?

Here is the error for cbind:
Error in value[3L] :
failed to rbind column 'tradeSeq' across DataFrame objects: duplicate 'row.names' are not allowed
In addition: Warning message:
non-unique values when setting 'row.names': ‘1’, ‘10’, ‘100’, ‘101’, ‘102’, ‘103’, ‘104’, ‘105’, ‘106’, ‘107’, ‘108’, ‘109’, ‘11’, ‘110’, ‘111’, ‘112’, ‘113’, ‘114’, ‘115’, ‘116’, ‘117’, ‘118’, ‘119’, ‘12’, ‘120’, ‘121’, ‘122’, ‘123’, ‘124’, ‘125’, ‘126’, ‘127’, ‘128’, ‘129’, ‘13’, ‘130’, ‘131’, ‘132’, ‘133’, ‘134’, ‘135’, ‘136’, ‘137’, ‘138’, ‘139’, ‘14’, ‘140’, ‘141’, ‘142’, ‘143’, ‘144’, ‘145’, ‘146’, ‘147’, ‘148’, ‘149’, ‘15’, ‘150’, ‘151’, ‘152’, ‘153’, ‘154’, ‘155’, ‘156’, ‘157’, ‘158’, ‘159’, ‘16’, ‘160’, ‘161’, ‘162’, ‘163’, ‘164’, ‘165’, ‘166’, ‘167’, ‘168’, ‘169’, ‘17’, ‘170’, ‘171’, ‘172’, ‘173’, ‘174’, ‘175’, ‘176’, ‘177’, ‘178’ [... truncated]