Question: Analysis Of Imputated Genetic Data Using Plink
3
gravatar for Jim
9.1 years ago by
Jim30
Jim30 wrote:

PLINK now takes imputation data for analysis by using --dosage. I suppose it automatically runs additive model. I would like to analyze data with recessive model and have tried to put --recessive in the commands. But the result turns out to be the same. So I would like to know how to correctly run recessive model for imputation data. thank you.

ADD COMMENTlink modified 9.1 years ago by Genotepes950 • written 9.1 years ago by Jim30
2
gravatar for Genotepes
9.1 years ago by
Genotepes950
Nantes (France)
Genotepes950 wrote:

As Larry says,

the dosage data makes it nearly impossible to test additive of dominant model. I guess you could put a threshold on the allele dosage and classify as homozyguous or not, but then you need to go out of plink - possible in R.

If you have genotype probabilities, not dosage - as most of the imputation programs can give if you ask them, then SNPTEST could run a recessive model by running a logistic regression where the likelihood is weighted by the genotype probability. In my experience, you need quite a lot of individuals to have enough rare homozyguous and get results - and the SNP needs to be frequent.

So anyways, to help you, we need some info about your input data.

Christian

ADD COMMENTlink written 9.1 years ago by Genotepes950
1
gravatar for Larry_Parnell
9.1 years ago by
Larry_Parnell16k
Boston, MA USA
Larry_Parnell16k wrote:

The very nature of using imputed data, which are a form of dosage effects correlates with an additive model, dictates that testing with a recessive model is not possible. In order to run the recessive model, it is necessary to use only the original genotypes and no imputed data.

Once you have preliminary associations from the original genotypes, you can consider a type of fine mapping to further define the SNPs in a recessive model. This is not simple as the imputed must be treated by setting a threshold (which I myself do not understand) before you can combine with the original genotypes. Nonetheless, this is possible.

ADD COMMENTlink written 9.1 years ago by Larry_Parnell16k
Please log in to add an answer.

Help
Access

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