data(Selection)
data(infos)
#Find probesets for EGR1
pbst_EGR1 = infos[infos$hgnc_symbol=="EGR1", "affy_hg_u133_plus_2"]
gene_IDs = infos[match(Selection@name, infos$affy_hg_u133_plus_), "hgnc_symbol"]
data(network)
#A nv value can chosen using the cutoff function
nv=.11
EGR1<-which(is.element(Selection@name,pbst_EGR1))
P<-position(network,nv=nv)
geneNeighborhood(network,targets=EGR1,nv=nv,ini=P,
label_v=gene_IDs)
Run the code above in your browser using DataLab