RAIN（草稿）

if (!require("BiocManager", quietly = TRUE))
  install.packages("BiocManager")

BiocManager::install("rain")
library(rain)
browseVignettes("rain")
help("browseVignettes")

### R code from vignette source 'rain.Rnw'

###################################################
### code chunk number 1: style-Sweave
###################################################
BiocManager::install("BiocStyle")
BiocStyle::latex()


###################################################
### code chunk number 2: rain.Rnw:35-43
###################################################
set.seed(123)
times <- c(1: 24) * 2
sin <- 1 + 0.5 * sin(times / 24 * 2 * pi) + rnorm(24, 0, 0.3)
saw <- rep(13:24 / 18 , 2) + rnorm(24, 0, 0.3)
measure <- cbind(sin, saw)
require('lattice')
xyplot(t(measure)~rep(times, each=2) | c('sin', 'saw'),
       layout = c(1, 2), type = 'o', xlab = 'time', ylab = 'value', cex.lab = 0.6)


###################################################
### code chunk number 3: rain.Rnw:49-55
###################################################
require(rain)
rainresult <- rain(measure, period=24, 
                   deltat=2, peak.border=c(0.1,0.9),
                   verbose=FALSE
)
rainresult