n <- 500
eou <- rnorm(n, 0, 1)
eou1 <- .9*eou + rnorm(n, 0, 0.1)
eou2 <- .95*eou + rnorm(n, 0, 0.08)
eou3 <- .85*eou + rnorm(n, 0, 0.12)
use <- eou + rnorm(n, 0, 0.1)
use1 <- .95*use + rnorm(n, 0, 0.08)
use2 <- .9*use + rnorm(n, 0, 0.1)
use3 <- .85*use + rnorm(n, 0, 0.12)
int <- eou + use + rnorm(n, 0, 0.05)
int1 <- int + rnorm(n, 0, 0.01)
int2 <- int + rnorm(n, 0, 0.03)
int3 <- int + rnorm(n, 0, 0.02)

data24 <- data.frame(cbind(eou1, eou2, eou3, use1, use2, use3, int1, int2, int3))

names(data24)=c('eou1', 'eou2', 'eou3', 'use1', 'use2', 'use3', 'int1', 'int2', 'int3')

write.csv(data24, 'dataset24.csv', row.names=F)
