Enter the body of text here
Code should be placed in three backticks as shown below
peakonly_DMSO2 <- readPeakFile("path_to_peak_bed_DMSO")
geneonly_DMSO <- seq2gene(peakonly_DMSO2, tssRegion = c(-1000, 1000), flankDistance = 3000, TxDb=txdb)
pathway2only_DMSO <- enrichPathway(geneonly_DMSO)
dotplot(pathway2only_DMSO)
i get dotplot graph and in the first annotation i get that there is 130 gene that involve, and i want to know what is this gene, how can i extract it? if i write:
pathway2only_DMSO[1]$geneID
gene_list <- strsplit(pathway2only_DMSO[1]$geneID, "/")[[1]]
i get:
this list:
gene_list
[1] "2475" "842" "3065" "5525" "80351" "1937" "572" "84447" "4041" "9846" "8322"
[12] "867" "7450" "4040" "3845" "8496" "2065" "9612" "2100" "207" "83464" "5604"
[23] "8312" "7249" "930" "5432" "9611" "6774" "5296" "7311" "5714" "6868" "10018"
[34] "182" "6238" "7529" "22916" "2261" "55109" "9542" "5159" "57670" "11160" "11064"
[45] "79109" "5293" "1185" "3932" "5906" "89780" "818" "10718" "8036" "6934" "5682"
[56] "408" "2120" "51411" "5908" "4254" "283455" "1024" "8660" "5529" "5663" "9321"
[67] "4140" "4088" "8826" "23164" "399687" "4763" "5717" "6777" "10197" "10014" "3690"
[78] "5573" "2885" "57674" "4089" "5605" "208" "238" "6654" "27436" "27248" "7039"
[89] "6772" "2066" "598" "5594" "613" "9978" "10533" "1499" "22906" "79718" "5290"
[100] "152831" "5156" "27123" "2247" "55534" "6500" "8841" "2246" "815" "10817" "5528"
[111] "2309" "2534" "2099" "9734" "1956" "1523" "27044" "8805" "80853" "673" "8454"
[122] "8658" "3084" "23189" "6907" "22866" "55869" "5519" "7314" "471"
so how can i know the full name gene from this R names?
thank you very nutch for any answer!!!!
thanks!!!!!!!!!!!!!!!