We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
1 parent b499f5c commit da88754Copy full SHA for da88754
ascets_resources.R
@@ -85,6 +85,12 @@ make_arm_call <- function(df, thresh) {
85
max_alt_type <- df[df$alt_frac == max_alt,]$alt
86
87
# assign call based on alteration threshold
88
+ if(length(max_alt_type) > 1 & ("DEL" %in% max_alt_type & "AMP" %in% max_alt_type)) { # edge case where two types affect exactly the same fraction of an arm
89
+ return(data.frame(sample = as.character(df$sample[1]), ARM = as.character(df$arm[1]), CALL = "NC"))
90
+ } else if (length(max_alt_type) > 1) {
91
+ max_alt_type <- max_alt_type[max_alt_type != "NEUTRAL"]
92
+ }
93
+
94
if(max_alt_type == "NEUTRAL") {
95
if(max_alt >= thresh) call <- max_alt_type
96
else call <- "NC"
0 commit comments