Question: Best implementation of negative binomial association test for GWAS data?
gravatar for Vincent Laufer
4.3 years ago by
Vincent Laufer1.0k
United States
Vincent Laufer1.0k wrote:

I am being asked to implement a negative binomial association test on Immunochip data.

I ran the logistic model on a dichotimization of the response variable, but we did not find any associations at genome-wide levels of significance, so I am now being ask to consider the response variable as a continuous variable and test again using the negative binomial.

Does anyone know which are the best implementations of such a test for genetic data? Currently the files are in Plink format, although I can convert them if necessary.

Thank you.

ADD COMMENTlink modified 19 months ago by Biostar ♦♦ 20 • written 4.3 years ago by Vincent Laufer1.0k
gravatar for brentp
4.3 years ago by
Salt Lake City, UT
brentp23k wrote:

Negative binomial is for modeling count data. What are your counts?

If you dichotomized the response variable as you say, then logistic regression makes sense. If you want to re-try on the continuous response variable, then just run a linear model.

Is your model like:

     response ~ genotype + age + ...


ADD COMMENTlink written 4.3 years ago by brentp23k

yes, it is:

radiographic severity ~ genotype + age + sex + PC1 + ... 

so I gather from your answer that you would just go with a standard R package, in other words, there is not a significantly better test thats adapted to genetic data?

thank you.

ADD REPLYlink written 4.3 years ago by Vincent Laufer1.0k
gravatar for jstats
3.9 years ago by
New Zealand
jstats0 wrote:

Not sure if you still want this, but I have R code to perform a neg binomial GWAS. Our data had lots of zero counts, but we felt that going to the full zero inflated models was a step too far.


ADD COMMENTlink written 3.9 years ago by jstats0

Hi Jstats,


Just seeing this. May I ask you why you thought it was going too far?

Also, I would take the R code for the association testing although I already wrote my own ..  but it would be good for comparison.

However, if you can, what I would be more interested in is your procedure in model selection and decision making process in using negative binomial regression instead of ZINB. I dont have access to any code written in R that walks me through that part of it ...

Thanks very much.

ADD REPLYlink modified 3.8 years ago • written 3.8 years ago by Vincent Laufer1.0k
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.3.0
Traffic: 1429 users visited in the last hour