Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- > head(df)
- POSITION MEANDEPTH CHROM
- 1 0:10000 0 chr1
- 2 10000:20000 0 chr1
- 3 20000:30000 0 chr1
- 4 30000:40000 0 chr1
- 5 40000:50000 0 chr1
- 6 50000:60000 0 chr1
- > tail(df)
- POSITION MEANDEPTH CHROM
- 308834 57170000:57180000 0 chrY
- 308835 57180000:57190000 0 chrY
- 308836 57190000:57200000 0 chrY
- 308837 57200000:57210000 0 chrY
- 308838 57210000:57220000 0 chrY
- 308839 57220000:57230000 0 chrY
- > levels(df$CHROM)
- [1] "chr1" "chr10" "chr11" "chr12" "chr13" "chr14" "chr15" "chr16" "chr17" "chr18" "chr19" "chr2" "chr20" "chr21" "chr22" "chr3" "chr4"
- [18] "chr5" "chr6" "chr7" "chr8" "chr9" "chrM" "chrX" "chrY"
- # RE_ORDER CHROMS
- chrom_order<-c('chr1','chr2','chr3','chr4','chr5','chr6','chr7','chr8','chr9','chr10','chr11',
- 'chr12','chr13','chr14','chr15','chr16','chr17','chr18','chr19','chr20','chr21','chr22','chrX','chrM')
- df<- df[match(chrom_order, df$CHROM),]
- > head(df)
- POSITION MEANDEPTH CHROM
- 1 0:10000 0 chr1
- 128716 0:10000 0 chr2
- 169134 0:10000 0 chr3
- 188964 0:10000 0 chr4
- 207986 0:10000 0 chr5
- 226140 0:10000 0 chr6
- library(dplyr)
- df %>%
- slice(match(CHROM, chrom_order))
Add Comment
Please, Sign In to add comment