This is a test version of Biostars. For the public version, visit https://www.biostars.org.
plink mergelist put many individuals into same binary file set

We have some difficulties about using 23andMe data on plink. We have raw 23andMe snp data for each individual (total of 756). We converted each individual file to binary fileset .bed,.bim, .fam via plink.

plink --23file $file --out $a

After creating binary set for each individual, we converted them into .map and .bed format.

plink --bfile binary_filesets --a2-allele ../All.vcf 4 3 '#' --recode --out $a

Then, we want to collect info of all individuals into same .bed file( or binary fileset.)

plink --file genom0 --merge-list ../filesets3.txt --make-bed --out ../combined/combdata

We assume --merge-list does this job, however, it does overwrite one individual into another. It doesn't create many individuals x info (matrix-like) data. How could we obtain one combined binary fileset for many individuals?

plink merge bed
• 3,563 views • link
updated by Ram Ram • written by gulcek • gulcek

Hello! I am also new to Plink. I encountered the same problem. How did you solve it?

• 0 views • link
written by Dawn • Dawn

Did you follow the answer of the developer chrchang523 below? Also please leave the answer field for actual answers and use the comment buttom instead.

• 0 views • link
written by ATpoint ATpoint

Yes, I want to follow the answer of the developer chrchang523. But since I am very new to Plink, I don't know how to assign a different ID to each sample. Can you give me some instructions? Thank you.

• 0 views • link
written by Dawn • Dawn

I don't know how to assign a different ID to each sample.

This intention was not stated in your original question. Please try to explain again exactly what you want to do. State:

  • what data you have
  • what you want to do
  • what you have already tried and the error messages
• 0 views • link
written by Kevin Blighe Kevin Blighe

1 answer

  1. --a2-allele is pointless with --recode, since the .ped/.map format does not track reference alleles. You must use --make-bed instead of --recode.
  2. You need to assign a different ID to each sample (this is easy to do during the --23file step). Otherwise --merge-list assumes all the genotype calls are supposed to be for the same sample.
• 0 views • link
updated by Ram Ram • written by chrchang523 chrchang523

Log in to answer this question.

More posts like this