Admixture cv error
2
0
Entering edit mode
8.6 years ago
s.bourne • 0

Hi all,

When doing Admixture and outputting cv error rate for a range of k values, one of the k values keeps on outputting "nan" as the error rate. It works for the other k values, but not for this one. Has this happened to anyone else?

Cheers,
Steve

Structure CV-Error-Rates Admixture • 5.3k views
ADD COMMENT
0
Entering edit mode

Hi Steve, did you ever figure this out?? I get the same issue every time I run K=3 on my data.

Open to suggestions!

ADD REPLY
0
Entering edit mode

Same here. Here's some sample output from the bottom of the file. I'm not so worried about the numbers as I am about the fact that the job isn't finishing. It's been running for four hours and all the other versions are done, including ones with higher K values.

586 (QN/Block)  Elapsed: 70.414 Loglikelihood: -nan     (delta): -nan
587 (QN/Block)  Elapsed: 70.526 Loglikelihood: -nan     (delta): -nan
588 (QN/Block)  Elapsed: 70.537 Loglikelihood: -nan     (delta): -nan
589 (QN/Block)  Elapsed: 70.47  Loglikelihood: -nan     (delta): -nan
ADD REPLY
0
Entering edit mode

I'm running into the same issue. I am receiving 'nan' for Loglikelihood and delta for K=6 only. K=1-5 works fine as does K=7. This is quite odd. Any thoughts?

Chris

ADD REPLY
0
Entering edit mode
7.8 years ago
chengcw • 0

You may want to run with admixture32 instead. Hope it helps.

ADD COMMENT
0
Entering edit mode
18 months ago
Declan • 0

I had this issue today and was able to resolve it by changing the random seed Admixture uses by adding -s time (which sets the random seed based on the current time) to my command so it looks something like this:

admixture --cv -s time inputFile.ped 6

I'm not sure why this works, but so far it has fixed the issue for me. I previously had tried restarting the K=6 run several times and the seed was always the same, so it always gave me the nan Loglikelihood and delta values. I figured trying a different starting point might fix things.

-Declan

Edit: While continuing to run Admixture with other data, the issue still sometimes persists with this solution in place. However, stopping the run and restarting to get a new seed seems to help. It's a bit hit or miss.

ADD COMMENT

Login before adding your answer.

Traffic: 1354 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6