Question: including information indicating upregulation/downregulation of genes in goana()
gravatar for jhj89
23 months ago by
jhj8910 wrote:


To put it simply, I have a list of genes with information about whether they are upregulated or downregulated. What should the format of the input be when using goana() such that I can include the information about upregulation or downregulation of each of the genes? I understand that goana() can accept a vector of genes or MArrayLM fit object.

Thank you.

ADD COMMENTlink modified 23 months ago by James W. MacDonald47k • written 23 months ago by jhj8910
gravatar for James W. MacDonald
23 months ago by
United States
James W. MacDonald47k wrote:

It's not necessarily clear from the help for goana, but there's this:


      de: a vector of Entrez Gene IDs, or a list of such vectors, or an
          'MArrayLM' fit object.

And the trick is the 'a list of such vectors' part, which is what is extracted if you feed goana an MArrayLM object. So if you do

deGenes <- list(Up = <vector of up-regulated Entrez Gene IDs>, Down = <vector of down-regulated Entrez Gene IDs>)

result <- goana(deGenes, universe = <universe of Entrez Gene IDs from your experiment>, <other arguments go here>)

Then you will get tests for up and down regulated genes.

ADD COMMENTlink written 23 months ago by James W. MacDonald47k

The example code on the help page for goana might also help. See the example starting ## Default usage.

ADD REPLYlink modified 23 months ago • written 23 months ago by Gordon Smyth34k

Dear James and Gordon,

Thank you so much! I got it to work.

ADD REPLYlink written 23 months ago by jhj8910
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 123 users visited in the last hour