Advertisement
Guest User

Untitled

a guest
Dec 9th, 2016
80
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 2.07 KB | None | 0 0
  1. PERSON_ID EVENT
  2. 1 10000000001 LST_Mammo_EPI
  3. 2 10000000001 EPI
  4. 3 10000000002 BC1
  5. 4 10000000002 R_B
  6. 5 10000000002 BC2
  7. 6 10000000002 DEATH
  8. 7 10000000002 EPI
  9. 8 10000000004 FST_Mammo_FUP1
  10. 9 10000000004 OV1_FUP1
  11. 10 10000000004 BC1
  12. 11 10000000004 LST_Mammo_EPI
  13. 12 10000000004 EPI
  14. 13 10000000004 OV2_FUP1
  15. 14 10000000004 LST_Mammo_FUP1
  16. 15 10000000004 FUP1
  17. 16 10000000004 BC2
  18. 17 10000000004 FUP2
  19. 18 10000000005 LST_Mammo_EPI
  20. 19 10000000005 EPI
  21. 20 10000000007 BC1
  22. 21 10000000007 BC2
  23. 22 10000000007 DEATH
  24. 23 10000000010 DEATH
  25. 24 10000000826 FST_Mammo_FUP1
  26. 25 10000000826 BC1
  27. 26 10000000826 L_B
  28. 27 10000000826 LST_Mammo_EPI
  29. 28 10000000826 LST_Mammo_FUP1
  30. 29 10000000826 BC2
  31. 30 10000000826 R_B
  32. 31 10000000826 EPI
  33. 32 10000000826 OV1_FUP1
  34. 33 10000000826 OV2_FUP1
  35. 34 10000000826 OV2_FUP2
  36. 35 10000000826 FUP1
  37. 36 10000000826 FUP2
  38. 37 10000000827 DEATH
  39. 38 10000000830 BC1
  40. 39 10000000830 DEATH
  41. 40 10000000844 LST_Mammo_EPI
  42.  
  43. node <- function(proxis, fulle, database){
  44.  
  45. result <- NULL
  46.  
  47. for (ID in unique(database$PERSON_ID)){
  48. person <- subset(database, database$PERSON_ID == ID)
  49. a <- which(person $EVENT == fulle)
  50.  
  51. if (length(a) == 0 ){
  52. next} else {
  53.  
  54. person <- person[c(1:a),]
  55. b <- which(person$EVENT == proxis)
  56.  
  57. if (length(b) == 0) {
  58. next} else{
  59.  
  60. person <- person[c(b:a),]
  61. a <- which(person$EVENT == "OV")
  62. b <- which(person$EVENT == "BC1")
  63. c <- which(person$EVENT == "BC2")
  64. if (min(a,b,c) != 0) {next
  65. }else{
  66. result <- rbind(result,person)
  67.  
  68. }}}}
  69.  
  70. result
  71. }
  72.  
  73. hello <- node("BC1","BC2",test)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement