Guest User

Untitled

a guest
Dec 12th, 2018
74
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.78 KB | None | 0 0
  1. $gw6.00033
  2. X Sample_Name Chr_No Copy_No Frequence
  3. 1 gw6.00033 chr1: cn=0 1
  4. 12 gw6.00033 chr2: cn=0 1
  5. 16 gw6.00033 chr4: cn=0 1
  6. 20 gw6.00033 chr6: cn=0 1
  7.  
  8. $gw6.0006
  9. X Sample_Name Chr_No Copy_No Frequence
  10. 26 gw6.0006 chr1: cn=0 1
  11. 29 gw6.0006 chr10 cn=0 3
  12. 31 gw6.0006 chr11 cn=0 2
  13. 34 gw6.0006 chr13 cn=0 1
  14. 37 gw6.0006 chr15 cn=0 1
  15. 38 gw6.0006 chr16 cn=0 1
  16. 41 gw6.0006 chr2: cn=0 1
  17. 47 gw6.0006 chr3: cn=0 1
  18. 57 gw6.0006 chr8: cn=0 2
  19.  
  20. sp <- split(cn0, cn0$sample_name)
  21. N <- 22
  22. sp <- lapply(sp, function(DF){
  23. M <- as.numeric(sub("[^[:digit:]]+([[:digit:]]+)[^[:digit:]]*", "\1", DF$Chr_no))
  24. Chr_no <- sub("[[:digit:]]+[^[:digit:]]*$", "", DF$Chr_no)
  25. Chr_no <- paste0(Chr_no, 1:N)
  26. Chr_no <- ifelse(nchar(Chr_no) == 4, paste0(Chr_no, ":"), Chr_no)
  27. sample_name <- rep(DF$sample_name[1], length(Chr_no))
  28. res <- data.frame(sample_name, Chr_no)
  29. res$Frequence <- 0
  30. res$Frequence[1:N %in% M] <- DF$Frequence
  31. res
  32. })
  33. abc <- as.data.frame(t(sapply(sp, '[[', 'Frequence')))
  34. names(abc) <- sp[[1]]$Chr_no
  35.  
  36. $gw6.00033
  37. Sample_Name Chr_No Frequence
  38. 1 gw6.00033 chr1: 1
  39. 2 gw6.00033 chr2: 1
  40. 3 gw6.00033 chr3: 0
  41. 4 gw6.00033 chr4: 1
  42. 5 gw6.00033 chr5: 0
  43. 6 gw6.00033 chr6: 1
  44. 7 gw6.00033 chr7: 0
  45. 8 gw6.00033 chr8: 0
  46. 9 gw6.00033 chr9: 0
  47.  
  48.  
  49. $gw6.0006
  50. Sample_Name Chr_No Frequence
  51. 1 gw6.0006 chr1: 1
  52. 2 gw6.0006 chr2: 3?
  53. 3 gw6.0006 chr3: 2?
  54. 4 gw6.0006 chr4: 0
  55. 5 gw6.0006 chr5: 0
  56. 6 gw6.0006 chr6: 0
  57. 7 gw6.0006 chr7: 0
  58. 8 gw6.0006 chr8: 1?
Add Comment
Please, Sign In to add comment