R/Filter33Motifs.R

In jakeyeung/TissueCiradianAnalysis: Analysis of oscillatory RNASeq data

Filter33Motifs <- function(pres_motifs){
  motif_of_int_clean=c("V_E4BP4_01",
                       "V_RORA2_01",
                       "V_CREBATF_Q6",
                       "BHLHE40",
                       "V_AR_Q6_01",
                       "V_FOXO1_Q5",
                       "SP1",
                       "V_ZF5_01",
                       "V_PPARG_01",
                       "V_AP2_Q3",
                       "V_ZIC2_04",
                       "V_PU1_Q4",
                       "V_RXRA_03",
                       "V_ERR1_Q2_01",
                       "V_MITF_Q6",
                       "HNF1A_M161",
                       "V_NFY_01",
                       "V_SREBP1_02",
                       "V_GABPAGABPB_Q6",
                       "Jdp2_M131",
                       "V_SMAD_Q6",
                       "V_HNF4A_02",
                       "V_NKX62_Q2",
                       "V_EKLF_Q5",
                       "V_GATA3_05",
                       "V_OCT1_01",
                       "V_LXRDR4_Q3",
                       "V_MAF_Q6_01",
                       "HSF1_M61",
                       "V_MEF2A_Q6",
                       "V_FXR_Q2",
                       "V_MTF1_02",
                       "EGR1_M3")
  x33family <- pres_motifs[, match(motif_of_int_clean ,colnames(pres_motifs))]
  
  colnames(x33family)<-c("D_Box",
                         "RRE",
                         "CREB",
                         "E_Box",
                         "GR",
                         "FOX",
                         "SP1", 
                         "ZF5",
                         "PPAR",  
                         "AP2", 
                         "ZIC2", 
                         "PU1", 
                         "RXR", 
                         "ERE", 
                         "MITF", 
                         "HNF1", 
                         "NFY", 
                         "SREB", 
                         "GABP", 
                         "JDP2", 
                         "SMAD",
                         "HNF4",
                         "NKX6", 
                         "KLF", 
                         "GATA3", 
                         "OCT1", 
                         "LXR",
                         "MAF", 
                         "HSF", 
                         "MEF2A", 
                         "FXR", 
                         "MTF1", 
                         "EGR1")  
  return(x33family)
}