Skip to content
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

Assigment of specific haplogroups #2

Open
cdquinto opened this issue Sep 24, 2018 · 1 comment
Open

Assigment of specific haplogroups #2

cdquinto opened this issue Sep 24, 2018 · 1 comment

Comments

@cdquinto
Copy link

Hello, I'm using this program to assign haplogroups to some 1000 Genomes samples, but I am only getting the major haplogroups.

According to the paper: "YHap can be customized to either only assign to within the major haplogroups (A through to T), or all possible haplogroups at the finest level of classification." However I don't see where I can change this in the program.

Any help is greatly appreciated!

@Griffan
Copy link
Owner

Griffan commented Sep 24, 2018

Hi, I should have written readme more in detail. Currently, the best guideline is to follow the "autorun.advanced.pl" in rundir, which is a wrapper to search in branches you specified to.
Regarding your question, you can find in the script:

sub runJava()
{
    my ($model)=@_;
    my $cmdline="nohup sh ../../runYhap.sh ".$pop." ".$model." 1>&2 "." 2>>".$pop.".e.log";
    print STDERR "executing $cmdline ...\n";
    return system($cmdline);
}

the variable $model is where you should specify the specific branch you want to search, the more detailed branches you specified, the more the detailed assignment will be.
at the same time, for each run, the less complicated the model is, the more accurate the result will be.
you can find that in the "autorun.advanced.pl" script, the search started from ABCDEF these main branches, and in each next round, the model will be limited within the sub-branches of branches of previous round.

Hope this resolve your problem, feel free to contact me. I'm glad to help you successfully get your result.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants