set.seed(1234) signal = c(-5.65, -5.56, -2.62, -1.20, -1.01, -0.90, -0.15, 1.65, 1.94, 3.57) n = 25 noise = rnorm(n, 0, 1) Y = c(signal, noise) dput(list(Y=Y), "DataMultipleComparisons.txt") beta = rep(0, times=n) gamma = rep(0, times=n) dput(list(beta=beta), "beta.txt") dput(list(gamma=gamma), "gamma.txt")