inst/scripts/Nonlinear/dopler.R

In gbonte/gbcode: Code from the handbook "Statistical foundations of machine learning"

dopler<-function(x){
  20*sqrt(x*(1-x))*sin(2*pi*1.05/(x+0.05))
}


dataset.dopler<-function(N,sigma=1){
  set.seed(0)
  x<-sort(runif(N,min=0.12,max=1))
  y<-dopler(x)+rnorm(N,sd=sigma)
  x.ts<-sort(runif(N,min=0.12,max=1))
  y.ts<-dopler(x.ts)
  list(x=x,y=y,x.ts=x.ts,y.ts=y.ts)
}