# Repeat until empirical correlation is about 0.14 N = 747 X = rnorm(N) Y = -X + 6.5*rnorm(N) cc = cor(X,Y) cat(sprintf("%g\n",cc)) # Scale X & Y to correspond to the variables cited X = X-min(X); X = 80*X/max(X) Y = Y-min(Y); Y = 6*Y/max(Y); Y = Y-mean(Y); png(filename="FakeMusicData1.png", width=600, height=600) plot(X,Y, type="p", pch="x", col="red", xlab = "Empathy Quotient", ylab = "MUSIC Dimension \"Mellow\"") dev.off()