What Does Nan Mean In Plink'S Sex Check?
1
1
Entering edit mode
12.3 years ago
nnlnn ▴ 60
        FID         IID       PEDSEX       SNPSEX       STATUS            F
  aaa100000   aaa100000            2            0      PROBLEM          nan
  bbb100000   bbb100000            2            0      PROBLEM       0.4208

I got the above results by doing sex check as

plink --noweb --bfile mydat --check-sex --out mydat

I wonder what nan means in the output. Seems like it is because inbreeding coefficient F cannot be computed but why is that? Thanks for your help.

plink • 6.4k views
ADD COMMENT
2
Entering edit mode

nan or NaN or NA in computing generally refer to something that is "Not A Number". When running statistical tests, NaN can be returned when the statistics cannot be computed, perhaps due to missing values in the data or inappropriately coded data.

ADD REPLY
2
Entering edit mode

Check the X chromosome genotype calls for that individual. Missing data often hoses things.

ADD REPLY
0
Entering edit mode

Thanks for the hint. Yes, those with nan F values had completely missing genotypes on the X chromosome.

ADD REPLY
4
Entering edit mode
12.3 years ago
pagel ▴ 90

I don't have any experience with this package in particular, but I believe you are correct: it cannot be computed. nan is usually written as NaN and means "not a number." Most likely this results from a division by zero error. Is the inbreeding coefficient of that first individual not calculable because one of its parents is not in the database?

ADD COMMENT
1
Entering edit mode

Or anything that created an underflow/overflow issue where a number got rounded to 0 due to it being a very small probability for instance.

ADD REPLY

Login before adding your answer.

Traffic: 2681 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6