Skip to content

error in normalize_binless #20

@mtihme

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?

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