Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- require(dplyr)
- require('magrittr')
- in_file = 'WT_JC4.pausing.2018_05_31.csv'
- df = read.delim(in_file, header=TRUE, sep='\t')
- maxes <- group_by(df, gene) %>% summarize(max=pi=max(pausing_index))
- write.table(maxes, 'out.tmp', row.names=FALSE, sep='\t', quote=F)
Add Comment
Please, Sign In to add comment