fitted.mix: Compute Mixture Model Fitted Values

In mixdist: Finite Mixture Distribution Models

Description

fitted.mix is a function which computes fitted values from objects returned by the modeling function mix. It is called via the generic function fitted.

Usage

## S3 method for class 'mix'
fitted(object, digits = NULL, ...)

Arguments

object	an object of class "mix", usually, the results returned by the model fitting function mix.
digits	a specified number of decimal places to be reserved.
...	other arguments.

Value

List with the following components:

mixed	the estimated mixed data, that is, the fitted numbers of observations falling into each interval.
joint	the estimated joint data, that is, the fitted numbers of observations from each component falling into every interval.
conditional	the estimated conditional data to be returned if usecondit of object is TRUE, which are the fitted numbers of observations from given intervals belonging to each component.
conditprob	the estimated conditional probabilities of observations from given interval belonging to each component.

See Also

mix for fitting mixture distributions.

Examples

data(pike65)
data(pikepar)
fit1 <- mix(pike65, pikepar, "lnorm", mixconstr(consigma = "CCV"), emsteps = 3)
fitted(fit1)
data(pike65sg)
fit2 <- mix(pike65sg, pikepar, "gamma", mixconstr(consigma = "CCV"), usecondit = TRUE)
fitted(fit2, digits = 2)

Example output

$mixed
 [1]  3.721598 11.913448 17.346125 13.796960 13.094478 26.715039 49.750928
 [8] 65.142476 63.745096 52.341610 40.850745 32.441799 26.035717 20.804333
[15] 16.741087 13.718550 11.348665  9.347472  7.653015  6.273290  5.157858
[22]  4.207481  3.343608  2.545647  4.962974

$joint
              [,1]         [,2]         [,3]         [,4]         [,5]
 [1,] 3.721559e+00 3.873151e-05 7.593157e-11 7.287403e-18 2.199903e-26
 [2,] 1.190974e+01 3.711410e-03 4.425021e-08 2.441086e-14 4.227216e-22
 [3,] 1.723380e+01 1.123216e-01 6.906035e-06 1.869170e-11 1.587396e-18
 [4,] 1.243291e+01 1.363671e+00 3.787348e-04 4.411015e-09 1.611421e-15
 [5,] 5.157334e+00 7.928223e+00 8.920763e-03 4.006669e-07 5.636155e-13
 [6,] 1.374700e+00 2.523531e+01 1.050132e-01 1.658353e-05 8.173938e-11
 [7,] 2.569587e-01 4.879910e+01 6.945143e-01 3.563888e-04 5.676154e-09
 [8,] 3.609028e-02 6.227199e+01 2.829985e+00 4.407369e-03 2.114481e-07
 [9,] 4.023885e-03 5.605589e+01 7.651122e+00 3.405591e-02 4.628361e-06
[10,] 3.721534e-04 3.758670e+01 1.457868e+01 1.757971e-01 6.409831e-05
[11,] 2.958013e-05 1.963291e+01 2.057667e+01 6.405421e-01 5.968675e-04
[12,] 2.079468e-06 8.290812e+00 2.242241e+01 1.724641e+00 3.930864e-03
[13,] 1.323598e-07 2.919589e+00 1.953105e+01 3.565972e+00 1.910396e-02
[14,] 7.776192e-09 8.799186e-01 1.400368e+01 5.849708e+00 7.102464e-02
[15,] 4.284293e-10 2.319943e-01 8.473245e+00 7.827553e+00 2.082948e-01
[16,] 2.243164e-11 5.451235e-02 4.419613e+00 8.749675e+00 4.947506e-01
[17,] 1.128527e-12 1.159759e-02 2.023810e+00 8.339548e+00 9.737098e-01
[18,] 5.787317e-14 2.264523e-03 8.264864e-01 6.899133e+00 1.619588e+00
[19,] 0.000000e+00 4.105571e-04 3.051357e-01 5.031070e+00 2.316398e+00
[20,] 0.000000e+00 6.980983e-05 1.030571e-01 3.277951e+00 2.892211e+00
[21,] 0.000000e+00 1.122997e-05 3.217138e-02 1.930903e+00 3.194773e+00
[22,] 0.000000e+00 1.722017e-06 9.366758e-03 1.039097e+00 3.159016e+00
[23,] 0.000000e+00 2.533668e-07 2.563771e-03 5.155619e-01 2.825482e+00
[24,] 0.000000e+00 3.597595e-08 6.643120e-04 2.377769e-01 2.307206e+00
[25,] 0.000000e+00 5.719545e-09 2.140011e-04 1.696008e-01 4.793159e+00

