This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Merge Two Plink Transposed Files

I want to look at concordance by SNP using Plink, however my two datasets are both in Plink TPED/TFAM format. the documentation on Plink merge covers regular ped and map files. Is there a way to merge these datsets as is or do I need to convert to ped/map before merging.

plink gwas
• 10,002 views • link
updated by Pablo Marin-Garcia Pablo Marin-Garcia • written by Zach Stednick Zach Stednick

How about using "--recode" to make it SNP major format and then do a regular plink merge?

• 0 views • link
written by jvijai jvijai

Thank you, this works ... I am not sure why I did not think of this earlier.

• 0 views • link
written by Zach Stednick Zach Stednick

1 answer

You can not use the tped in the merge. In order to save time and disk space I would recommend you to convert to bed and then do the merge-comparison.

Having

  • myfile1.tped
  • myfile2.tped

I would do the comparison as follow:

plink --noweb --tfile myfile1 --make-bed --out myfile1
plink --noweb --tfile myfile2 --make-bed --out myfile2

# then the comparison
plink --noweb --bfile myfile1 --bmerge myfile2.bed myfile2.bim byfile2.fam --merge-mode 6 --out myfile1vs2
• 1 views • link
written by Pablo Marin-Garcia Pablo Marin-Garcia

Hi, I tried this command mentioned above to calculate the concordance rate plink --noweb --bfile file1 --bmerge file2.bed file2.bim file2.fam --merge-mode 6 --out myfile1vs2 but it shows "Concordance rate is nan" -> what does this mean ? Am i going wrong somewhere ?

• 0 views • link
written by P.NJ • P.NJ

nan means "not a number" Are you sure that you have any SNP-ind genotype in common between datasets?. It could be that it can not calculate the 'rate' of concordance out of total in common because no genotypes in common so the rate is a division by 0. This is only speculation, I have never seen this error before.

• 0 views • link
written by Pablo Marin-Garcia Pablo Marin-Garcia

Log in to answer this question.

More posts like this