-
Notifications
You must be signed in to change notification settings - Fork 1
Description
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?