Guest User

Untitled

a guest
Nov 18th, 2018
79
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.82 KB | None | 0 0
  1. US1ORLA0076 <- read_csv("US1ORLA0076_cf.csv")
  2. head(US1ORLA0076)
  3.  
  4. water_date PRCP prcp_log
  5. <date> <dbl> <dbl>
  6. 1 2006-12-22 0.09 0.0899
  7. 2 2006-12-23 0.75 0.693
  8. 3 2006-12-24 1.63 1.26
  9. 4 2006-12-25 0.06 0.0600
  10. 5 2006-12-26 0.36 0.353
  11. 6 2006-12-27 0.63 0.594
  12.  
  13. library(mFilter)
  14.  
  15. US1ORLA0076_cffilter <- cffilter(US1ORLA0076$prcp_log,pl=180,pu=365,root=FALSE,drift=FALSE,
  16. type=c("asymmetric"),
  17. nfix=NULL,theta=1)
  18.  
  19. head(US1ORLA0076_cffilter$trend)
  20. [,1]
  21. [1,] 0.1482724
  22. [2,] 0.7501137
  23. [3,] 1.3202868
  24. [4,] 0.1139883
  25. [5,] 0.4051551
  26. [6,] 0.6453462
  27.  
  28. head(US1ORLA0076_cffilter$cycle)
  29. [,1]
  30. [1,] -0.05839342
  31. [2,] -0.05696651
  32. [3,] -0.05550995
  33. [4,] -0.05402422
  34. [5,] -0.05250982
  35. [6,] -0.05096727
  36.  
  37. plot(US1ORLA0076_cffilter)
Add Comment
Please, Sign In to add comment