Site Tools

lab:projects:13applying_student_t_distribution_on_observed_variables_or_latent_variables:t-distribution_on_random_effects
## A robust model with student t-distribution on latent factors

model{
# Model specification for linear growth curve model
for (i in 1:500){
LS[i,1:2]~dmnorm(muLS[i,1:2], Inv_cov[1:2,1:2])

muLS[i,1]<-bL[1]
muLS[i,2]<-bS[1]

for (t in 1:5){
y[i, t] ~ dt(muY[i,t], Inv_Sig_e2, d)
muY[i,t]<-LS[i,1]+LS[i,2]*(t-1)
}
}

#Priors for model parameter
for (i in 1:1){
bL[i] ~ dnorm(0, 1.0E-6)
bS[i] ~ dnorm(0, 1.0E-6)
}

d~dunif(0,100)

Inv_cov[1:2,1:2]~dwish(R[1:2,1:2], 2)
R[1,1]<-1
R[2,2]<-1
R[2,1]<-R[1,2]
R[1,2]<-0
Inv_Sig_e2 ~ dgamma(.001, .001)

Cov[1:2,1:2]<-inverse(Inv_cov[1:2,1:2])
para[1]<-bL[1]
para[2]<-bS[1]
para[3] <- Cov[1,1]
para[4] <- Cov[2,2]
para[5] <- Cov[1,2]
para[6] <- 1/Inv_Sig_e2
para[7] <- d
}

# The (naive) starting values for model parameters.
list(bL=c(6), bS = c(1), Inv_Sig_e2 = 1, d=5,
Inv_cov= structure(.Data = c(1,0,0,1),.Dim=c(2,2)))

