This is a test version of Biostars. For the public version, visit https://www.biostars.org.
combat-seq: group variable

I'm using combat-seq to perform a metanalysis, but I don't understand very well what the group variable means... could you help me out with an example? what do they call a biological variable? are the biological variables the controls?

In summary, how can I specify to combat seq in a command line that my first samples are healthy controls, followed by 3 different batches of cancer samples?

combat-seq

If you haven’t, I recommend to read this (and perhaps the preceding and following sections).

ok, thank you very much, it worked, but it is giving me a list.. where do I get the raw counts table?

1 answer

It has a argument batch to define categorical batches to remove. It has an argument group for biological variables of interest to preserve. Alternatively, you can use covar_mod to provide a design matrix (model.matrix()) containing your experimental design WITHOUT the batch variable. The latter two arguments are used to preserve the biological effects. That obviously only works if batch is not nested with any of it.

As always when I read "meta-analysis", I say that it's not a meta-analysis to combine data into a single quantitative analysis, as they're not independent anymore. Meta would be to analyze separately and aggregate results, like Fisher, Stouffer, RRA, intersect, whatever.

Log in to answer this question.