Skip to content

Commit

Permalink
Update BISCUT_peak_finding.R
Browse files Browse the repository at this point in the history
Removing problematic and unnecessary print statements
  • Loading branch information
shahab-sarmashghi authored Aug 10, 2023
1 parent b35ffac commit 71bb37a
Showing 1 changed file with 4 additions and 4 deletions.
8 changes: 4 additions & 4 deletions BISCUT_peak_finding.R
Original file line number Diff line number Diff line change
Expand Up @@ -490,15 +490,15 @@ do_arm_gistic <- function(arm,tt,direc,telcent,mode,prefix=NULL,leftright=NULL,d
#print(probelist)
sorted_probelist<-probelist[order(probelist)]
left_closest_index <- which(abs(sorted_probelist-left_boundary_coord)==min(abs(sorted_probelist-left_boundary_coord)))
print(left_closest_index)
#print(left_closest_index)
if(left_boundary_coord!= sorted_probelist[left_closest_index]) {
left_boundary_coord <- sorted_probelist[left_closest_index-1][1]
} else {
left_boundary_coord <- sorted_probelist[left_closest_index][1]
}

right_closest_index <- which(abs(sorted_probelist-right_boundary_coord)==min(abs(sorted_probelist-right_boundary_coord)))
print(right_closest_index)
#print(right_closest_index)
if(right_boundary_coord!= sorted_probelist[right_closest_index]) {
right_boundary_coord <- sorted_probelist[right_closest_index+1][1]
} else {
Expand Down Expand Up @@ -590,8 +590,8 @@ do_arm_gistic <- function(arm,tt,direc,telcent,mode,prefix=NULL,leftright=NULL,d
if(mode=='overlap') {
rightdf <- df1[vector_of_interest>left_boundary,]
leftdf <- df1[vector_of_interest<right_boundary,]
print(prior_peaks)
print(left_boundary, right_boundary)
#print(prior_peaks)
#print(left_boundary, right_boundary)
if(check_overlaps(prior_peaks,c(left_boundary,right_boundary))) {
stop('repeat found!')
}
Expand Down

0 comments on commit 71bb37a

Please sign in to comment.