n <- 1000
peou <- rnorm(n)
peou1 <- peou * 1.1 + rnorm(n, 0, 0.25)
peou2 <- peou * 0.9 + rnorm(n, 0, 0.15)
peou3 <- peou * 0.75 + rnorm(n, 0, 0.20)
peou.copy <- rnorm(n)
peou1.copy <- peou.copy * 1.1 + rnorm(n, 0, 0.25)
peou2.copy <- peou.copy * 0.9 + rnorm(n, 0, 0.15)
peou3.copy <- peou.copy * 0.75 + rnorm(n, 0, 0.20)
my.data <- data.frame(cbind(peou1, peou2, peou3, peou1.copy, peou2.copy, peou3.copy))
write.csv(my.data, '~/Desktop/dataset12.csv', row.names=F)

