data("GH")
sel <- which(GH$COMPS=="V")
amp <- list()
dees <- list()
stamps <- list()
for( i in 1:3)
{
amp[[i]] <- GH$JSTR[[sel[i]]]
dees[i] <- GH$dt[sel[i]]
stamps[i] <- paste(GH$STNS[sel[i]], GH$COMPS[sel[i]])
}
a <- list(y=amp, dt=dees, stamps=stamps)
f1 <- 0.1
f2 <- floor(0.33*(1/a$dt[[1]]))
speccol <- c('red', 'blue', 'purple')
MTM.drive(a, f1, f2, COL=speccol, PLOT=TRUE)
Run the code above in your browser using DataLab