inst/cwt.md

In nathaniel-mahieu/centWaveP: CWT Based Chromatographic Peak Detection

Wavelet Scale Correspondance

library(centWaveP)

for (i in seq(1, 50, 5)) {
  wav = centWaveP:::return.wavelet('mexh') %>% { centWaveP:::scale.wavelet(i, .)$y }
  plot(wav)
  cat(sum(wav > 0), i, "\n")
  }

## 2 1

## 20 6

## 26 11

## 32 16

## 42 21

## 52 26

## 62 31

## 72 36

## 82 41

## 92 46

for (i in seq(1, 50, 5)) {
  wav = centWaveP:::return.wavelet('nmexh') %>% { centWaveP:::scale.wavelet(i, .)$y }
  plot(wav)
  cat(sum(wav < 0), i, "\n")
  }

## 2 1

## 20 6

## 26 11

## 32 16

## 42 21

## 52 26

## 62 31

## 72 36

## 82 41

## 92 46

eic = { dnorm(seq(-6, 6, by =0.1)) } %>% { ./max(.) }
plot(eic)

for (i in seq(1, 50, 5)) {
  wav = cwt(eic, i, 'mexh')
  print(plot(wav, main = i))
  cat(sum(wav > 0), i, "\n")
}

## NULL
## 21 1

## NULL
## 23 6

## NULL
## 29 11

## NULL
## 38 16

## NULL
## 47 21

## NULL
## 56 26

## NULL
## 65 31

## NULL
## 75 36

## NULL
## 84 41

## NULL
## 94 46

eic = { dnorm(seq(-6, 6, by =0.1)) + dnorm(seq(-6, 6, by =0.1), mean = 3) } %>% { ./max(.) }
plot(eic)

for (i in seq(1, 50, 5)) {
  wav = cwt(eic, i, 'nmexh')
  print(plot(wav, main = i))
  cat(sum(wav < 0), i, "\n")
}

## NULL
## 35 1

## NULL
## 38 6

## NULL
## 53 11

## NULL
## 59 16

## NULL
## 61 21

## NULL
## 67 26

## NULL
## 74 31

## NULL
## 82 36

## NULL
## 90 41

## NULL
## 95 46

nathaniel-mahieu/centWaveP documentation built on May 23, 2019, 12:19 p.m.