error in normalize_binless

[mtihme]

Hi, I am trying to run optimized_binless protocol. I successfully run preprocessing.R, At the third iteration of normalize_binless I have a reproducible error:

Iteration 1

Residuals Exposures log-likelihood = 22.44619 Dispersion fit: alpha 108913.2 log-likelihood = -19171238 Genomic group 1 : 10 steps lambda_iota 0.01666102 lambda_rho 0.0151357 group 2 : 7 steps lambda_iota 0.01875217 lambda_rho 0.02100009 log-likelihood = -239971.1 Decay group 1 : 3 steps lambda_diag 0.1814284 log-likelihood = -111541.3

Iteration 2

Residuals Exposures log-likelihood = 33.22606 Dispersion fit: alpha 0.01 log-likelihood = -4687710 Genomic group 1 : 42 steps lambda_iota 1.693895e-05 lambda_rho 1.143178e-05 group 2 : 25 steps lambda_iota 0.0003145 lambda_rho 0.0001278387 log-likelihood = -348224061259 Decay group 1 : 1 steps lambda_diag 0.1268076 log-likelihood = -244909

Iteration 3

Residuals Exposures log-likelihood = 15.92615 Dispersion Error in { : task 7 failed - "missing value where TRUE/FALSE needed" Calls: normalize_binless ... system.time -> -> %dopar% -> Timing stopped at: 36.34 18.49 12.62 Execution halted

Could you please help me to solve the error?

[mtihme]

I change base.res from 5000 to 10000 And get the error at 5-th iteration

Iteration 5

Residuals Exposures log-likelihood = 15.87611 Dispersion Error in { : task 1 failed - "missing value where TRUE/FALSE needed" Calls: normalize_binless ... system.time -> -> %dopar% -> Timing stopped at: 36.01 26.1 13.29 Execution halted