Nothing
# vim:textwidth=80:expandtab:shiftwidth=4:softtabstop=4
library(oce)
test_that("as.lisst()", {
set.seed(1333334L)
t <- seq(0, 6, 1/15) * 3600 + as.POSIXct("2012-01-01 00:00:00", tz="UTC")
n <- length(t)
pre <- 5 + sin(as.numeric(t - t[1]) / 12.4 / 3600 * 2 * pi) + rnorm(n, sd=0.01)
dpdt <- c(0, diff(pre))
tem <- 10 + 5 * sin(as.numeric(t - t[1]) / 24 / 3600 * 2 * pi) + cumsum(rnorm(n, sd=0.2))
C <- (dpdt + rnorm(n, sd=0.1) + cumsum(rnorm(n, sd=0.5)))^2 * 2
sd <- rep(1, length.out=32) + (1:32) / 100
data <- matrix(nrow=n, ncol=42)
for (i in 1:32) {
fake <- abs(C * (1 + i / 5) + cumsum(rnorm(n, sd=sd[i]))) / 100
data[, i] <- fake
}
data[, 33] <- rep(0, n) # lts
data[, 34] <- rep(4, n) + 0.01 * cumsum(rnorm(n, 0.05)) # voltage
data[, 35] <- rep(0.07, n) # aux
data[, 36] <- runif(n, 3.9, 4.1) # lrs
data[, 37] <- pre
data[, 38] <-tem
tt <- as.POSIXlt(t)
data[, 39] <- 100 * tt$yday + tt$hour
data[, 40] <- 100 * tt$min + tt$sec
data[, 41] <- abs((pre - min(pre)) / diff(range(pre)) + cumsum(rnorm(n, sd=0.05))) # transmission
data[, 42] <- 40 - 20*data[, 41] # beam
lisst <- as.lisst(data, filename="(constructed)", year=2012, "UTC")
expect_equal(lisst[["beam"]], 40-20*data[, 41])
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.