Skip to content
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 in normalize_binless #20

Open
mtihme opened this issue Jun 14, 2020 · 1 comment
Open

error in normalize_binless #20

mtihme opened this issue Jun 14, 2020 · 1 comment

Comments

@mtihme
Copy link

mtihme commented Jun 14, 2020

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
Copy link
Author

mtihme commented Jun 14, 2020

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

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant