Section 8 Write Tables
8.1 Table S1 - Paralog Pairs
library(openxlsx)
#Reload Protein Complexes
DrerioComplex <- ComplexSet(organism = "drerio",identifier = "ENSEMBL",grepOnCol = F)
HsapieComplex <- ComplexSet(organism = "hsapiens",identifier = "SYMBOL",grepOnCol = F)
#Attach info on Complexes on each dataset
AllCorrelations$complex <- GetPairComplex(AllCorrelations[,c('p.1','p.2')],
DrerioComplex)
AllCorrelations.tr$complex <- GetPairComplex(AllCorrelations.tr[,c('p.1','p.2')],
HsapieComplex)
AllCorrelations.pr$complex <- GetPairComplex(AllCorrelations.pr[,c('p.1','p.2')],
HsapieComplex)
#Create a list
l <- list('Embryo Development'= RNASeqData.exprs,
'Tissues - Transcriptome' = em_transcr.exprs,
'Tissues - Proteome' = em_proteome.exprs,
'Paralog pairs (Embryo)' = AllCorrelations,
'Paralog pairs (Transc.)' = AllCorrelations.tr,
'Paralog pairs (Prot.)' = AllCorrelations.pr)
#Write table
openxlsx::write.xlsx(l,'../out/tables/Table_S1.xlsx')8.2 Table S2 - Protein Complexes
#Create a list
l <- list('Embryo Development'= Complexes_Table_Drerio,
'Tissues - Transcriptome' = Complexes_Table_Hsap_T,
'Tissues - Proteome' = Complexes_Table_Hsap_P)
#Write table
openxlsx::write.xlsx(l,'../out/tables/Table_S2.xlsx')8.3 Table S3 - GOEnrichment
#Combined GOPlot
tab1 <- GOEnrich.var.Pairs.dev$table;tab1$cond <- "Embryo Development"
tab2 <- GOEnrich.var.Pairs.tr$table;tab2$cond <- "Tissues - Transcriptome"
tab3 <- GOEnrich.var.Pairs.pr$table;tab3$cond <- "Tissues - Proteome"
GOTable_tab <- rbind(tab1,tab2,tab3)
l <- list('GOEnrichment '= GOTable_tab)
#Write table
openxlsx::write.xlsx(l,'../out/tables/Table_S3.xlsx')