if(file.exists('SRAmetadb.sqlite')) {
library(SRAdb)
sra_dbname <- 'SRAmetadb.sqlite'
sra_con <- dbConnect(dbDriver("SQLite"), sra_dbname)
## create a graphNEL object from SRA accessions, which are full text search results of terms 'primary thyroid cell line'
acc <- getSRA (search_terms ='primary thyroid cell line', out_types=c('sra'), sra_con=sra_con, acc_only=TRUE)
g <- entityGraph(acc)
## Not run:
# library(Rgraphviz)
# attrs <- getDefaultAttrs(list(node=list(fillcolor='lightblue', shape='ellipse')))
# plot(g, attrs= attrs)
# ## End(Not run)
} else {
print("use getSRAdbFile() to get a copy of the SRAmetadb.sqlite file
and then rerun the example")
}
Run the code above in your browser using DataLab