-
Notifications
You must be signed in to change notification settings - Fork 33
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
pmap and .out blank #113
Comments
Hello,
Yes I believe the reason you don’t get any output is because all your
genetic distances are equal. It is an oversight that selscan is allowed to
function this way without errors or warnings.
When you use —pmap (and provide vcf files) you don’t need to provide any
map file (—pmap is a true/false flag). With —pmap set, selscan will take
the physical distance information from the vcf file.
Hope this helps,
Zachary
Le ven. 17 mai 2024 à 10:16, Tom159357 ***@***.***> a écrit :
… Dear Szpiech
Hello,I encountered some issues when using selscan to calculate XP-EHH。
①I don't quite understand the map file.My map file is generated like
this:vcftools --vcf SNP_merge.vcf --plink --out SNP.The generated file is:1
1:146 0 146.My command is:./selscan --xpehh --vcf --vcf-ref --map --out
.The third column of the map file is all 0, running without any errors, but
the output. xpeh.out is empty,is this because the genetic distance of the.
map file is all 0? Why is the third column set to 0 and this command still
running, but the. xpeh. out file is empty and the. xpeh. log only has
partial bit information (WANING: Reached chrome edge before EHH retired
below 0.05 Skipping calculation at position 146 id: 1:146)?
)②--pmap : Use physical map instead of a genetic map.What is the format of
the pmap file, what is it generated from, and how many columns of
information are there.I want to use pmap instead of map.When I use: selscan
-- xpehh -- vcf-- vcf ref-- pmap-- out, an error will be reported: ERROR:
command line flag chr2pmap.map not recognized
Thank you very much!Looking forward to your reply.
—
Reply to this email directly, view it on GitHub
<#113>, or unsubscribe
<https://github.com/notifications/unsubscribe-auth/ABAKRQXSUTUBXUOZWCK4233ZCYGLDAVCNFSM6AAAAABH4HK4USVHI2DSMVQWIX3LMV43ASLTON2WKOZSGMYDEOBXGM3DKNQ>
.
You are receiving this because you are subscribed to this thread.Message
ID: ***@***.***>
|
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Dear Szpiech
Hello,I encountered some issues when using selscan to calculate XP-EHH。
①I don't quite understand the map file.My map file is generated like this:vcftools --vcf SNP_merge.vcf --plink --out SNP.The generated file is:1 1:146 0 146.My command is:./selscan --xpehh --vcf --vcf-ref --map --out .The third column of the map file is all 0, running without any errors, but the output. xpeh.out is empty,is this because the genetic distance of the. map file is all 0? Why is the third column set to 0 and this command still running, but the. xpeh. out file is empty and the. xpeh. log only has partial bit information (WANING: Reached chrome edge before EHH retired below 0.05 Skipping calculation at position 146 id: 1:146)?
)②--pmap : Use physical map instead of a genetic map.What is the format of the pmap file, what is it generated from, and how many columns of information are there.I want to use pmap instead of map.When I use: selscan -- xpehh -- vcf-- vcf ref-- pmap-- out, an error will be reported: ERROR: command line flag chr2pmap.map not recognized
Thank you very much!Looking forward to your reply.
The text was updated successfully, but these errors were encountered: