Question: merlin .ped file
0
gravatar for User6891
3.1 years ago by
User6891240
Europe
User6891240 wrote:

Hi,

I'm creating a .ped file for Merlin to do linkage analysis for SNParray data.

I have a family of 12 members, but for some siblings, I only have one parent available, so in the .ped file, the father ID might be '8', but the mother id '0' for unknown. 

When I do this, pedstats is giving me an error: Parent named 0 for Person 5 in Family 1 is missing

Any idea on how to solve this?

 

pedigree merlin linkage • 1.6k views
ADD COMMENTlink modified 3.1 years ago by dhibar40 • written 3.1 years ago by User6891240
2
gravatar for dhibar
3.1 years ago by
dhibar40
dhibar40 wrote:

You have to include an entry for each parent in the Ped file even if they are missing. Also, even if a parent is missing you should give them a non-zero ID if the other parent isn't missing. So if the parent (father) is missing from family 1 it currently looks like this:

Family1 Sib1 Father0 Mother8 1

And you need to change it such that the missing father for the group still gets a unique ID:

Family1 Sib1 Father7 Mother8 1

In addition, you need to add the father to your pedigree file:

Family1 Father7 0 0 1

Notice that both parents are missing for Father7 so you can set both to 0 and there isn't the expectation that they will be listed in the Ped file somewhere.

ADD COMMENTlink written 3.1 years ago by dhibar40

But my .ped file, used as input for merlin, also contains the marker genotypes (from the SNP-array). But I have no data for Father7, so can I then just add the first 5 columns and leave the rest empty?

ADD REPLYlink written 3.1 years ago by User6891240

You can provide a "linking file" that lists the full set of subjects like mentioned above with the Ped file that has your genotypes. In merlin you would do this by separating them with a comma:

merlin -p familystructure.ped,pheno.ped -d familystructure.dat,pheno.dat  [...]

ADD REPLYlink written 3.1 years ago by dhibar40

I'm sorry but that's not completely clear for me.

So the pheno.ped should be my original file with the individuals for whom I have marker data?

And then the familystructure.ped should be a .ped file with only 5 columns in which there are lines for all individuals even those without marker data?

pheno.dat is then my original .dat file. But what needs to be inside familystructure.dat?

ADD REPLYlink written 3.1 years ago by User6891240

Problem solved:

pheno.ped = all genotyped individuals (6 columns + marker data)

familystructure.ped = all individuals (6 columns)

familystructure.dat = only one line, with disease condition (A condition)

ADD REPLYlink written 3.1 years ago by User6891240
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: 1400 users visited in the last hour