-
Notifications
You must be signed in to change notification settings - Fork 92
/
Copy pathTideDataR.Rmd
49 lines (39 loc) · 972 Bytes
/
TideDataR.Rmd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
---
title: "TideDataR"
output: github_document
---
Prepare tide data in R.
Read data.
```{r}
library(wrapr)
files <- list.files('tide_data',
pattern='.*.csv') %.>%
sort(.)
dat_list <- lapply(
files,
function(fi) {
read.csv(paste0('tide_data/', fi),
strip.white = TRUE,
na.strings = c("", "NA", "-"),
stringsAsFactors = FALSE)
})
tides = do.call(rbind, dat_list)
tides['tide_feet'] <- tides['Verified..ft.']
missing_v = is.na(tides['tide_feet'])
tides['tide_feet'][missing_v] <- tides['Preliminary..ft.'][missing_v]
summary(tides['tide_feet'])
```
```{r}
tides['dt'] <- as.POSIXct(paste(tides[['Date']], tides[['Time..GMT.']]),
tz = "UTC")
head(tides)
```
```{r}
deltas <- difftime(tides[['dt']][-1],
tides[['dt']][-length(tides[['dt']])],
units='secs')
summary(as.numeric(deltas))
```
```{r}
saveRDS(tides, 'tides.RDS')
```