3sp_sim.Rdump

N <- 10
S <- 3
y <- 
structure(c(20.6660644730747, 38.9571084611488, 56.5619403130382, 6.4059324634517, 4.29931160236809, 3.81825520764806,
3.58332230847073, 3.56196269004428, 3.56196269004428, 3.56196269004428, 9.59863583423739, 18.8595349749282,
27.8737279344202, 50.9363800001624, 63.8989259358062, 18.668875206334, 14.4903503859552, 13.9870103346858,
13.9870103344666, 13.9870103344666, 9.60662147720212, 21.7509909684246, 34.7078183237834, 49.3240608787249,
65.1347019616715, 97.5128695860177, 121.926327305574, 142.45102697527, 152.007583839223, 152.007583839223),
.Dim = c(10, 3))
T <- 100
ts <- 
c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29,
30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58,
59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87,
88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100)
J <- 2
R <- 
structure(c(60, 60, 60, 60, 60, 60, 60, 60, 60, 60, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100),
.Dim = c(10, 2))
mu_true <- 
c(0.5, 0.75, 1)