hakyimlab/PrediXcan

How to handle the missing dosages?

Closed this issue · 2 comments

Hi Predixcan guys,
How to handle the missing dosages, while creating the predixcan format dosage files ? If its 'NA', the script throws error saying cannot convert string to float NA.
Thanks
Veera

Hi Veera

The best option is to impute the genotype and not have the NA. A very
simple though not optimal is to impute with the mean of the population.
HTH
Haky

On Thursday, November 3, 2016, Veera notifications@github.com wrote:

Hi Predixcan guys,
How to handle the missing dosages, while creating the predixcan format
dosage files ? If its 'NA', the script throws error saying cannot convert
string to float NA.
Thanks
Veera


You are receiving this because you are subscribed to this thread.
Reply to this email directly, view it on GitHub
#25, or mute the thread
https://github.com/notifications/unsubscribe-auth/AC2oub8JQbww-ZJejcgqDlE4QJ44P9Svks5q6dv-gaJpZM4KoVOW
.

Hi Haky,
Thanks for the quick response. I guess that's what the convert_plink_to_dosage.py does, right ? (I see that the script replaces NA with MAF*2.

-Veera