list(  y = structure(.Data = c( 3.858, 5.526, 5.928, 5.549, 5.959,
6.145, 7.168, 7.515, 7.782, 6.608, 5.49, 3.841, 4.849, 6.546, 6.088, 8.598,
6.752, 5.042, 6.676, 6.276, 6.29, 7.012, 6.517, 7.526, 6.836, 9.766, 7.234,
7.58, 7.236, 5.101, 5.03, 9.004, 5.92, 8.211, 7.335, 7.143, 5.207, 5.628,
6.238, 7.082, 7.528, 6.818, 7.632, 7.115, 7.149, 4.986, 3.791, 4.078, 6.215,
6.708, 6.255, 7.753, 8.143, 7.751, 9.903, 4.631, 5.725, 6.794, 7.659, 6.889,
5.915, 6.219, 6.23, 6.51, 6.935, 8.473, 7.49, 7.305, 8.928, 6.891, 7.982,
10.35, 8.338, 11.28, 7.794, 6.437, 6.721, 6.274, 5.586, 5.742, 4.883, 6.912,
6.542, 7.21, 6.384, 4.728, 4.178, 4.338, 4.377, 5.868, 7.364, 5.768, 6.873,
8.991, 5.385, 7.02, 5.414, 8.139, 4.959, 7.218, 3.837, 3.075, 4.929, 4.072,
5.971, 6.033, 5.9, 6.284, 6.782, 6.731, 5.462, 5.471, 6.167, 5.068, 5.607,
9.644, 7.907, 6.777, 9.155, 7.659, 8.442, 8.662, 10.46, 10.22, 8.819, 3.543,
3.497, 3.459, 4.316, 4.001, 2.35, 3.883, 4.26, 7.685, 5.312, 3.598, 3.609,
4.34, 3.881, 4.415, 5.241, 3.871, 5.547, 6.131, 4.7, 3.402, 3.974, 3.074,
2.724, 2.515, 4.339, 6.848, 6.64, 5.761, 7.462, 4.496, 5.86, 6.184, 6.007,
5.89, 3.587, 5, 7.817, 4.921, 5.257, 8.484, 8.059, 7.903, 7.571, 8.874, 6.587,
6.675, 7.783, 7.217, 7.129, 4.78, 6.19, 6.19, 6.736, 6.505, 5.912, 6.054,
6.312, 6.659, 6.172, 7.757, 8.453, 8.181, 9.22, 8.036, 7.041, 4.99, 9.161,
9.405, 9.186, 4.741, 5.859, 6.821, 5.454, 5.182, 4.792, 4.216, 2.561, 3.352,
4.23, 5.821, 7.994, 7.625, 5.762, 6.814, 4.774, 6.864, 6.265, 7.244, 7.412,
6.93, 6.784, 6.093, 8.518, 7.354, 2.74, 3.266, 3.163, 4.272, 5.121, 5.305,
5.616, 5.122, 6.788, 5.966, 6.127, 4.747, 4.347, 5.751, 5.197, 7.186, 6.301,
8.075, 8.347, 9.619, 9.574, 9.136, 8.743, 10.27, 10.31, 6.136, 8.286, 7.811,
7.021, 7.605, 3.545, 5.588, 6.039, 6.079, 5.598, 6.972, 7.926, 8.085, 6.146,
7.82, 7.019, 5.841, 5.136, 8.295, 7.466, 5.378, 4.047, 5.192, 6.051, 4.945,
5.816, 6.266, 6.135, 5.233, 5.969, 6.665, 7.223, 7.001, 6.53, 7.816, 5.457,
5.55, 6.43, 8.348, 5.774, 6.436, 8.449, 11.37, 8.387, 10.87, 4.708, 4.646,
4.207, 4.622, 4.038, 3.047, 4.645, 3.486, 4.208, 3.337, 3.83, 4.869, 5.69,
7.626, 6.462, 5.479, 5.305, 5.059, 5.598, 6.379, 6.516, 7.515, 7.189, 6.31,
8.111, 3.843, 5.014, 4.231, 3.132, 4.704, 5.862, 4.172, 3.377, 3.43, 6.553,
4.833, 4.803, 5.282, 4.174, 5.482, 6.28, 6.613, 5.798, 7.076, 9.189, 8.122,
6.383, 5.021, 10.37, 8.611, 7.393, 8.754, 9.178, 8.121, 9.682, 5.431, 5.121,
4.955, 10.71, 6.539, 8.276, 8.386, 8.898, 8.46, 8.8, 5.82, 4.31, 6.659, 4.874,
6.432, 6.133, 7.868, 5.894, 7.995, 8.315, 7.06, 7.219, 6.702, 8.14, 8.668,
9.631, 10.5, 9.688, 8.511, 7.708, 8.49, 5.18, 7.686, 7.96, 6.384, 2.736, 4.832,
4.775, 4.776, 4.603, 5.586, 6.939, 8.266, 9.077, 7.647, 4.91, 7.741, 7.441,
3.595, 9.304, 7.14, 6.894, 6.698, 7.391, 6.149, 3.422, 5.974, 7.718, 6.164,
8.259, 2.298, 2.997, 3.874, 3.394, 4.371, 7.731, 7.845, 7.959, 9.907, 10.06,
5.344, 3.697, 4.796, 4.879, 5.389, 4.252, 5.026, 4.925, 6.91, 6.137, 4.043,
5.64, 3.18, 3.971, 4.707, 3.714, 2.588, 3.186, 3.192, 5.178, 7.671, 8.271,
8.206, 9.515, 9.381, 5.99, 6.282, 7.158, 6.938, 7.265, 5.461, 5.948, 5.852,
5.252, 4.557, 3.148, 8.62, 8.952, 8.049, 12.59, 4.058, 5.811, 5.568, 6.923,
6.91, 3.878, 2.56, 2.938, 3.936, 5.695, 7.435, 7.378, 8.228, 6.519, 7.651,
4.691, 7.011, 6.903, 7.321, 7.476, 9.136, 6.657, 6.069, 7.784, 5.122, 3.568,
3.413, 5.325, 4.429, 5.445, 5.809, 3.066, 6.069, 7.719, 11.72, 4.975, 4.65,
5.273, 6.268, 5.279, 5.72, 7.334, 8.314, 8.378, 8.187, 2.909, 3.172, 4.06,
3.051, 3.321, 5.622, 5.177, 5.148, 4.722, 5.708, 6.748, 8.377, 8.934, 7.071,
8.501, 5.715, 7.211, 6.248, 6.206, 6.731, 7.092, 7.295, 7.064, 9.578, 9.745,
2.459, 4.676, 3.097, 2.141, 3.349, 7.149, 7.005, 7.579, 8.31, 7.709, 6.918,
7.441, 7.31, 8.055, 8.343, 7.289, 8.046, 7.495, 8.414, 9.223, 3.535, 7.578,
6.331, 6.677, 5.599, 4.993, 10.42, 8.805, 2.372, 5.703, 8.953, 8.326, 10.94,
10.25, 11.32, 7.115, 7.018, 7.802, 8.159, 8.154, 3.998, 4.206, 4.368, 3.877,
4.729, 4.907, 6.061, 5.586, 6.867, 7.057, 8.165, 8.249, 8.649, 8.292, 7.351,
5.962, 6.35, 8.594, 7.32, 8.778, 7.055, 7.004, 8.283, 10.18, 7.879, 7.203,
5.372, 6.231, 9.611, 8.102, 7.129, 6.362, 5.74, 6.86, 7.279, 7.564, 7.35,
8.103, 7.221, 8.232, 6.883, 5.625, 6.027, 3.921, 6.944, 5.889, 5.264, 5.778,
5.709, 5.08, 3.679, 3.866, 5.813, 5.779, 5.605, 2.698, 4.354, 4.787, 5.329,
5.995, 8.852, 6.852, 6.928, 8.107, 7.997, 5.748, 5.986, 6.445, 5.705, 7.877,
7.406, 7.283, 8.253, 7.843, 9.136, 6.462, 6.208, 7.677, 8.853, 7.911, 10.25,
11.07, 9.652, 10.45, 10.53, 3.553, 4.788, 4.568, 4.253, 3.985, 6.441, 6.189,
6.645, 7.36, 7.929, 6.033, 6.341, 6.687, 6.874, 7.269, 7.019, 5.905, 6.481,
6.425, 5.661, 4.977, 6.278, 4.409, 6.872, 7.055, 3.625, 6.144, 6.416, 8.879,
7.34, 5.431, 4.76, 7.672, 4.94, 7.672, 8.316, 6.917, 9.166, 9.242, 9.474,
9.339, 8.118, 8.63, 9.606, 7.806, 6.056, 6.415, 6.997, 6.244, 6.797, 6.914,
5.536, 5.084, 6.664, 6.7, 9.792, 6.14, 2.188, 4.341, 6.845, 6.907, 6.237,
5.698, 6.484, 6.014, 6.512, 6.486, 6.482, 6.685, 7.443, 5.217, 4.826, 5.037,
4.746, 5.685, 1.106, 4.26, 2.873, 4.983, 5.756, 6.63, 5.652, 6.019, 4.021,
9.453, 6.344, 7.434, 6.827, 7.456, 8.307, 3.887, 3.515, 5.23, 4.537, 8.933,
4.195, 5.317, 8.697, 6.69, 5.831, 6.572, 7.953, 7.744, 6.418, 7.836, 5.19,
6.109, 7.539, 7.685, 8.139, 8.646, 9.929, 7.055, 7.151, 9.289, 5.103, 6.189,
5.74, 6.185, 5.051, 5.325, 6.621, 6.513, 6.799, 7.513, 6.901, 6.206, 8.114,
7.709, 8.265, 7.779, 9.334, 9.701, 9.84, 9.654, 5.748, 6.178, 6.835, 6.82,
6.578, 5.204, 5.884, 4.681, 5.03, 7.197, 6.069, 5.962, 6.212, 6.593, 6.852,
7.643, 10.53, 10.06, 3.334, 6.845, 0.7468, 2.157, 5.238, 4.784, 2.393, 6.52,
8.792, 4.031, 8.875, 7.816, 7.293, 2.091, 6.839, 7.82, 5.211, 3.636, 2.51, 3.4,
5.576, 5.288, 5.943, 4.502, 7.686, 7.025, 8.057, 4.932, 4.922, 5.533, 5.906,
6.101, 4.776, 5.26, 6.856, 6.178, 6.536, 6.701, 4.981, 5.25, 6.936, 7.523,
6.688, 6.567, 6.534, 5.802, 5.39, 4.415, 1.889, 6.23, 4.699, 3.203, 7.37,
7.022, 7.616, 6.638, 7.494, 7.504, 9.343, 9.064, 9.62, 8.448, 7.638, 7.385,
7.921, 6.632, 7.943, 7.308, 9.011, 3.197, 10.16, 9.951, 3.567, 4.479, 3.725,
5.388, 4.497, 8.539, 8.011, 9.234, 8.666, 9.011, 4.847, 5.502, 4.276, 6.53,
6.246, 6.996, 5.593, 7.126, 7.168, 7.031, 5.811, 7.628, 4.573, 7.501, 9.406,
5.546, 4.771, 5.148, 6.296, 6.61, 5.33, 6.051, 3.882, 7.412, 7.714, 4.275,
7.856, 6.194, 7.263, 7.667, 7.583, 6.669, 7.933, 7.262, 9.659, 8.649, 7.617,
7.501, 7.945, 7.692, 7.191, 7.501, 6.088, 6.415, 7.323, 7.041, 7.571, 6.123,
6.572, 9.09, 5.847, 7.103, 6.147, 8.849, 7.557, 4.863, 6.518, 6.744, 7.131,
7.181, 5.003, 5.223, 6.035, 5.558, 4.581, 5.724, 6.566, 5.51, 6.64, 6.476,
7.217, 7.318, 8.152, 6.472, 7.917, 7.798, 7.712, 8.383, 8.072, 9.297, 6.824,
5.543, 6.688, 6.612, 7.724, 5.054, 3.254, 5.653, 5.761, 6.555, 3.445, 4.571,
6.94, 4.704, 4.474, 2.971, 6.101, 6.748, 2.58, 8.673, 3.565, 5.233, 4.836,
4.723, 4.659, 6.464, 5.616, 6.041, 6.571, 6.2, 7.029, 6.524, 6.506, 5.902,
7.439, 5.103, 5.216, 4.322, 4.005, 5.271, 3.06, 3.595, 5.257, 5.402, 5.483,
5.604, 6.658, 6.889, 7.751, 8.581, 6.303, 5.145, 6.126, 7.902, 7.508, 4.625,
5.176, 5.195, 6.028, 5.149, 4.207, 4.339, 5.175, 5.651, 4.742, 9.698, 10.22,
8.582, 8.921, 10.96, 5.987, 7.46, 5.859, 6.601, 6.076, 4.46, 5.906, 4.499,
4.294, 7.611, 7.498, 6.176, 8.385, 5.962, 6.228, 6.003, 6.386, 6.4, 7.856,
6.429, 8.249, 9.415, 6.852, 8.861, 10.48, 4.366, 6.506, 7.386, 7.964, 5.117,
6.157, 5.933, 5.221, 8.333, 6.185, 4.887, 5.499, 5.287, 6.52, 6.202, 4.431,
4.878, 4.954, 4.988, 5.621, 7.563, 4.398, 6.196, 4.21, 7.046, 6.433, 6.488,
5.908, 5.907, 7.149, 5.339, 6.036, 5.802, 6.856, 5.135, 4.265, 4.086, 4.109,
6.328, 6.094, 2.742, 2.881, 3.839, 3.607, 3.809, 6.679, 7.974, 9.397, 7.735,
6.652, 6.641, 6.743, 7.224, 8.819, 8.316, 8.86, 10.9, 10.63, 10.53, 11.28,
3.35, 3.895, 4.911, 5.638, 5.614, 6.402, 4.9, 7.421, 6.138, 7.295, 5.135,
6.975, 5.037, 5.843, 5.433, 6.292, 4.963, 5.077, 6.034, 5.323, 7.71, 8.061,
7.773, 8.294, 9.015, 5.684, 7.454, 6.329, 7.28, 7.125, 5.867, 5.299, 4.831,
5.876, 5.136, 4.996, 3.758, 6.307, 6.636, 6.819, 11.3, 10.82, 11.8, 10.89,
12.27, 5.336, 4.835, 5.329, 5.434, 5.786, 6.514, 5.436, 5.223, 7.793, 6.701,
8.172, 1.607, 3.695, 5.03, 8.036, 6.28, 5.113, 6.356, 7.76, 7.277, 6.542,
6.641, 8.128, 6.997, 7.69, 4.459, 3.138, 4.295, 3.736, 4.417, 6.185, 6.322,
6.363, 6.771, 6.148, 6.758, 6.792, 6.905, 7.007, 7.197, 6.692, 6.253, 5.814,
5.707, 5.828, 5.401, 6.322, 6.015, 6.311, 5.779, 3.902, 4.298, 2.586, 3.792,
3.493, 8.271, 5.771, 8.079, 15.39, 10.08, 5.96, 7.189, 7.711, 7.805, 9.98,
5.836, 5.311, 6.55, 6.008, 5.84, 6.985, 6.754, 8.616, 9.284, 8.862, 4.479,
7.127, 7.276, 8.765, 5.694, 5.7, 8.061, 6.615, 8.698, 9.284, 8.179, 7.283,
7.736, 8.029, 7.755, 5.5, 6.013, 6.19, 7.092, 6.801, 4.282, 3.932, 3.473, 3.46,
3.932, 2.676, 4.777, 3.37, 4.279, 4.445, 5.725, 4.851, 4.58, 5.98, 5.08, 6.357,
8.648, 7.236, 9.811, 10.55, 2.148, 7.451, 8.129, 7.699, 7.798, 7.136, 1.375,
5.271, 5.447, 7.036, 2.569, 2.309, 4.339, 3.422, 2.766, 4.873, 4.087, 3.648,
4.394, 5.228, 8.647, 8.29, 9.997, 8.958, 9.008, 6.909, 7.605, 8.312, 9.143,
6.064, 7.729, 8.921, 7.796, 3.932, 8.125, 3.858, 4.961, 4.744, 5.054, 5.323,
7.248, 5.479, 7.255, 6.75, 8.962, 1.998, 1.872, 2.445, 2.922, 4.106, 7.282,
7.158, 6.866, 8.458, 8.071, 11.03, 7.636, 9.018, 13.21, 12.56, 4.76, 5.863,
6.142, 6.196, 7.787, 6.623, 6.742, 9.097, 7.111, 7.033, 6.352, 7.176, 7.123,
8.197, 7.655, 6.78, 6.694, 6.562, 5.901, 6.668, 6.094, 5.526, 6.145, 7.413,
6.656, 8.64, 6.341, 5.441, 5.223, 8.157, 6.921, 7.693, 8.142, 8.211, 8.587,
7.184, 8.622, 8.398, 9.413, 9.323, 7.953, 5.593, 8.098, 7.654, 7.565, 6.256,
5.825, 7.313, 6.981, 6.511, 7.607, 8.949, 8.075, 9.148, 12.34, 4.679, 5.216,
5.783, 6.412, 6.349, 2.99, 3.227, 3.702, 4.948, 4.727, 8.471, 4.997, 5.624,
7.697, 11.58, 3.149, 4.015, 3.231, 3.806, 2.803, 4.779, 5.126, 5.19, 6.118,
5.167, 4.521, 6.065, 6.803, 5.867, 5.963, 5.983, 6.524, 6.41, 6.776, 6.435,
7.016, 4.066, 4.9, 4.872, 5.654, 5.719, 4.586, 7.102, 6.877, 5.722, 4.909,
4.532, 7.938, 7.823, 7.824, 6.839, 7.46, 7.16, 8.274, 8.225, 6.172, 6.735,
5.956, 5.422, 5.878, 5.581, 6.388, 7.606, 6.05, 5.789, 6.068, 6.279, 4.046,
7.339, 7.588, 7.662, 8.735, 8.412, 8.348, 8.365, 5.784, 5.564, 7.823, 6.926,
7.075, 5.503, 6.862, 7.419, 8.031, 8.062, 5.486, 4.724, 5.819, 6.948, 7.522,
5.961, 9.516, 9.755, 8.525, 8.24, 6.315, 6.479, 9.976, 9.494, 11.79, 3.432,
3.654, 4.721, 5.383, 3.776, 5.107, 4.938, 4.268, 5.405, 6.003, 2.221, 4.282,
4.589, 5.198, 6.697, 7.333, 8.904, 7.758, 9.524, 8.526, 7.523, 7.299, 7.085,
8.318, 6.852, 7.622, 8.402, 7.495, 7.647, 7.444, 6.431, 6.246, 6.782, 6.978,
7.565, 1.784, 2.358, 2.016, 1.766, 1.956, 5.784, 5.845, 6.352, 5.816, 5.759,
4.798, 4.218, 4.421, 5.451, 5.163, 7.282, 9.043, 7.192, 2.343, 3.572, 7.301,
8.002, 4.685, 6.632, 7.133, 6.572, 6.293, 5.604, 6.861, 5.391, 6.473, 7.589,
7.415, 7.737, 7.159, 5.389, 8.577, 4.672, 8.241, 10.4, 9.007, 9.401, 10.23,
11.64, 12.92, 5.751, 5.434, 3.199, 7.059, 3.778, 6.293, 6.405, 6.774, 6.188,
5.64, 5.779, 5.42, 5.823, 6.93, 7.383, 6.205, 5.685, 5.795, 6.22, 5.91, 7.546,
6.372, 7.934, 9.347, 9.418, 10.11, 10.04, 7.147, 20.59, 7.695, 5.091, 5.54,
2.702, 4.725, 3.933, 3.988, 5.813, 4.855, 5.631, 4.822, 7.251, 6.893, 7.207,
6.659, 7.394, 8.199, 8.426, 8.913, 8.276, 9.102, 6.177, 5.043, 3.822, 6.878,
6.02, 3.698, 3.599, 3.732, 4.231, 3.278, 5.928, 6.457, 6.031, 7.472, 7.293,
3.084, 4.429, 4.236, 3.67, 2.503, 2.453, 4.155, 5.43, 3.52, 4.495, 7.123, 7.68,
6.818, 6.966, 6.563, 7.745, 6.953, 8.628, 8.231, 6.828, 2.191, 3.85, 4.366,
3.984, 4.788, 4.447, 4.736, 4.329, 6.579, 4.406, 5.97, 6.326, 7.564, 8.222,
8.055, 7.565, 7.136, 5.377, 6.219, 6.646, 6.532, 8.403, 8.083, 7.22, 9.511,
3.396, 5.285, 6.356, 8.104, 5.739, 9.549, 7.706, 8.2, 8.947, 10.45, 4.981,
4.927, 6.032, 7.07, 7.651, 7.05, 6.557, 6.137, 6.96, 7.295, 7.838, 6.673,
9.052, 9.441, 8.707, 5.345, 5.711, 6.153, 6.169, 6.493, 7.488, 7.746, 7.812,
5.994, 8.129, 5.019, 7.132, 7.923, 7.35, 8.197, 3.076, 5.054, 3.948, 4.338,
4.6, 2.565, 3.902, 3.85, 6.226, 3.583, 3.704, 3.177, 4.065, 3.853, 4.602,
5.833, 7.863, 8.661, 8.308, 8.523, 3.722, 5.018, 4.456, 6.026, 6.245, 7.875,
4.679, 8.467, 8.445, 8.397, 4.971, 6.365, 6.893, 7.6, 6.336, 6.532, 9.058,
9.051, 8.245, 6.998, 7.364, 2.547, 5.639, 7.139, 5.202, 7.517, 8.35, 7.577,
9.506, 10.38, 5.455, 5.842, 8.374, 6.221, 6.517, 7.193, 8.509, 8.734, 9.154,
9.53, 5.757, 5.315, 6.944, 6.109, 6.141, 4.873, 4.859, 6.145, 7.293, 6.994,
3.756, 4.943, 4.524, 5.439, 5.777, 6.224, 6.361, 6.393, 5.453, 6.266, 6.846,
8.138, 6.54, 8.322, 8.382, 6.733, 7.003, 4.786, 5.749, 6.685, 5.265, 4.431,
5.329, 6.567, 5.311, -0.3581, 5.737, 4.634, 3.114, 9.494, 7.227, 7.356, 7.643,
8.448, 8.791, 7.538, 8.216, 8.618, 10.42, 9.5, 5.274, 5.518, 7.777, 6.539,
5.957, 3.88, 4.227, 4.809, 4.451, 5.603, 9.948, 11.02, 10.97, 11.8, 9.583,
7.287, 6.728, 6.723, 6.736, 6.875, 3.352, 6.444, 3.66, 5.925, 5.543, 7.212,
7.471, 9.594, 9.799, 8.405, 7.038, 6.181, 8.681, 7.441, 7.731, 7.253, 7.169,
7.799, 8.145, 7.357, 6.8, -1.434, 12.03, 4.352, 4.215, 5.161, 6.56, 5.629,
7.243, 6.405, 6.576, 8.231, 8.411, 8.316, 10.44, 5.19, 5.316, 7.209, 6.557,
7.775, 4.692, 5.133, 5.957, 5.916, 6.44, 4.31, 5.57, 5.31, 6.698, 6.399, 6.73,
6.79, 7.969, 8.1, 8.346, 7.495, 7.454, 6.961, 7.005, 7.449, 3.963, 4.365,
4.151, 4.454, 5.453, 4.917, 6.649, 5.02, 4.13, 6.823, 6.419, 7.004, 5.291,
6.885, 7.872, 6.468, 7.546, 6.475, 6.88, 8.402, 5.347, 5.219, 5.298, 6.187,
4.621, 7.288, 7.433, 10.94, 8.383, 10.51, 7.21, 4.771, 6.828, 6.2, 6.023,
5.239, 4.623, 5.173, 4.131, 4.89, -0.5749, 6.414, 9.2, 3.507, 5.7, 9.201,
9.488, 11.5, 11.85, 11.14, 6.726, 6.484, 6.199, 7.166, 7.771, 5.885, 5.407,
5.928, 5.561, 6.279, 3.668, 4.838, 4.079, 5.412, 6.115, 7.295, 8.262, 9.142,
7.847, 8.513, 5.196, 4.96, 4.994, 4.523, 3.994, 2.557, 7.596, 8.013, 6.469,
7.064, 4.138, 5.285, 3.73, 3.533, 7.846, 3.594, 3.746, 5.876, 4.295, 4.667,
5.321, 5.071, 5.237, 5.563, 6.172, 6.268, 5.806, 6.503, 6.03, 5.929, 4.025,
2.315, 2.399, 2.675, 3.661, 10.57, 10.09, 10.29, 11.22, 10.47, 6.951, 7.397,
8.602, 8.441, 9.115, 5.29, 5.419, 5.007, 6.208, 4.965, 7.273, 6.826, 6.381,
8.494, 7.822, 8.023, 8.336, 8.674, 8.912, 9.728, 7.524, 6.004, 8.512, 6.37,
7.233, 6.511, 3.62, 4.196, 6.557, 5.965, 5.25, 5.589, 6.149, 6.094, 6.138,
4.633, 5.338, 5.918, 5.875, 5.436, 7.402, 5.466, 7.074, 5.551, 5.746, 5.592,
6.479, 8.126, 7.395, 8.688, 8.511, 8.035, 7.554, 7.959, 8.684, 7.256, 4.476,
6.477, 4.928, 6.79, 4.327, 6.839, 8.447, 3.665, 5.941, 4.37, 4.465, 4.387,
5.888, 5.403, 5.979, 6.596, 6.894, 7.787, 3.905, 5.321, 5.901, 6.118, 7.073,
5.784, 9.318, 7.15, 5.813, 7.476, 7.154, 7.661, 9.01, 7.376, 9.038, 8.584,
3.315, 5.09, 5.179, 2.791, 7.328, 4.592, 4.818, 4.287, 8.401, 5.664, 9.621,
8.029, 10.75, 6.658, 6.896, 8.01, 7.241, 6.974, 8.567, 7.948, 6.137, 4.404,
4.56, 4.608, 6.989, 6.208, 6.931, 6.968, 7.336, 7.1, 4.473, 4.146, 3.482,
4.169, 5.15, 5.694, 3.706, 4.099, 5.835, 2.178, 4.102, 5.6, 5.551, 5.946, 5.98,
5.199, 5.359, 5.503, 5.448, 6.772, 7.333, 6.578, 6.601, 8.574, 7.818, 4.523,
4.223, 5.719, 5.893, 7.212, 9.77, 10.67, 10.95, 9.244, 11.23, 4.722, 2.944,
4.505, 5.602, 4.689, 7.194, 6.78, 7.903, 9.028, 8.618, 2.904, 3.021, 2.549,
3.285, 3.805, 5.905, 6.004, 7.585, 8.41, 8.312, 8.888, 8.297, 9.523, 9.189,
9.601, 5.126, 4.367, 6.361, 6.148, 6.389, 2.124, 3.615, 6.298, 6.646, 6.664,
3.751, 4.91, 5.147, 4.51, 4.467, 3.889, 3.419, 6.289, 5.553, 6.247, 6.249,
5.338, 5.784, 5.975, 5.879, 7.825, 6.28, 6.081, 6.336, 8.148, 2.079, 6.612,
4.72, 3.646, 4.182, 4.873, 5.145, 5.687, 6.505, 6.2, 4.845, 4.056, 4.523,
3.722, 3.662, 2.963, 3.875, 5.504, 5.535, 5.2, 9.468, 8.86, 9.459, 9.307,
9.122, 8.294, 8.467, 8.009, 8.505, 9.709, 5.358, 6.512, 3.459, 4.746, 6.021,
9.634, 9.5, 9.852, 10.82, 10.38, 4.887, 5.769, 6.389, 8.833, 5.394, 5.389,
8.16, 12.04, 9.56, 10.13, 5.221, 5.086, 4.903, 6.401, 5.914, 5.042, 5.77,
5.302, 6.589, 6.093, 3.322, 4.922, 6.393, 7.022, 4.507, 4.875, 6.941, 6.028,
8.349, 7.237, 6.384, 7.391, 6.534, 7.932, 7.35, 6.597, 10.43, 5.494, 7.65,
7.467, 4.753, 6.009, 5.762, 5.618, 5.965, 6.513, 6.588, 7.507, 7.013, 7.14,
7.066, 7.374, 8.601, 7.83, 9.618, 4.315, 5.714, 3.023, 5.537, 5.317, 6.933,
7.806, 7.598, 8.017, 9.077, 9.388, 6.48, 9.938, 6.041, 8.752, 6.525, 7.179,
7.209, 7.978, 8.229, 3.931, 9.302, 3.56, 11.41, 7.371, 7.041, 8.327, 7.109,
8.574, 8.066, 4.29, 4.729, 4.568, 4.133, 4.363, 8.415, 9.079, 8.162, 9.974,
9.423, 8.315, 8.9, 8.833, 9.191, 7.952, 6.464, 6.594, 5.638, 6.41, 5.059,
7.074, 6.069, 7.125, 8.451, 6.927, 6.775, 7.648, 8.872, 8.361, 8.911, 3.473,
3.388, 4.789, 5.447, 6.258, 7.215, 8.049, 8.764, 7.308, 6.579, 4.153, 4.57,
4.841, 5.537, 5.807, 9.329, 5.775, 3.885, 6.085, 6.556, 6.266, 7.252, 6.696,
8.145, 8.486, 6.546, 5.312, 5.699, 6.391, 6.323, 7.176, 6.09, 8.218, 6.666,
7.507, 4.414, 5.702, 6.541, 5.616, 5.528, 6.584, 6.523, 6.92, 5.64, 6.544,
3.917, 1.925, 4.058, 2.579, 5.357, 5.088, 5.426, 6.182, 6.499, 5.495, 5.065,
4.386, 6.119, 5.384, 5.369, 7.886, 5.818, 5.732, 7.204, 6.154, 7.256, 7.058,
8.391, 9.67, 10.59, 6.57, 6.168, 6.904, 7.697, 7.762, 6.982, 5.596, 5.372,
8.193, 4.991, 3.762, 2.575, 4.051, 5.2, 4.786, 6.516, 7.757, 8.336, 9.421,
9.267, 8.612, 9.025, 8.973, 9.269, 10.1, 5.904, 11.22, 5.087, 4.178, 2.675,
4.819, 6.604, 7.587, 6.392, 7.194, 8.092, 8.251, 8.586, 8.566, 7.406), .Dim=c(
500, 5 )) )