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