Skip to content
Jump to bottom

correct maxpos calculation #53

Merged
merged 2 commits into from
Jan 10, 2019
Merged

correct maxpos calculation #53

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
b128196
correct maxpos calculation
JannikHamp Jan 10, 2019
5343f60
fixed for correct contains_maxpos calculation
JannikHamp Jan 10, 2019
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Jump to
Jump to file
Failed to load files.
Diff view
Diff view
6 changes: 3 additions & 3 deletions bin/1.2_filter_motifs/compareBed_runinfo.R
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -8,7 +8,7 @@
# adds a column with a flag "contains_maxpos"
# and creates a file with information of the bedtool comparison

# parameters: Parameters are not named. Respect the parameter order.
# parameters: Parameters are not named, respect the parameter order.
# min: minimum footprint size threshold
# max: maximum footprint size threshold
# input_raw: unfiltered BED-file
Expand Down Expand Up @@ -52,10 +52,10 @@ data_filtered[[7]] <- data_filtered[[3]] - data_filtered[[2]]

# adding column "contains_maxpos", containing flag (0 or 1)
# max_pos is the position of maximum score of a footprint
data_filtered[[8]] <- data_filtered[[8]] - data_filtered[[2]]
data_filtered[[8]] <- data_filtered[[8]] + data_filtered[[2]]
data_filtered <- cbind(data_filtered, contains_maxpos = 0)
data_filtered$contains_maxpos[intersect(which(data_filtered[[2]] <= data_filtered[[8]]), which(data_filtered[[3]] > data_filtered[[8]]))] = 1
data_filtered[[8]] <- data_filtered[[8]] + data_filtered[[2]]
data_filtered[[8]] <- data_filtered[[8]] - data_filtered[[2]]

fwrite(data_filtered, output, col.names=FALSE, quote = FALSE, sep = '\t')

Expand Down