Skip to content
Snippets Groups Projects
Commit db955384 authored by Sascha Meiers's avatar Sascha Meiers
Browse files

IMPORTANT UPDATE: don't really *force* bi-allelic, just give a prior aginst other states

parent 480a3cdf
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
utils/mosaiClassifier/makeSVcalls.R
+ 1
− 1
View file @ db955384
...@@ -56,7 +56,7 @@ makeSVCallSimple <- function(probs, llr_thr = 1) { ...@@ -56,7 +56,7 @@ makeSVCallSimple <- function(probs, llr_thr = 1) {
return(probs[sv_call_name != "ref_hom" & llr_to_ref > llr_thr]) return(probs[sv_call_name != "ref_hom" & llr_to_ref > llr_thr])
} }
forceBiallelic <- function(probs, penalize_factor=0) forceBiallelic <- function(probs, penalize_factor = 0.1)
{ {
# Add reference probability as an extra column (same as in makeSimpleSVCalls) # Add reference probability as an extra column (same as in makeSimpleSVCalls)
probs[, probs[,
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment