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
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
|
\name{BabySS}
\docType{data}
\alias{BabySS}
\title{Physiological data time series.}
\description{
Two linked medical time series containing 2048 observations sampled every 16 seconds recorded from 21:17:59 to 06:27:18. Both these time series were recorded from the same 66 day old infant by Prof. Peter Fleming, Dr Andrew Sawczenko and Jeanine Young of the Institute of Child Health, Royal Hospital for Sick Children, Bristol. \code{BabyECG}, is a record of the infant's heart rate (in beats per minute). \code{BabySS} is a record of the infant's sleep state on a scale of 1 to 4 as determined by a trained expert monitoring EEG (brain) and EOG (eye-movement). The sleep state codes are 1=quiet sleep, 2=between quiet and active sleep, 3=active sleep, 4=awake.
}
\format{
The \code{BabyECG} time series is a nice examples of a non-stationary time series whose spectral (time-scale) properties vary over time. The function \code{\link{ewspec}} can be used to anaylse this time series to inspect the variation in the power of the series over time and scales.
The \code{BabySS} time series is a useful independent time series that can be associated with changing power in the \code{BabyECG} series. See the discussion in Nason, von Sachs and Kroisandt.
}
\source{Institute of Child Health, Royal Hospital for Sick Children, Bristol.}
\section{RELEASE}{
Version 3.9 Copyright Guy Nason 1998 }
\references{
Nason, G.P., von Sachs, R. and Kroisandt, G. (1998). Wavelet processes and adaptive estimation of the evolutionary wavelet spectrum. \emph{Technical Report}, Department of Mathematics University of Bristol/ Fachbereich Mathematik, Kaiserslautern.
}
\section{SEE ALSO}{
\code{\link{ewspec}}
}
\examples{
data(BabyECG)
data(BabySS)
#
# Plot the BabyECG data with BabySS overlaid
#
# Note the following code does some clever scaling to get the two
# time series overlaid.
#
myhrs <- c(22, 23, 24, 25, 26, 27, 28, 29, 30)
mylab <- c("22", "23", "00", "01", "02", "03", "04", "05", "06")
initsecs <- 59 + 60 * (17 + 60 * 21)
mysecs <- (myhrs * 3600)
secsat <- (mysecs - initsecs)/16
mxy <- max(BabyECG)
mny <- min(BabyECG)
ro <- range(BabySS)
no <- ((mxy - mny) * (BabySS - ro[1]))/(ro[2] - ro[1]) + mny
rc <- 0:4
nc <- ((mxy - mny) * (rc - ro[1]))/(ro[2] - ro[1]) + mny
\dontrun{plot(1:length(BabyECG), BabyECG, xaxt = "n", type = "l", xlab =
"Time (hours)", ylab = "Heart rate (beats per minute)")}
\dontrun{lines(1:length(BabyECG), no, lty = 3)}
\dontrun{axis(1, at = secsat, labels = mylab)}
\dontrun{axis(4, at = nc, labels = as.character(rc))}
#
# Sleep state is the right hand axis
#
#
}
\keyword{datasets}
\author{G P Nason}
|
