Detailed description |
GenABEL v. 1.7-4 (February 22, 2013)
When exporting a subset of individuals from gwaa.data with export.merlin, the genotypes in ped files are incorrect.
Additionally, ped files produced by analogous commands
export.merlin(df[ids, snps], datafile=data_file_name, pedfile=ped_file_name, mapfile=map_file_name, order=T, fixstrand="+", stepids=1)
and
export.merlin(df[ids, snps], datafile=data_file_name, pedfile=ped_file_name, mapfile=map_file_name, order=T, fixstrand="+", stepids=100)
are not identical.
When stepids=1, the ped file seems to be correct.
This behavior is observed only with subset of ids. When gwaa.data is exported completely, the everything is ok. |
|