I've been simulating some RNAseq data using your package and I started getting errors like the following.
Error in .Call2("C_solve_user_SEW", refwidths, start, end, width, translate.negative.coord, :
solving row 206369: negative values are not allowed in 'width'
Calls: simulate_experiment ... make_IRanges_from_windows_args -> solveUserSEW -> .Call2
Execution halted
I am using distr="custom" and the only thing changing between successful and failed simulations is the shape of the logspline distribution. Stopping rlogspline from producing negative fragment sizes appears to have solved the problem for me.
Let me know if you're interested in the pull request and whether you'd like a reporducible example.
Hi @alyssafrazee ,
I've been simulating some RNAseq data using your package and I started getting errors like the following.
Error in .Call2("C_solve_user_SEW", refwidths, start, end, width, translate.negative.coord, : solving row 206369: negative values are not allowed in 'width' Calls: simulate_experiment ... make_IRanges_from_windows_args -> solveUserSEW -> .Call2 Execution halted
I am using distr="custom" and the only thing changing between successful and failed simulations is the shape of the logspline distribution. Stopping rlogspline from producing negative fragment sizes appears to have solved the problem for me.
Let me know if you're interested in the pull request and whether you'd like a reporducible example.
Kind regards,
Scott