syntheke / bayesR

Bayesian hierarchical model for complex trait analysis
45 stars 26 forks source link

Missing phenotypic values in .fam file #3

Closed zhiqiuhu closed 6 years ago

zhiqiuhu commented 8 years ago

Based on the 'BayesRmanual-0.75.pdf' file, which saying the missing phenotypes to be coded as ‘NA’, I believe missing values are allowed in the .fam file.

However, as I included missing values in the fam (denoted as NA), BayesRv2 generated the following error message

forrtl: severe (24): end-of-file during read, unit -5, file Internal List-Directed Read Image PC Routine Line Source
BayesRv2 000000010EA2E15A Unknown Unknown Unknown BayesRv2 000000010EA56B97 Unknown Unknown Unknown BayesRv2 000000010EA55582 Unknown Unknown Unknown BayesRv2 000000010E9FE812 Unknown Unknown Unknown BayesRv2 000000010EA0D16D Unknown Unknown Unknown

BayesRv2 000000010E9F74AE Unknown Unknown Unknown

I tried to generate a new file by keeping only the individuals with valid phenotypic values and BayesRv2 works fine. To bypass the issue, I simply generated new plink binary files for each trait. But it would be great if anyone would check whether it is a common issue or BayesR.

FYI, my system is OSX.

Thanks