privefl2 comments

Repositories
Issues
Comments

Results 4 comments of


                                            privefl2

Unexcepted output with modkit adjust-mods --convert h m

Thank you for the detailed explanation. I was wrongly assuming the intervals to be centered aroung q/255 instead of (q+0.5) / 256. So, when we sum p1 + p2 =...

Unexcepted output with modkit adjust-mods --convert h m

But, how would you be able to update this? Either you use q1+q2+1 and you are upward-biased, or use q1+q2 and are downward-biased, right?

Unexcepted output with modkit adjust-mods --convert h m

If you use q1+q2, then the probability is downward-biased by 0.5/256. Maybe the best solution is to use q1+q2 when 127, but maybe it is too much of a hassle...

Question: position in extract full

Sorry, the full code is this: ```r library(furrr) NCORES % # all are "+" so haven't tested with mod_strand == "-" yet mutate(pos = ifelse(mod_strand == ref_strand, pos, ifelse( ref_strand...