Guest User

Untitled

a guest
Apr 24th, 2018
77
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.30 KB | None | 0 0
  1. > head(df)
  2. POSITION MEANDEPTH CHROM
  3. 1 0:10000 0 chr1
  4. 2 10000:20000 0 chr1
  5. 3 20000:30000 0 chr1
  6. 4 30000:40000 0 chr1
  7. 5 40000:50000 0 chr1
  8. 6 50000:60000 0 chr1
  9. > tail(df)
  10. POSITION MEANDEPTH CHROM
  11. 308834 57170000:57180000 0 chrY
  12. 308835 57180000:57190000 0 chrY
  13. 308836 57190000:57200000 0 chrY
  14. 308837 57200000:57210000 0 chrY
  15. 308838 57210000:57220000 0 chrY
  16. 308839 57220000:57230000 0 chrY
  17.  
  18. > levels(df$CHROM)
  19. [1] "chr1" "chr10" "chr11" "chr12" "chr13" "chr14" "chr15" "chr16" "chr17" "chr18" "chr19" "chr2" "chr20" "chr21" "chr22" "chr3" "chr4"
  20. [18] "chr5" "chr6" "chr7" "chr8" "chr9" "chrM" "chrX" "chrY"
  21.  
  22. # RE_ORDER CHROMS
  23. chrom_order<-c('chr1','chr2','chr3','chr4','chr5','chr6','chr7','chr8','chr9','chr10','chr11',
  24. 'chr12','chr13','chr14','chr15','chr16','chr17','chr18','chr19','chr20','chr21','chr22','chrX','chrM')
  25.  
  26. df<- df[match(chrom_order, df$CHROM),]
  27.  
  28. > head(df)
  29. POSITION MEANDEPTH CHROM
  30. 1 0:10000 0 chr1
  31. 128716 0:10000 0 chr2
  32. 169134 0:10000 0 chr3
  33. 188964 0:10000 0 chr4
  34. 207986 0:10000 0 chr5
  35. 226140 0:10000 0 chr6
  36.  
  37. library(dplyr)
  38. df %>%
  39. slice(match(CHROM, chrom_order))
Add Comment
Please, Sign In to add comment