Advertisement
Guest User

Untitled

a guest
Apr 17th, 2015
212
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.58 KB | None | 0 0
  1. > jj
  2. Qtr1 Qtr2 Qtr3 Qtr4
  3. 1960 0.71 0.63 0.85 0.44
  4. 1961 0.61 0.69 0.92 0.55
  5. 1962 0.72 0.77 0.92 0.60
  6. 1963 0.83 0.80 1.00 0.77
  7. 1964 0.92 1.00 1.24 1.00
  8. 1965 1.16 1.30 1.45 1.25
  9. 1966 1.26 1.38 1.86 1.56
  10. 1967 1.53 1.59 1.83 1.86
  11. 1968 1.53 2.07 2.34 2.25
  12. 1969 2.16 2.43 2.70 2.25
  13. 1970 2.79 3.42 3.69 3.60
  14. 1971 3.60 4.32 4.32 4.05
  15. 1972 4.86 5.04 5.04 4.41
  16. 1973 5.58 5.85 6.57 5.31
  17. 1974 6.03 6.39 6.93 5.85
  18. 1975 6.93 7.74 7.83 6.12
  19. 1976 7.74 8.91 8.28 6.84
  20. 1977 9.54 10.26 9.54 8.73
  21. 1978 11.88 12.06 12.15 8.91
  22. 1979 14.04 12.96 14.85 9.99
  23. 1980 16.20 14.67 16.02 11.61
  24.  
  25. k = c(.5, 1, 1, 1, .5)
  26.  
  27. fjj1 = filter(jj, sides=2,k)
  28. max_threshold = fjj1 * 1.2
  29. min_threshold = fjj1 * .8
  30.  
  31. min_outliers <- jj < min_threshold
  32. max_outlieres <- jj > max_threshold
  33.  
  34. min_outliers[min_outliers==FALSE] <- NA
  35. max_outliers[max_outliers==FALSE] <-- NA
  36.  
  37. plot(jj)
  38. lines(max_threshold, col="red", lty="dashed")
  39. lines(min_threshold, col="green", lty="dashed")
  40. lines(fjj1, col="blue", lty="dashed")
  41.  
  42. > min_outliers
  43. Qtr1 Qtr2 Qtr3 Qtr4
  44. 1960 NA NA NA TRUE
  45. 1961 NA NA NA TRUE
  46. 1962 NA NA NA TRUE
  47. 1963 NA NA NA NA
  48. 1964 NA NA NA NA
  49. 1965 NA NA NA NA
  50. 1966 NA NA NA NA
  51. 1967 NA NA NA NA
  52. 1968 NA NA NA NA
  53. 1969 NA NA NA NA
  54. 1970 NA NA NA NA
  55. 1971 NA NA NA NA
  56. 1972 NA NA NA NA
  57. 1973 NA NA NA NA
  58. 1974 NA NA NA NA
  59. 1975 NA NA NA NA
  60. 1976 NA NA NA TRUE
  61. 1977 NA NA NA NA
  62. 1978 NA NA NA TRUE
  63. 1979 NA NA NA TRUE
  64. 1980 NA NA NA NA
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement