data(prostateIds)
expandedWelsh <- expandedProfile(welsh01EntrezIDs[1:100], onto="ANY",
level=2, orgPackage="org.Hs.eg.db")
expandedSingh <- expandedProfile(singh01EntrezIDs[1:100], onto="ANY",
level=2, orgPackage="org.Hs.eg.db")
commonGenes <- intersect(welsh01EntrezIDs[1:100], singh01EntrezIDs[1:100])
commonExpanded <- expandedProfile(commonGenes, onto="ANY", level=2,
orgPackage="org.Hs.eg.db")
### FUnciona si fem:
equivMF <-equivalentGOProfiles (expandedWelsh[["MF"]],
qm = expandedSingh[["MF"]],
pqn0= commonExpanded[["MF"]])
equivsList <- lapply(1:length(expandedWelsh),
function (onto){
equivalentGOProfiles (expandedWelsh[[onto]],
qm = expandedSingh[[onto]],
pqn0= commonExpanded[[onto]])
}
)
Run the code above in your browser using DataLab