Biostar Beta. Not for public use.
Question: GOplot and making circle_dat out of microarray and DAVID data
Entering edit mode

I am going to use GOplot package of R to visualize my ontology results out of microarray analysis of cancerous and normal tissues. I have problem with using this code:

circ <- circle_dat(EC$david, EC$genelist)

I have two text files. One for microarray with columns namely: ID, adj.P.Val, logFC, AveExpr, B, and t. Another file is for david data with columns of: Category, ID, Term, PValue, and Genes. I import these two files into R-studio with these codes:

genelist <- read.csv("genelist.txt", sep = "\t", header = TRUE)
david <- read.csv("BP1.txt", sep = "\t", header = TRUE)

But when I am going to build the circle_dat out of these with this code:

circ <- circle_dat(david, genelist)

I get the error of:

Error in data.frame(category = rep(as.character(terms$category), count),  : 
  arguments imply differing number of rows: 8765, 0

How should I form the circle_dat for my data. By the way, I am just a biologist and does not have any strong background in programming, so please kindly give me a code for this because I cannot understand the programming explanations. I have already saw this link and could not understand anything :( Thank you.

ADD COMMENTlink 12 months ago farzaneafzali • 50 • updated 12 months ago MasMarius • 10
Entering edit mode

Looks like your data has 'Category' whereas the function looks for 'category' (small 'c').

You should also look into your objects to see how they look. It is a good habit to have. For example:


As a side note, I posted a tutorial for plotting DAVID results: Clustering of DAVID gene enrichment results from gene expression studies

ADD REPLYlink 12 months ago
Kevin Blighe
Entering edit mode

Thanks a lot Kevin. I've already figured out how to plot my data via GOplot and now I am going to make a tutorial for Youtube. I will post the link later here so everyone can use it especially biologists. By the way, I loved your way of plotting, thank you for introducing that to me.

ADD REPLYlink 12 months ago
• 50
Entering edit mode

@farzaneafzali could you also share it by written? It would be much appreciated. Thanks in advance!

ADD REPLYlink 12 months ago
• 10

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.0