The function join clusters that have the minimum ratio of overlap as specified by user.
consolidate(emfit, minover)
| emfit | an object of class |
|---|---|
| minover | a single |
a list containing information about the updated
clusters obtained from a mixture model estimation:
mu a numeric vector representing the mean
for each component. If there is more than one component, the
kth element is the mean of the kth component of the
mixture model.
pro a vector whose kth component is the mixing
proportion for the kth component of the mixture model. If
missing, equal proportions are assumed.
z a numeric matrix whose [i,k]th entry
is the probability that observation i in the test data belongs
to the kth class.
groups a matrix of double, used as
integer, with the number of rows corresponding
to the current number of clusters while the number of columns is
corresponding to the initial number of clusters. The presence of
1 in position [i,k] indicates that the initial
ith cluster is now part of the new kth cluster.
ngroups a numeric, used as an integer, giving the final
number of clusters.
sigmasq a numeric vector giving the common
variance for each component in the mixture model "E".
Alexander Krasnitz, Guoli Sun
## Load Mclust object data(EMexample) ## Group clusters that have at least 0.4% of overlap ## The inital object has 5 clusters while the return object has only ## 4 clusters CNprep:::consolidate(EMexample, minover=0.004)#> $mu #> 1 2 3 5 #> -0.233274516 -0.081386468 0.004774756 0.245866821 #> #> $pro #> [1] 0.09326425 0.29018379 0.45852088 0.15803109 #> #> $z #> [,1] [,2] [,3] [,4] #> [1,] 1.518326e-120 2.430697e-24 1.000000e+00 3.216325e-67 #> [2,] 3.595740e-126 9.516728e-27 1.000000e+00 1.082023e-62 #> [3,] 1.575231e-121 9.218454e-25 1.000000e+00 1.991578e-66 #> [4,] 2.685068e-133 8.488695e-30 1.000000e+00 5.877146e-57 #> [5,] 8.480066e-123 2.640361e-25 1.000000e+00 2.090796e-65 #> [6,] 3.060337e-122 4.572648e-25 1.000000e+00 7.443863e-66 #> [7,] 3.187284e-118 2.395315e-23 1.000000e+00 4.353225e-69 #> [8,] 1.000000e+00 3.496168e-23 2.767048e-52 0.000000e+00 #> [9,] 1.000000e+00 6.425580e-24 2.616218e-53 0.000000e+00 #> [10,] 1.000000e+00 4.436161e-32 1.134045e-64 0.000000e+00 #> [11,] 1.000000e+00 3.620265e-32 8.545341e-65 0.000000e+00 #> [12,] 1.000000e+00 1.228463e-26 4.289744e-57 0.000000e+00 #> [13,] 1.000000e+00 3.571670e-30 5.108677e-62 0.000000e+00 #> [14,] 1.000000e+00 5.952846e-26 3.861121e-56 0.000000e+00 #> [15,] 1.000000e+00 1.886382e-30 2.100312e-62 0.000000e+00 #> [16,] 1.000000e+00 3.498971e-25 4.547266e-55 0.000000e+00 #> [17,] 1.000000e+00 2.736723e-28 2.148049e-59 0.000000e+00 #> [18,] 1.000000e+00 1.319747e-31 5.174836e-64 0.000000e+00 #> [19,] 1.000000e+00 4.232806e-32 1.062319e-64 0.000000e+00 #> [20,] 1.000000e+00 7.019656e-30 1.308861e-61 0.000000e+00 #> [21,] 6.074535e-120 4.399429e-24 1.000000e+00 1.053935e-67 #> [22,] 8.595568e-121 1.905453e-24 1.000000e+00 5.083848e-67 #> [23,] 8.595568e-121 1.905453e-24 1.000000e+00 5.083848e-67 #> [24,] 8.114711e-120 4.979782e-24 1.000000e+00 8.348634e-68 #> [25,] 1.005949e-123 1.060459e-25 1.000000e+00 1.162298e-64 #> [26,] 7.594533e-125 3.510316e-26 1.000000e+00 9.294927e-64 #> [27,] 8.541281e-123 2.648500e-25 1.000000e+00 2.078729e-65 #> [28,] 1.678264e-117 4.876116e-23 1.000000e+00 1.143593e-69 #> [29,] 1.504708e-122 3.374704e-25 1.000000e+00 1.317952e-65 #> [30,] 3.218946e-121 1.251608e-24 1.000000e+00 1.120583e-66 #> [31,] 2.091672e-122 3.885479e-25 1.000000e+00 1.011102e-65 #> [32,] 1.500923e-123 1.258511e-25 1.000000e+00 8.423168e-65 #> [33,] 2.254461e-124 5.591801e-26 1.000000e+00 3.872522e-64 #> [34,] 7.112440e-119 1.260727e-23 1.000000e+00 1.455432e-68 #> [35,] 3.836090e-122 5.036793e-25 1.000000e+00 6.206432e-66 #> [36,] 1.504708e-122 3.374704e-25 1.000000e+00 1.317952e-65 #> [37,] 4.948763e-124 7.828227e-26 1.000000e+00 2.056967e-64 #> [38,] 8.541281e-123 2.648500e-25 1.000000e+00 2.078729e-65 #> [39,] 1.678264e-117 4.876116e-23 1.000000e+00 1.143593e-69 #> [40,] 7.484915e-127 4.862093e-27 1.000000e+00 3.825762e-62 #> [41,] 2.813406e-125 2.295101e-26 1.000000e+00 2.066741e-63 #> [42,] 3.075372e-119 8.806631e-24 1.000000e+00 2.857581e-68 #> [43,] 5.259733e-120 4.136479e-24 1.000000e+00 1.183442e-67 #> [44,] 3.801921e-118 2.583051e-23 1.000000e+00 3.777336e-69 #> [45,] 2.091672e-122 3.885479e-25 1.000000e+00 1.011102e-65 #> [46,] 1.194917e-125 1.590978e-26 1.000000e+00 4.116687e-63 #> [47,] 1.000000e+00 3.748095e-28 3.328240e-59 0.000000e+00 #> [48,] 1.000000e+00 1.498617e-34 4.107990e-68 0.000000e+00 #> [49,] 1.000000e+00 5.817924e-38 7.317193e-73 0.000000e+00 #> [50,] 1.000000e+00 9.356656e-30 1.952848e-61 0.000000e+00 #> [51,] 1.000000e+00 5.709093e-32 1.611308e-64 0.000000e+00 #> [52,] 1.000000e+00 8.698358e-31 7.147891e-63 0.000000e+00 #> [53,] 1.000000e+00 1.590717e-33 1.101719e-66 0.000000e+00 #> [54,] 1.000000e+00 2.248514e-23 1.496590e-52 0.000000e+00 #> [55,] 1.000000e+00 4.076632e-28 3.741310e-59 0.000000e+00 #> [56,] 1.000000e+00 3.527140e-26 1.863046e-56 0.000000e+00 #> [57,] 1.000000e+00 5.155694e-36 3.767232e-70 0.000000e+00 #> [58,] 1.000000e+00 2.573032e-35 3.532812e-69 0.000000e+00 #> [59,] 1.000000e+00 5.806512e-35 1.097195e-68 0.000000e+00 #> [60,] 1.000000e+00 6.543596e-23 6.623004e-52 0.000000e+00 #> [61,] 1.467280e-114 8.847813e-22 1.000000e+00 4.910550e-72 #> [62,] 9.141573e-113 5.184668e-21 1.000000e+00 1.766387e-73 #> [63,] 1.366585e-115 3.204160e-22 1.000000e+00 3.316463e-71 #> [64,] 2.296091e-124 5.635754e-26 1.000000e+00 3.815922e-64 #> [65,] 1.591224e-117 4.766251e-23 1.000000e+00 1.193667e-69 #> [66,] 1.106904e-113 2.100717e-21 1.000000e+00 9.658889e-73 #> [67,] 1.391110e-112 6.205058e-21 1.000000e+00 1.259956e-73 #> [68,] 1.141207e-112 5.700947e-21 1.000000e+00 1.477603e-73 #> [69,] 9.520600e-115 7.352816e-22 1.000000e+00 6.954916e-72 #> [70,] 2.387992e-117 5.670457e-23 1.000000e+00 8.610213e-70 #> [71,] 0.000000e+00 8.283441e-172 6.467252e-74 1.000000e+00 #> [72,] 0.000000e+00 6.346946e-173 1.208873e-74 1.000000e+00 #> [73,] 0.000000e+00 1.291930e-168 7.857038e-72 1.000000e+00 #> [74,] 0.000000e+00 7.185108e-170 1.191435e-72 1.000000e+00 #> [75,] 0.000000e+00 9.342815e-176 1.711901e-76 1.000000e+00 #> [76,] 0.000000e+00 4.331409e-176 1.036403e-76 1.000000e+00 #> [77,] 0.000000e+00 9.231175e-170 1.403185e-72 1.000000e+00 #> [78,] 0.000000e+00 1.448733e-167 3.806957e-71 1.000000e+00 #> [79,] 0.000000e+00 2.817379e-167 5.876999e-71 1.000000e+00 #> [80,] 0.000000e+00 4.375713e-171 1.916981e-73 1.000000e+00 #> [81,] 0.000000e+00 7.331663e-173 1.328229e-74 1.000000e+00 #> [82,] 0.000000e+00 1.494710e-172 2.114619e-74 1.000000e+00 #> [83,] 0.000000e+00 8.988854e-167 1.253416e-70 1.000000e+00 #> [84,] 0.000000e+00 1.343972e-175 2.170553e-76 1.000000e+00 #> [85,] 0.000000e+00 7.907577e-170 1.268339e-72 1.000000e+00 #> [86,] 0.000000e+00 5.924443e-169 4.722941e-72 1.000000e+00 #> [87,] 0.000000e+00 2.817379e-167 5.876999e-71 1.000000e+00 #> [88,] 0.000000e+00 1.378949e-165 7.451823e-70 1.000000e+00 #> [89,] 0.000000e+00 6.309418e-168 2.212601e-71 1.000000e+00 #> [90,] 0.000000e+00 1.378949e-165 7.451823e-70 1.000000e+00 #> [91,] 0.000000e+00 2.792448e-175 3.498720e-76 1.000000e+00 #> [92,] 0.000000e+00 1.291930e-168 7.857038e-72 1.000000e+00 #> [93,] 0.000000e+00 7.317298e-172 5.964276e-74 1.000000e+00 #> [94,] 0.000000e+00 5.043903e-170 9.456940e-73 1.000000e+00 #> [95,] 0.000000e+00 3.132632e-170 6.929561e-73 1.000000e+00 #> [96,] 0.000000e+00 8.060812e-169 5.774544e-72 1.000000e+00 #> [97,] 0.000000e+00 1.080066e-169 1.554656e-72 1.000000e+00 #> [98,] 0.000000e+00 7.291707e-175 6.546995e-76 1.000000e+00 #> [99,] 0.000000e+00 8.988854e-167 1.253416e-70 1.000000e+00 #> [100,] 0.000000e+00 7.039123e-177 3.164937e-77 1.000000e+00 #> [101,] 0.000000e+00 3.543831e-175 4.087618e-76 1.000000e+00 #> [102,] 0.000000e+00 2.410948e-173 6.426008e-75 1.000000e+00 #> [103,] 0.000000e+00 3.184665e-171 1.557886e-73 1.000000e+00 #> [104,] 0.000000e+00 1.326825e-175 2.152433e-76 1.000000e+00 #> [105,] 0.000000e+00 6.496252e-171 2.481159e-73 1.000000e+00 #> [106,] 0.000000e+00 6.484529e-166 4.553511e-70 1.000000e+00 #> [107,] 0.000000e+00 6.407168e-164 9.133246e-69 1.000000e+00 #> [108,] 0.000000e+00 6.496252e-171 2.481159e-73 1.000000e+00 #> [109,] 0.000000e+00 7.440912e-167 1.107928e-70 1.000000e+00 #> [110,] 0.000000e+00 1.068310e-173 3.777185e-75 1.000000e+00 #> [111,] 0.000000e+00 2.308750e-166 2.320312e-70 1.000000e+00 #> [112,] 0.000000e+00 3.958594e-171 1.795619e-73 1.000000e+00 #> [113,] 0.000000e+00 2.820920e-171 1.439291e-73 1.000000e+00 #> [114,] 0.000000e+00 8.283441e-172 6.467252e-74 1.000000e+00 #> [115,] 0.000000e+00 1.911025e-166 2.050885e-70 1.000000e+00 #> [116,] 0.000000e+00 1.240386e-175 2.059822e-76 1.000000e+00 #> [117,] 0.000000e+00 8.283441e-172 6.467252e-74 1.000000e+00 #> [118,] 0.000000e+00 1.279341e-171 8.589345e-74 1.000000e+00 #> [119,] 0.000000e+00 7.331663e-173 1.328229e-74 1.000000e+00 #> [120,] 0.000000e+00 6.465973e-169 5.000478e-72 1.000000e+00 #> [121,] 0.000000e+00 1.800991e-168 9.759911e-72 1.000000e+00 #> [122,] 0.000000e+00 2.095404e-174 1.304156e-75 1.000000e+00 #> [123,] 0.000000e+00 9.691452e-166 5.919348e-70 1.000000e+00 #> [124,] 0.000000e+00 3.389203e-167 6.630520e-71 1.000000e+00 #> [125,] 0.000000e+00 5.217176e-171 2.150237e-73 1.000000e+00 #> [126,] 0.000000e+00 6.484529e-166 4.553511e-70 1.000000e+00 #> [127,] 0.000000e+00 6.465973e-169 5.000478e-72 1.000000e+00 #> [128,] 0.000000e+00 4.259299e-173 9.317323e-75 1.000000e+00 #> [129,] 0.000000e+00 7.331663e-173 1.328229e-74 1.000000e+00 #> [130,] 0.000000e+00 2.308750e-166 2.320312e-70 1.000000e+00 #> [131,] 0.000000e+00 3.954119e-169 3.627230e-72 1.000000e+00 #> [132,] 4.369643e-68 3.892733e-05 9.999611e-01 1.155221e-119 #> [133,] 2.530440e-64 4.471794e-04 9.995528e-01 1.995822e-124 #> [134,] 9.666337e-68 4.868799e-05 9.999513e-01 4.229181e-120 #> [135,] 3.117329e-63 9.071283e-04 9.990929e-01 8.306549e-126 #> [136,] 7.462601e-68 4.526434e-05 9.999547e-01 5.867828e-120 #> [137,] 4.103997e-66 1.400082e-04 9.998600e-01 3.679909e-122 #> [138,] 2.755606e-68 3.418463e-05 9.999658e-01 2.070521e-119 #> [139,] 5.268137e-68 4.103359e-05 9.999590e-01 9.117637e-120 #> [140,] 5.268137e-68 4.103359e-05 9.999590e-01 9.117637e-120 #> [141,] 1.263721e-65 1.922127e-04 9.998078e-01 8.863600e-123 #> [142,] 3.718600e-66 1.361713e-04 9.998638e-01 4.169112e-122 #> [143,] 3.376563e-66 1.325189e-04 9.998675e-01 4.710659e-122 #> [144,] 2.755606e-68 3.418463e-05 9.999658e-01 2.070521e-119 #> [145,] 4.994915e-65 2.831101e-04 9.997169e-01 1.556368e-123 #> [146,] 1.220112e-65 1.903202e-04 9.998097e-01 9.266468e-123 #> [147,] 2.150715e-65 2.232799e-04 9.997767e-01 4.521808e-123 #> [148,] 2.150715e-65 2.232799e-04 9.997767e-01 4.521808e-123 #> [149,] 1.292112e-66 1.010943e-04 9.998989e-01 1.588874e-121 #> [150,] 5.279335e-66 1.503043e-04 9.998497e-01 2.675503e-122 #> [151,] 3.620436e-65 2.585689e-04 9.997414e-01 2.338977e-123 #> [152,] 3.718600e-66 1.361713e-04 9.998638e-01 4.169112e-122 #> [153,] 5.268137e-68 4.103359e-05 9.999590e-01 9.117637e-120 #> [154,] 1.342037e-64 3.740134e-04 9.996260e-01 4.454313e-124 #> [155,] 4.367095e-69 2.034162e-05 9.999797e-01 2.131102e-118 #> [156,] 4.384099e-65 2.728947e-04 9.997271e-01 1.835760e-123 #> [157,] 4.370816e-64 5.216138e-04 9.994784e-01 9.991649e-125 #> [158,] 1.484489e-68 2.871633e-05 9.999713e-01 4.529769e-119 #> [159,] 2.277346e-65 2.269082e-04 9.997731e-01 4.205940e-123 #> [160,] 1.743493e-66 1.099998e-04 9.998900e-01 1.087428e-121 #> [161,] 2.918474e-66 1.271850e-04 9.998728e-01 5.665309e-122 #> [162,] 5.084797e-64 5.443266e-04 9.994557e-01 8.249938e-125 #> [163,] 1.530043e-69 1.513661e-05 9.999849e-01 8.036666e-118 #> [164,] 7.462601e-68 4.526434e-05 9.999547e-01 5.867828e-120 #> [165,] 5.220215e-68 4.092806e-05 9.999591e-01 9.223699e-120 #> [166,] 5.084797e-64 5.443266e-04 9.994557e-01 8.249938e-125 #> [167,] 1.457641e-65 2.001019e-04 9.997999e-01 7.398365e-123 #> [168,] 3.376563e-66 1.325189e-04 9.998675e-01 4.710659e-122 #> [169,] 3.868201e-68 3.761293e-05 9.999624e-01 1.347918e-119 #> [170,] 2.044620e-65 2.201199e-04 9.997799e-01 4.820817e-123 #> [171,] 4.384099e-65 2.728947e-04 9.997271e-01 1.835760e-123 #> [172,] 1.274315e-64 3.685969e-04 9.996314e-01 4.756051e-124 #> [173,] 9.248738e-67 9.200398e-05 9.999080e-01 2.425974e-121 #> [174,] 9.666337e-68 4.868799e-05 9.999513e-01 4.229181e-120 #> [175,] 2.093603e-65 2.215931e-04 9.997784e-01 4.678500e-123 #> [176,] 2.734770e-65 2.389173e-04 9.997611e-01 3.336156e-123 #> [177,] 2.179954e-67 6.122752e-05 9.999388e-01 1.510957e-120 #> [178,] 1.827860e-65 2.132781e-04 9.997867e-01 5.555511e-123 #> [179,] 3.393356e-68 3.625006e-05 9.999637e-01 1.590932e-119 #> [180,] 9.644386e-68 4.865680e-05 9.999513e-01 4.241368e-120 #> [181,] 1.889495e-69 1.606399e-05 9.999839e-01 6.153086e-118 #> [182,] 1.831686e-66 1.115400e-04 9.998885e-01 1.021589e-121 #> [183,] 1.056499e-63 6.688367e-04 9.993312e-01 3.268727e-125 #> [184,] 2.406718e-67 6.295895e-05 9.999370e-01 1.333084e-120 #> [185,] 6.279060e-65 3.019613e-04 9.996980e-01 1.165025e-123 #> [186,] 2.831578e-68 3.444763e-05 9.999656e-01 2.000464e-119 #> [187,] 5.183941e-66 1.495341e-04 9.998505e-01 2.737970e-122 #> [188,] 2.277346e-65 2.269082e-04 9.997731e-01 4.205940e-123 #> [189,] 1.827860e-65 2.132781e-04 9.997867e-01 5.555511e-123 #> [190,] 5.545518e-66 1.524021e-04 9.998476e-01 2.514007e-122 #> [191,] 3.620436e-65 2.585689e-04 9.997414e-01 2.338977e-123 #> [192,] 1.467051e-66 1.047769e-04 9.998952e-01 1.352988e-121 #> [193,] 1.268130e-84 8.300913e-10 1.000000e+00 9.077704e-99 #> [194,] 1.268130e-84 8.300913e-10 1.000000e+00 9.077704e-99 #> [195,] 1.452199e-90 1.099598e-11 1.000000e+00 6.833210e-92 #> [196,] 1.347569e-84 8.447443e-10 1.000000e+00 8.415800e-99 #> [197,] 1.347569e-84 8.447443e-10 1.000000e+00 8.415800e-99 #> [198,] 4.621315e-84 1.203349e-09 1.000000e+00 1.805887e-99 #> [199,] 3.800254e-83 2.195324e-09 1.000000e+00 1.284808e-100 #> [200,] 7.903682e-84 1.403017e-09 1.000000e+00 9.222675e-100 #> [201,] 5.593962e-82 4.708793e-09 1.000000e+00 4.344722e-102 #> [202,] 1.783820e-88 5.722856e-11 1.000000e+00 3.916493e-94 #> [203,] 2.467169e-86 2.626240e-10 1.000000e+00 1.168800e-96 #> [204,] 3.631416e-88 7.191354e-11 1.000000e+00 1.740423e-94 #> [205,] 5.235316e-86 3.281901e-10 1.000000e+00 4.667475e-97 #> [206,] 9.938561e-86 3.962778e-10 1.000000e+00 2.126033e-97 #> [207,] 5.142054e-88 8.031672e-11 1.000000e+00 1.165705e-94 #> [208,] 2.399234e-80 1.362563e-08 1.000000e+00 3.772348e-104 #> [209,] 1.347569e-84 8.447443e-10 1.000000e+00 8.415800e-99 #> [210,] 1.562516e-85 4.523808e-10 1.000000e+00 1.217860e-97 #> [211,] 1.347569e-84 8.447443e-10 1.000000e+00 8.415800e-99 #> [212,] 5.593962e-82 4.708793e-09 1.000000e+00 4.344722e-102 #> [213,] 4.756398e-84 1.213318e-09 1.000000e+00 1.741941e-99 #> [214,] 4.745796e-79 3.163154e-08 1.000000e+00 8.661807e-106 #> [215,] 1.035232e-82 2.918667e-09 1.000000e+00 3.641601e-101 #> [216,] 1.000000e+00 6.606254e-60 2.036694e-103 0.000000e+00 #> [217,] 1.000000e+00 4.879277e-61 5.411455e-105 0.000000e+00 #> [218,] 1.000000e+00 6.606254e-60 2.036694e-103 0.000000e+00 #> [219,] 1.000000e+00 1.749200e-57 4.813826e-100 0.000000e+00 #> [220,] 1.000000e+00 3.670668e-63 5.976557e-108 0.000000e+00 #> [221,] 1.000000e+00 2.946826e-62 1.086443e-106 0.000000e+00 #> [222,] 4.186924e-71 5.490670e-06 9.999945e-01 7.637855e-116 #> [223,] 1.000000e+00 7.589947e-45 1.898869e-82 0.000000e+00 #> [224,] 1.107239e-89 2.258404e-11 1.000000e+00 8.314726e-93 #> [225,] 1.739927e-86 2.366580e-10 1.000000e+00 1.786038e-96 #> [226,] 2.301382e-91 5.561213e-12 1.000000e+00 4.211190e-91 #> [227,] 3.712437e-89 3.410759e-11 1.000000e+00 2.252283e-93 #> [228,] 4.814686e-94 4.844545e-13 1.000000e+00 1.131933e-88 #> [229,] 5.591864e-87 1.681616e-10 1.000000e+00 7.013222e-96 #> [230,] 2.100543e-94 3.444708e-13 1.000000e+00 2.305584e-88 #> [231,] 6.899423e-82 4.996912e-09 1.000000e+00 3.334712e-102 #> [232,] 2.712247e-82 3.835660e-09 1.000000e+00 1.082334e-101 #> [233,] 7.856415e-81 9.940959e-09 1.000000e+00 1.546324e-103 #> [234,] 1.000000e+00 4.952239e-63 9.068564e-108 0.000000e+00 #> [235,] 1.000000e+00 1.485520e-49 5.280679e-89 0.000000e+00 #> [236,] 3.761667e-109 1.825214e-19 1.000000e+00 2.180392e-76 #> [237,] 1.756439e-105 6.780483e-18 1.000000e+00 2.428752e-79 #> [238,] 4.192788e-108 5.121130e-19 1.000000e+00 3.132188e-77 #> [239,] 7.161270e-106 4.619148e-18 1.000000e+00 5.000724e-79 #> [240,] 5.438087e-110 7.978274e-20 1.000000e+00 1.033862e-75 #> [241,] 2.471519e-105 7.847269e-18 1.000000e+00 1.844915e-79 #> [242,] 1.461947e-110 4.547634e-20 1.000000e+00 2.975576e-75 #> [243,] 1.013985e-108 2.789883e-19 1.000000e+00 9.816715e-77 #> [244,] 2.431365e-107 1.086388e-18 1.000000e+00 7.612126e-78 #> [245,] 4.262542e-105 9.907902e-18 1.000000e+00 1.189656e-79 #> [246,] 4.262542e-105 9.907902e-18 1.000000e+00 1.189656e-79 #> [247,] 2.471519e-105 7.847269e-18 1.000000e+00 1.844915e-79 #> [248,] 9.550994e-107 1.950840e-18 1.000000e+00 2.530882e-78 #> [249,] 9.550994e-107 1.950840e-18 1.000000e+00 2.530882e-78 #> [250,] 7.161270e-106 4.619148e-18 1.000000e+00 5.000724e-79 #> [251,] 1.112379e-104 1.493450e-17 1.000000e+00 5.495845e-80 #> [252,] 5.436183e-111 2.978122e-20 1.000000e+00 6.596435e-75 #> [253,] 7.161270e-106 4.619148e-18 1.000000e+00 5.000724e-79 #> [254,] 5.371677e-106 4.084453e-18 1.000000e+00 6.303088e-79 #> [255,] 2.008257e-104 1.922830e-17 1.000000e+00 3.415534e-80 #> [256,] 1.461947e-110 4.547634e-20 1.000000e+00 2.975576e-75 #> [257,] 1.013985e-108 2.789883e-19 1.000000e+00 9.816715e-77 #> [258,] 1.166837e-106 2.125342e-18 1.000000e+00 2.154168e-78 #> [259,] 1.756439e-105 6.780483e-18 1.000000e+00 2.428752e-79 #> [260,] 9.740136e-109 2.742283e-19 1.000000e+00 1.013964e-76 #> [261,] 3.472859e-107 1.265420e-18 1.000000e+00 5.713331e-78 #> [262,] 3.780903e-110 6.829120e-20 1.000000e+00 1.385128e-75 #> [263,] 1.756439e-105 6.780483e-18 1.000000e+00 2.428752e-79 #> [264,] 5.593369e-105 1.112919e-17 1.000000e+00 9.559174e-80 #> [265,] 1.583396e-107 9.042509e-19 1.000000e+00 1.075007e-77 #> [266,] 5.436183e-111 2.978122e-20 1.000000e+00 6.596435e-75 #> [267,] 5.341685e-107 1.521395e-18 1.000000e+00 4.040102e-78 #> [268,] 1.461947e-110 4.547634e-20 1.000000e+00 2.975576e-75 #> [269,] 6.016379e-110 8.330844e-20 1.000000e+00 9.531118e-76 #> [270,] 3.761667e-109 1.825214e-19 1.000000e+00 2.180392e-76 #> [271,] 6.813198e-107 1.688323e-18 1.000000e+00 3.321560e-78 #> [272,] 2.046512e-109 1.406673e-19 1.000000e+00 3.558622e-76 #> [273,] 7.495186e-104 3.377446e-17 1.000000e+00 1.182723e-80 #> [274,] 4.192788e-108 5.121130e-19 1.000000e+00 3.132188e-77 #> [275,] 3.272035e-37 9.999989e-01 1.060857e-06 2.835052e-172 #> [276,] 1.078011e-36 9.999993e-01 6.644819e-07 2.172274e-173 #> [277,] 8.398456e-38 9.999982e-01 1.808842e-06 5.309917e-171 #> [278,] 1.048662e-37 9.999983e-01 1.657913e-06 3.290865e-171 #> [279,] 7.371203e-38 9.999981e-01 1.903852e-06 7.033441e-171 #> [280,] 2.233282e-36 9.999995e-01 4.993062e-07 4.522526e-174 #> [281,] 2.487209e-36 9.999995e-01 4.786477e-07 3.586035e-174 #> [282,] 3.628252e-38 9.999975e-01 2.514328e-06 3.239151e-170 #> [283,] 6.156720e-39 9.999950e-01 5.042966e-06 1.479774e-168 #> [284,] 3.762826e-37 9.999990e-01 1.004249e-06 2.097910e-172 #> [285,] 1.255505e-38 9.999962e-01 3.812929e-06 3.187313e-169 #> [286,] 2.771410e-36 9.999995e-01 4.587530e-07 2.840369e-174 #> [287,] 3.180413e-38 9.999974e-01 2.647715e-06 4.302318e-170 #> [288,] 8.398456e-38 9.999982e-01 1.808842e-06 5.309917e-171 #> [289,] 8.174763e-37 9.999993e-01 7.406725e-07 3.942590e-173 #> [290,] 8.041095e-38 9.999982e-01 1.839969e-06 5.831432e-171 #> [291,] 1.103165e-38 9.999960e-01 4.011451e-06 4.211760e-169 #> [292,] 4.911658e-38 9.999978e-01 2.232609e-06 1.686708e-170 #> [293,] 2.195281e-36 9.999995e-01 5.026798e-07 4.692886e-174 #> [294,] 3.536741e-38 9.999975e-01 2.539657e-06 3.422430e-170 #> [295,] 1.078011e-36 9.999993e-01 6.644819e-07 2.172274e-173 #> [296,] 1.110729e-37 9.999984e-01 1.620926e-06 2.907402e-171 #> [297,] 5.236649e-38 9.999978e-01 2.177182e-06 1.469225e-170 #> [298,] 6.545433e-39 9.999951e-01 4.923266e-06 1.296904e-168 #> [299,] 4.925850e-39 9.999945e-01 5.504204e-06 2.392787e-168 #> [300,] 8.505936e-38 9.999982e-01 1.799840e-06 5.166409e-171 #> [301,] 3.412952e-38 9.999974e-01 2.575409e-06 3.695491e-170 #> [302,] 3.032785e-37 9.999989e-01 1.092939e-06 3.338955e-172 #> [303,] 7.953619e-36 9.999997e-01 3.033363e-07 2.930053e-175 #> [304,] 5.475081e-38 9.999979e-01 2.139476e-06 1.334828e-170 #> [305,] 5.039078e-37 9.999991e-01 8.955167e-07 1.118170e-172 #> [306,] 2.481567e-36 9.999995e-01 4.790744e-07 3.603623e-174 #> [307,] 3.679694e-38 9.999975e-01 2.500477e-06 3.142373e-170 #> [308,] 3.791342e-38 9.999975e-01 2.471322e-06 2.946380e-170 #> [309,] 4.517285e-39 9.999943e-01 5.694416e-06 2.883524e-168 #> [310,] 4.911658e-38 9.999978e-01 2.232609e-06 1.686708e-170 #> [311,] 2.224873e-37 9.999988e-01 1.234198e-06 6.508458e-172 #> [312,] 1.393621e-36 9.999994e-01 6.007930e-07 1.249205e-173 #> [313,] 9.103727e-38 9.999982e-01 1.752507e-06 4.463087e-171 #> [314,] 3.855384e-38 9.999975e-01 2.455133e-06 2.841941e-170 #> [315,] 6.792155e-38 9.999980e-01 1.965959e-06 8.389212e-171 #> [316,] 1.289139e-37 9.999985e-01 1.528902e-06 2.109222e-171 #> [317,] 1.794113e-37 9.999987e-01 1.342933e-06 1.034747e-171 #> [318,] 1.794113e-37 9.999987e-01 1.342933e-06 1.034747e-171 #> [319,] 2.295246e-36 9.999995e-01 4.939730e-07 4.263558e-174 #> [320,] 2.012050e-39 9.999922e-01 7.821055e-06 1.646980e-167 #> [321,] 3.272035e-37 9.999989e-01 1.060857e-06 2.835052e-172 #> [322,] 8.026829e-39 9.999955e-01 4.544508e-06 8.356056e-169 #> [323,] 1.363173e-38 9.999963e-01 3.691800e-06 2.669541e-169 #> [324,] 1.771327e-39 9.999918e-01 8.222033e-06 2.167314e-167 #> [325,] 3.412952e-38 9.999974e-01 2.575409e-06 3.695491e-170 #> [326,] 3.615296e-37 9.999990e-01 1.020133e-06 2.286716e-172 #> [327,] 7.430984e-40 9.999884e-01 1.156117e-05 1.408423e-166 #> [328,] 9.103727e-38 9.999982e-01 1.752507e-06 4.463087e-171 #> [329,] 3.032785e-37 9.999989e-01 1.092939e-06 3.338955e-172 #> [330,] 1.030332e-36 9.999993e-01 6.763815e-07 2.394655e-173 #> [331,] 6.822481e-37 9.999992e-01 7.951344e-07 5.820842e-173 #> [332,] 4.911658e-38 9.999978e-01 2.232609e-06 1.686708e-170 #> [333,] 1.648444e-38 9.999966e-01 3.426560e-06 1.772701e-169 #> [334,] 2.651407e-37 9.999988e-01 1.152118e-06 4.460284e-172 #> [335,] 4.911658e-38 9.999978e-01 2.232609e-06 1.686708e-170 #> [336,] 1.289139e-37 9.999985e-01 1.528902e-06 2.109222e-171 #> [337,] 7.793898e-37 9.999992e-01 7.546689e-07 4.369434e-173 #> [338,] 7.446140e-37 9.999992e-01 7.683068e-07 4.820990e-173 #> [339,] 2.915634e-38 9.999973e-01 2.739577e-06 5.188459e-170 #> [340,] 4.103907e-36 9.999996e-01 3.932604e-07 1.219068e-174 #> [341,] 2.651407e-37 9.999988e-01 1.152118e-06 4.460284e-172 #> [342,] 1.653347e-39 9.999916e-01 8.447432e-06 2.514304e-167 #> [343,] 4.872912e-38 9.999978e-01 2.239558e-06 1.715736e-170 #> [344,] 3.412952e-38 9.999974e-01 2.575409e-06 3.695491e-170 #> [345,] 6.816158e-38 9.999980e-01 1.963240e-06 8.325691e-171 #> [346,] 3.791342e-38 9.999975e-01 2.471322e-06 2.946380e-170 #> [347,] 6.606534e-37 9.999992e-01 8.052329e-07 6.238530e-173 #> [348,] 6.816158e-38 9.999980e-01 1.963240e-06 8.325691e-171 #> [349,] 2.669222e-38 9.999972e-01 2.836158e-06 6.275704e-170 #> [350,] 6.816158e-38 9.999980e-01 1.963240e-06 8.325691e-171 #> [351,] 1.916511e-38 9.999968e-01 3.229851e-06 1.281287e-169 #> [352,] 3.454557e-38 9.999974e-01 2.563194e-06 3.600266e-170 #> [353,] 2.885203e-36 9.999995e-01 4.515667e-07 2.604488e-174 #> [354,] 1.289139e-37 9.999985e-01 1.528902e-06 2.109222e-171 #> [355,] 3.762826e-37 9.999990e-01 1.004249e-06 2.097910e-172 #> [356,] 2.960285e-37 9.999989e-01 1.103365e-06 3.517637e-172 #> [357,] 3.062136e-38 9.999973e-01 2.687382e-06 4.668363e-170 #> [358,] 3.749005e-39 9.999939e-01 6.126553e-06 4.308835e-168 #> [359,] 3.791342e-38 9.999975e-01 2.471322e-06 2.946380e-170 #> [360,] 9.242921e-38 9.999983e-01 1.742104e-06 4.319532e-171 #> [361,] 3.762826e-37 9.999990e-01 1.004249e-06 2.097910e-172 #> [362,] 9.242921e-38 9.999983e-01 1.742104e-06 4.319532e-171 #> [363,] 3.615296e-37 9.999990e-01 1.020133e-06 2.286716e-172 #> [364,] 3.791342e-38 9.999975e-01 2.471322e-06 2.946380e-170 #> [365,] 9.663583e-39 9.999958e-01 4.225362e-06 5.602169e-169 #> [366,] 5.299114e-39 9.999947e-01 5.348693e-06 2.044357e-168 #> [367,] 1.839916e-37 9.999987e-01 1.329715e-06 9.800438e-172 #> [368,] 1.143754e-38 9.999960e-01 3.954980e-06 3.896313e-169 #> [369,] 3.069486e-39 9.999934e-01 6.626647e-06 6.629552e-168 #> [370,] 3.046115e-38 9.999973e-01 2.692918e-06 4.721423e-170 #> [371,] 2.339504e-38 9.999970e-01 2.986745e-06 8.337478e-170 #> [372,] 8.505936e-38 9.999982e-01 1.799840e-06 5.166409e-171 #> [373,] 1.078011e-36 9.999993e-01 6.644819e-07 2.172274e-173 #> [374,] 2.528322e-39 9.999928e-01 7.150635e-06 1.006862e-167 #> [375,] 6.156720e-39 9.999950e-01 5.042966e-06 1.479774e-168 #> [376,] 2.651407e-37 9.999988e-01 1.152118e-06 4.460284e-172 #> [377,] 3.412952e-38 9.999974e-01 2.575409e-06 3.695491e-170 #> [378,] 7.818156e-38 9.999981e-01 1.860380e-06 6.195613e-171 #> [379,] 1.339101e-38 9.999963e-01 3.717700e-06 2.774010e-169 #> [380,] 1.443952e-39 9.999911e-01 8.908415e-06 3.366122e-167 #> [381,] 3.032785e-37 9.999989e-01 1.092939e-06 3.338955e-172 #> [382,] 1.289139e-37 9.999985e-01 1.528902e-06 2.109222e-171 #> [383,] 1.916511e-38 9.999968e-01 3.229851e-06 1.281287e-169 #> [384,] 3.311022e-38 9.999974e-01 2.606232e-06 3.944971e-170 #> [385,] 8.041095e-38 9.999982e-01 1.839969e-06 5.831432e-171 #> [386,] 5.279942e-38 9.999978e-01 2.170160e-06 1.443392e-170 #> #> $groups #> [,1] [,2] [,3] [,4] [,5] #> [1,] 1 0 0 0 0 #> [2,] 0 1 0 0 0 #> [3,] 0 0 1 1 0 #> [4,] 0 0 0 0 1 #> #> $ngroups #> [1] 4 #> #> $sigmasq #> [1] 0.0001367895 0.0001367895 0.0008660118 0.0001367895 #>