Skip to content

error in normalize_binless #20

Description

@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?

Activity

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

      Milestone

      No milestone

      Relationships

      None yet

      Development

      No branches or pull requests

      Issue actions