# # # R代码从装饰图案的useHomology来源。Rnw“# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #代码块1号:loadAndListPackages # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #库(DBI)库(hom.Hs.inp.db) ls(包:hom.Hs.inp.db) # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #代码块2号:listMapping # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # as.list (hom.Hs.inpMUSMU[1:4]) # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #代码块3号:mapGeneExample # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #加载生物注释数据为人类库(org.Hs.eg.db) # entrex基因ID和运用蛋白质ID基因“MSX2”选择(org.Hs.eg象征。db,键=“MSX2”,列= c (“ENTREZID”、“ENSEMBLPROT”), keytype =“象征”)#使用inparanoid包让老鼠基因被认为是#相当于“ENSP00000239243”选择(hom.Hs.inp运用蛋白质ID。db,键=“ENSP00000239243”,列=“MUS_MUSCULUS keytype =“HOMO_SAPIENS”) #加载生物注释数据鼠标图书馆(org.Mm.eg.db) # entrez基因ID和“ENSMUSP00000021922”选择(org.Mm.eg的象征。db,键=“ENSMUSP00000021922”,列= c (“ENTREZID”、“符号”),keytype =“ENSEMBLPROT”) # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #代码块数量4:seedPairExample # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # mget (“ENSP00000301011”, hom.Hs.inpMUSMU) # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #代码块5号:seedPairExample2 # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #人类连接数据库mycon <——hom.Hs.inp_dbconn() #列出所有可用的表的db头(dbListTables (mycon)) #列出表中的列的感兴趣的dbListFields (mycon, mus_musculus) # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #代码块6号:seedPairExample3 # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #做一个查询,让我们看看哪些clust_id我们需要sql < -“SELECT *从mus_musculus inp_id =“ENSP00000301011”;“#retrieve the data dataOut <- dbGetQuery(mycon, sql) dataOut ################################################### ### code chunk number 7: seedPairExample4 ################################################### #make a query that will let us see all the data that is affiliated with a clust id sql <- "SELECT * FROM mus_musculus WHERE clust_id = '1731';" #retrieve the data dataOut <- dbGetQuery(mycon, sql) dataOut ################################################### ### code chunk number 8: useHomology.Rnw:250-251 ################################################### toLatex(sessionInfo())