Entering edit mode
12.1 years ago
Tyro
•
0
I am trying to input a file using the compound-genotypes option (as well as several other options), but I get the following error message from PLINK:
ERROR: Problem with compound genotype [ ] should be two characters long
My command line is:
plink --1 --no-fid --no-parents --no-sex --compound-genotypes --map3 --file test
The PED file is:
1 1 gg gg gg gg AA AA AA AA AA gg
2 0 AA gg gg AA gg gg gg AA AA AA
3 0 AA AA AA gg AA gg AA gg gg AA
4 1 gg AA gg AA gg gg gg gg AA AA
5 1 gg AA gg AA AA AA gg gg AA AA
6 1 gg gg AA gg AA gg gg AA gg AA
7 1 gg gg gg gg AA gg gg AA AA AA
8 0 AA gg AA gg AA gg gg AA AA gg
9 0 AA AA gg gg gg AA gg gg gg gg
10 1 gg AA gg AA AA AA gg AA AA gg
The MAP file is:
0 X[0] 0
0 X[1] 0
0 X[2] 0
0 X[3] 0
0 X[4] 0
0 X[5] 0
0 X[6] 0
0 X[7] 0
0 X[8] 0
0 X[9] 0
I am running v1.07 on a Windows XP platform.
Does anyone know what I'm doing wrong? This is just a toy example; I need to convert many large data files to PED/MAP and I'm totally stuck now.
Thanks!
I've narrowed it down to just the --compound-genotypes option. If I remove that option, and add spaces between the alleles, it parses fine. Orthogonally, if I add the missing columns and get rid of the corresponding options, the --compound-genotypes still does not work.
Oops. In my last comment I should have made it clear that I still need help understanding why the --compound-genotypes option doesn't work.
Is your map file formatted correctly? I have never seen a plink map without positional data.