Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- US1ORLA0076 <- read_csv("US1ORLA0076_cf.csv")
- head(US1ORLA0076)
- water_date PRCP prcp_log
- <date> <dbl> <dbl>
- 1 2006-12-22 0.09 0.0899
- 2 2006-12-23 0.75 0.693
- 3 2006-12-24 1.63 1.26
- 4 2006-12-25 0.06 0.0600
- 5 2006-12-26 0.36 0.353
- 6 2006-12-27 0.63 0.594
- library(mFilter)
- US1ORLA0076_cffilter <- cffilter(US1ORLA0076$prcp_log,pl=180,pu=365,root=FALSE,drift=FALSE,
- type=c("asymmetric"),
- nfix=NULL,theta=1)
- head(US1ORLA0076_cffilter$trend)
- [,1]
- [1,] 0.1482724
- [2,] 0.7501137
- [3,] 1.3202868
- [4,] 0.1139883
- [5,] 0.4051551
- [6,] 0.6453462
- head(US1ORLA0076_cffilter$cycle)
- [,1]
- [1,] -0.05839342
- [2,] -0.05696651
- [3,] -0.05550995
- [4,] -0.05402422
- [5,] -0.05250982
- [6,] -0.05096727
- plot(US1ORLA0076_cffilter)
Add Comment
Please, Sign In to add comment