rpart
dat = data.frame(x1=rnorm(1000,mean=0,sd =3 ),
x2=c(rnorm(500,mean=2,sd=.5),rnorm(500,mean=-2,sd=.5)),
y=rep(0:1,times=c(500,500)))
plot(dat$x1,dat$x2)
rotate = matrix(c(cos(pi/4),sin(pi/4),-sin(pi/4),cos(pi/4)),byrow=F,nc=2)
m = dat[,1:2]
m = as.matrix(m)
mt = m %*% rotate
d = cbind(mt,rep(0:1,times=c(500,500)))
d = as.data.frame(d)
d$V3 = as.factor(d$V3)
library(rpart)
fit = rpart(V3 ~ . , data = d)