This is a test version of Biostars. For the public version, visit https://www.biostars.org.
GWAS & Logistic regression in R

Hi

I have my genotype data and my phenotype data, I made GWAS with GAPIT tool, but now I would like to use the logistic regression analysis in R, which software is the appropriate ? What would be the first step to work with the data?

Thanks.

r
• 4,472 views • link
updated by Jeremy Leipzig Jeremy Leipzig • written by carlos.maldonadomota • carlos.maldonadomota

1 answer

A logistic regression is where your dependent, outcome variable is dichotomous or binary (aka logical).

GAPIT is a mixed model linear regression, which means it tries to divide the regressors into fixed and random effects. This is super important in agriculture where most of all the GWAS signals end up being SNPs associated with some inbred line.

They correct for this by this assigning fixed effects to SNPs and call the population a random effect, or vice-versa.

GAPIT assumes a continuous linear outcome variable. There is such a thing as mixed model logistic regression. Actually, I wouldn't be surprised if GAPIT has a logistic method hiding somewhere.

So...is your outcome a number or a category?

• 0 views • link
written by Jeremy Leipzig Jeremy Leipzig

It's probably still ok to use GAPIT for this, you just have to round the responses to 0 or 1, assuming these responses are ordered in some meaningful way

• 0 views • link
written by Jeremy Leipzig Jeremy Leipzig

Log in to answer this question.

More posts like this