To provide a latent variable model estimation apparatus capable of implementing the model selection at high speed even if the number of model candidates increases exponentially as the latent state number and the kind of the observation probability increase. A variational probability calculating unit 71 calculates a variational probability by maximizing a reference value that is defined as a lower bound of an approximation amount, in which Laplace approximation of a marginalized log likelihood function is performed with respect to an estimator for a complete variable. A model estimation unit 72 estimates an optimum latent variable model by estimating the kind and a parameter of the observation probability with respect to each latent state. A convergence determination unit 73 determines whether a reference value, which is used by the variational probability calculating unit 71 to calculate the variational probability, converges.