Question: plink mergelist put many individuals into same binary file set
0
gravatar for gulcek
5.5 years ago by
gulcek20
Turkey
gulcek20 wrote:

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 • 1.8k views
ADD COMMENTlink modified 9 months ago by Dawn0 • written 5.5 years ago by gulcek20

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

ADD REPLYlink written 9 months ago by Dawn0

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.

ADD REPLYlink modified 9 months ago • written 9 months ago by ATpoint40k

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.

ADD REPLYlink written 9 months ago by Dawn0

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
ADD REPLYlink written 9 months ago by Kevin Blighe66k

New question: How to merge 23andme samples in Plink? Each .ped/.map file represent a sample.

ADD REPLYlink written 9 months ago by Kevin Blighe66k
1
gravatar for chrchang523
5.5 years ago by
chrchang5237.3k
United States
chrchang5237.3k wrote:

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.

ADD COMMENTlink written 5.5 years ago by chrchang5237.3k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1183 users visited in the last hour