shapeit5
shapeit5 copied to clipboard
phase_common failing at PBWT step
I'm using SHAPEIT5 to phase a large panel (300K+) using phase_common and it keeps stopping at the PBWT step? This is the relevant error:
WARNING: No PBWT states found [1000042 / w=76] / Using 100 random states
WARNING: No PBWT states found [1000042 / w=77] / Using 100 random states
ERROR: Haploid underflow impossible to recover for [1000039]
The parameters used are below:
Parameters:
- Seed : 15052011
- Threads : 8 threads
- MCMC : 15 iterations [5b + 1p + 1b + 1p + 1b + 1p + 5m]
- PBWT : [window = 4cM / depth = auto / modulo = auto / mac = 5 / missing = 0.1]
- HMM : [window = 4cM / Ne = 15000 / Recombination rates given by genetic map]
Please advise?