$conditprob
              [,1]         [,2]         [,3]         [,4]         [,5]
 [1,] 9.999896e-01 1.040723e-05 2.040295e-11 1.958138e-18 5.911177e-27
 [2,] 9.996885e-01 3.115311e-04 3.714308e-09 2.049017e-15 3.548272e-23
 [3,] 9.935243e-01 6.475312e-03 3.981313e-07 1.077572e-12 9.151298e-20
 [4,] 9.011340e-01 9.883850e-02 2.745060e-05 3.197092e-10 1.167954e-16
 [5,] 3.938556e-01 6.054631e-01 6.812615e-04 3.059816e-08 4.304223e-14
 [6,] 5.145791e-02 9.446106e-01 3.930866e-03 6.207562e-07 3.059677e-12
 [7,] 5.164903e-03 9.808681e-01 1.395983e-02 7.163461e-06 1.140914e-10
 [8,] 5.540207e-04 9.559353e-01 4.344301e-02 6.765738e-05 3.245933e-09
 [9,] 6.312463e-05 8.793757e-01 1.200268e-01 5.342514e-04 7.260733e-08
[10,] 7.110086e-06 7.181036e-01 2.785294e-01 3.358649e-03 1.224615e-06
[11,] 7.241026e-07 4.806009e-01 5.037037e-01 1.568006e-02 1.461093e-05
[12,] 6.409840e-08 2.555596e-01 6.911581e-01 5.316107e-02 1.211666e-04
[13,] 5.083777e-09 1.121378e-01 7.501638e-01 1.369646e-01 7.337599e-04
[14,] 3.737775e-10 4.229497e-02 6.731137e-01 2.811774e-01 3.413935e-03
[15,] 2.559148e-11 1.385778e-02 5.061347e-01 4.675654e-01 1.244213e-02
[16,] 1.635132e-12 3.973623e-03 3.221632e-01 6.377988e-01 3.606435e-02
[17,] 9.944137e-14 1.021935e-03 1.783302e-01 7.348483e-01 8.579950e-02
[18,] 6.191317e-15 2.422604e-04 8.841817e-02 7.380747e-01 1.732648e-01
[19,] 0.000000e+00 5.364645e-05 3.987131e-02 6.573972e-01 3.026779e-01
[20,] 0.000000e+00 1.112811e-05 1.642792e-02 5.225251e-01 4.610359e-01
[21,] 0.000000e+00 2.177255e-06 6.237353e-03 3.743613e-01 6.193991e-01
[22,] 0.000000e+00 4.092751e-07 2.226215e-03 2.469641e-01 7.508093e-01
[23,] 0.000000e+00 7.577648e-08 7.667679e-04 1.541933e-01 8.450399e-01
[24,] 0.000000e+00 1.413234e-08 2.609600e-04 9.340526e-02 9.063338e-01
[25,] 0.000000e+00 1.152443e-09 4.311953e-05 3.417322e-02 9.657837e-01

$mixed
 [1]  3.00 11.28 18.59 15.35 13.46 27.03 50.78 65.15 60.36 47.14 38.48 34.92
[13] 30.88 24.36 17.49 12.54  9.77  8.28  7.23  6.26  5.31  4.40  3.52  2.67
[25]  4.76

$joint
       [,1]  [,2]  [,3] [,4] [,5]
 [1,]  3.00  0.00  0.00 0.00 0.00
 [2,] 11.27  0.01  0.00 0.00 0.00
 [3,] 18.45  0.13  0.00 0.00 0.00
 [4,] 13.90  1.44  0.00 0.00 0.00
 [5,]  5.36  8.10  0.01 0.00 0.00
 [6,]  1.15 25.82  0.07 0.00 0.00
 [7,]  0.15 50.17  0.46 0.00 0.00
 [8,]  0.01 63.01  2.13 0.00 0.00
 [9,]  0.00 53.69  6.66 0.01 0.00
[10,]  0.00 32.34 14.77 0.03 0.00
[11,]  0.00 14.25 24.09 0.14 0.00
[12,]  0.00  4.74 29.72 0.47 0.00
[13,]  0.00  1.22 28.46 1.20 0.01
[14,]  0.00  0.25 21.63 2.44 0.04
[15,]  0.00  0.04 13.31 4.01 0.14
[16,]  0.00  0.01  6.74 5.45 0.34
[17,]  0.00  0.00  2.85 6.20 0.72
[18,]  0.00  0.00  1.02 5.98 1.28
[19,]  0.00  0.00  0.32 4.95 1.96
[20,]  0.00  0.00  0.08 3.56 2.61
[21,]  0.00  0.00  0.02 2.24 3.05
[22,]  0.00  0.00  0.00 1.25 3.15
[23,]  0.00  0.00  0.00 0.62 2.90
[24,]  0.00  0.00  0.00 0.27 2.39
[25,]  0.00  0.00  0.00 0.17 4.59

$conditional
      [,1] [,2]  [,3] [,4] [,5]
 [1,] 0.00 0.00  0.00 0.00 0.00
 [2,] 0.00 0.00  0.00 0.00 0.00
 [3,] 0.00 0.00  0.00 0.00 0.00
 [4,] 9.96 1.04  0.00 0.00 0.00
 [5,] 5.57 8.42  0.01 0.00 0.00
 [6,] 0.00 0.00  0.00 0.00 0.00
 [7,] 0.00 0.00  0.00 0.00 0.00
 [8,] 0.00 0.00  0.00 0.00 0.00
 [9,] 0.00 0.00  0.00 0.00 0.00
[10,] 0.00 0.00  0.00 0.00 0.00
[11,] 0.00 0.00  0.00 0.00 0.00
[12,] 0.00 4.88 30.63 0.48 0.00
[13,] 0.00 0.83 19.35 0.81 0.01
[14,] 0.00 0.00  0.00 0.00 0.00
[15,] 0.00 0.02  7.61 2.29 0.08
[16,] 0.00 0.01  6.45 5.22 0.33
[17,] 0.00 0.00  3.51 7.62 0.88
[18,] 0.00 0.00  1.36 7.94 1.70
[19,] 0.00 0.00  0.35 5.48 2.17
[20,] 0.00 0.00  0.04 1.71 1.25
[21,] 0.00 0.00  0.02 2.53 3.45
[22,] 0.00 0.00  0.01 1.70 4.30
[23,] 0.00 0.00  0.00 0.00 0.00
[24,] 0.00 0.00  0.00 0.00 0.00
[25,] 0.00 0.00  0.00 0.00 0.00

$conditprob
      [,1] [,2] [,3] [,4] [,5]
 [1,] 1.00 0.00 0.00 0.00 0.00
 [2,] 1.00 0.00 0.00 0.00 0.00
 [3,] 0.99 0.01 0.00 0.00 0.00
 [4,] 0.91 0.09 0.00 0.00 0.00
 [5,] 0.40 0.60 0.00 0.00 0.00
 [6,] 0.04 0.95 0.00 0.00 0.00
 [7,] 0.00 0.99 0.01 0.00 0.00
 [8,] 0.00 0.97 0.03 0.00 0.00
 [9,] 0.00 0.89 0.11 0.00 0.00
[10,] 0.00 0.69 0.31 0.00 0.00
[11,] 0.00 0.37 0.63 0.00 0.00
[12,] 0.00 0.14 0.85 0.01 0.00
[13,] 0.00 0.04 0.92 0.04 0.00
[14,] 0.00 0.01 0.89 0.10 0.00
[15,] 0.00 0.00 0.76 0.23 0.01
[16,] 0.00 0.00 0.54 0.43 0.03
[17,] 0.00 0.00 0.29 0.63 0.07
[18,] 0.00 0.00 0.12 0.72 0.15
[19,] 0.00 0.00 0.04 0.69 0.27
[20,] 0.00 0.00 0.01 0.57 0.42
[21,] 0.00 0.00 0.00 0.42 0.57
[22,] 0.00 0.00 0.00 0.28 0.72
[23,] 0.00 0.00 0.00 0.18 0.82
[24,] 0.00 0.00 0.00 0.10 0.90
[25,] 0.00 0.00 0.00 0.04 0.96

mixdist documentation built on May 2, 2019, 3:34 p.m.