-
Notifications
You must be signed in to change notification settings - Fork 6
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
ERROR: [__FUNCTION__] invalid allele frequencies in 1st iteration of ngsF-HMM #7
Comments
Could you try running with |
Hi, fgvieira, Thank you for the quick response. The option of "--freq_est 1" did run. It took over 3 days to finish. I am checking the inbreeding coefficient data now. Meanwhile I have some questions on these runs as welll as using ngsF-HMMplot.R.
Here is the command I used:
where:
I don't have a .path.gz file in my output. What is that file? How can I produce one? Thank you in advance for your help! |
Hi @wul88 and sorry for the late reply.
|
Hi, fgvieira, As running ngsF-HMM, I have following concerns. Please advise me.
|
|
Hello, fgvieira
I tried to run ngsF-HMM to find IBD of samples from different sub-populations. I used the Beagle GL files computed previously as below:
I have .beagle.gz and .mafs.gz files. I extracted first two columns from .mafs.gz file and removed the header, then saved as pos file. The nsites was computed by "cat $pos | wc -l " . The command to run ngsF-HMM is as follows:
WHERE
glf2 is the Beagle format glf from angsd run shown at the beginning, and inF_transitPar_118AveAuto.txt is a two-column file: 1st column is the inbreeding coefficient from running ngsF and the second column is 0.1 as transition parameter. one row per individual. ( I also used 0.1, 0.1 as shown in tutorial in different run, received the same error.)
AND
thisSeed=$RANDOM
glf2="all118_MafGlf2_allChrs.beagle.gz"
#maf="all118_MafGlf2_allChrs.mafs.gz"
pos="all118_MafGlf2_allChrs_noHeader.pos"
NSAMS=118
NSITES=$((
cat ${pos} | wc -l
))outDir=ngsFHMM/"rep_"$thisRep
out="all118depthSet_inFHMM_"$thisRep
( thisRep is replicate number. I ran 10 replicates but every one had the same error.)
--- one of the .err files ---->
thisRep=2
thisSeed=18950
glf2=all118_MafGlf2_allChrs.beagle.gz
pos=all118_MafGlf2_allChrs_noHeader.pos
NSAMS=118
outDir=ngsFHMM/rep_2
NSITES=2205512
out=all118depthSet_inFHMM_2
ngsF-HMM --geno all118_MafGlf2_allChrs.beagle.gz --lkl --pos all118_MafGlf2_allChrs_noHeader.pos --n_ind 118 --n_sites 2205512 --freq 0.1 --freq_est 2 --indF inF_transitPar_118AveAuto.txt --out ngsFHMM/rep_2/all118depthSet_inFHMM_2 --n_threads 8 --seed 18950 --min_epsilon 1e-9 --verbose 2 --log 1
Could you please give me some advice as what I did wrong? Thank you very much in advance!
The text was updated successfully, but these errors were encountered: