##Simple bivariate normal data
set.seed( 1 )
n = 15
d = 2
props=c( 0.6, 0.4 )
x <- matrix( rnorm( n*d ), ncol = d )
shiftvec <- ifelse( runif( n ) > props[ 1 ], 0, 1)
x[,1] <- x[,1] + shiftvec
interactive2D( x, shiftvec )
Run the code above in your browser using DataLab