-
Notifications
You must be signed in to change notification settings - Fork 0
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
unique markers #2
Comments
Hi,
The problem comes from having the genetic distance as all zero. I suggest
using the average recombination rate to estimate the genetic distances and
using that in the genetic distance column in the input file. Hope this
helps!
Best,
Amy
…On Wed, Apr 17, 2019 at 1:42 AM Dr Lauren C White ***@***.***> wrote:
Hi,
I've been trying to run CLAPPER, but I get an error that I cannot figure
out how to solve.
"Thre are more than one markers in the same position in filename.tped!
Every marker must have a unique position."
I have double checked that all positions are unique. However, I don't have
genetic distances in my tped file. This values of this column are all set
to 0. (in my options file #conditiononld is also set to 0). Could this be
the problem? Can you suggest a way around it?
Cheers
—
You are receiving this because you are subscribed to this thread.
Reply to this email directly, view it on GitHub
<#2>, or mute the thread
<https://github.com/notifications/unsubscribe-auth/AFboJ78w2r_D5f9_JfH5Re5d0Td_UlAHks5vht5cgaJpZM4c0ghw>
.
|
Thanks for the quick reply! And yes, that's really helpful. Cheers |
Hi there, I'm having the same issue. I am not familiar with cM distances and recombination. My PLINK files don't have recombination rates, so I have to add them. I fetched them from : http://bochet.gcc.biostat.washington.edu/beagle/genetic_maps/plink.GRCh38.map.zip and I converted them to a format PLINK likes, using the recombination rate average of 1.2: for i in {1..22}; do awk '{print $4, "1.2", $3*100}' plink.chr$i.GRCh38.map | sponge plink.chr$i.GRCh38.map; done This gives a file like this per chromosome (position, recomb_rate, Morgan distance):
I then annotate my file with I still get the same error. What am I doing wrong? Thanks, A |
So, it turns out that if you have sequencing data, you might have consecutive variants that are really close, and therefore they might have the same coordinates in Morgan. I went around this by removing variants that have the same coordinates. Another way would be to add a random small noise to the duplicates using R. I also ran into another error after that saying that my coordinates were decreasing. It turns out the chromosome map resets after every chromosome end, whereas I wonder if it would be possible to use chromosome and position instead of cM distance in a future release. Files annotated with cM distance are now becoming rarer and rarer, and recombination maps are hard to find for new builds. |
Hi,
I've been trying to run CLAPPER, but I get an error that I cannot figure out how to solve.
"Thre are more than one markers in the same position in filename.tped! Every marker must have a unique position."
I have double checked that all positions are unique. However, I don't have genetic distances in my tped file. This values of this column are all set to 0. (in my options file #conditiononld is also set to 0). Could this be the problem? Can you suggest a way around it?
Cheers
The text was updated successfully, but these errors were encountered: