FlexDir (version 1.0)

FD.stddev: Standard Deviation of the ML estimators of a Flexible Dirichlet

Description

Conditional Bootstrap evaluation of the standard errors of the maximum likelihood parameter estimates of a Flexible Dirichlet distribution.

Usage

FD.stddev(x, iter.bootstrap = 500)

Arguments

x

an object of class FDfitted, usually the result of FD.estimation.

iter.bootstrap

number of iterations of the Bootstrap.

References

Ongaro, A. and Migliorati, S. (2013) A generalization of the Dirichlet distribution. Journal of Multivariate Analysis, 114, 412--426. Migliorati, S., Ongaro, A. and Monti, G. S. (2016) A structured Dirichlet mixture model for compositional data: inferential and applicative issues. Statistics and Computing, doi:10.1007/s11222-016-9665-y.

See Also

FD.estimation, FD.aicbic, FD.barycenters

Examples

data <- FD.generate(n=20,a=c(12,7,15),p=c(0.3,0.4,0.3),t=8)
data
results <- FD.estimation(data, normalize=TRUE,iter.initial.SEM = 5,iter.final.EM = 10)
results
FD.stddev(results)