123 questions
- Bountied 0
- Unanswered
- Frequent
- Score
- Trending
- Week
- Month
- Unanswered (my tags)
0
votes
1
answer
116
views
Error in EM algorithm: Non-finite Function Value in Likelihood Estimation with Integrals
I am encountering a persistent issue in my R program while performing Monte Carlo simulations for Maximum Likelihood Estimation (MLE). The simulation involves drawing samples from a modified Weibull ...
0
votes
0
answers
312
views
Manually program EM in r to update multiple parameters and solve missing data
I am trying to use EM (Expectation-maximization) to fill in missing data in R, but am not sure how to model/code it for my specific case. I am generally trying to follow the example format used in ...
0
votes
0
answers
226
views
Fitting a Gaussian Mixture Model with known share of noise/outliers
Question was moved to stats.stackexchange
A Gaussian Mixture model is fitted by the Expectation-Maximization algorithm.
This fairly simple algorithm consists of two steps and the initialization.
...
0
votes
1
answer
124
views
Computing Mixing Coefficients (Weights) of Mixed Copula Model (Gumbel and Unstructured Student-t) using EM Algorithm in R
Hello Stack Overflow community,
I am currently working on a project involving the estimation of a mixed copula model composed of a Gumbel copula and an unstructured Student-t copula. I would like to ...
0
votes
0
answers
49
views
How do I handle this error "Error in par[1] : object of type 'closure' is not subsettable"
I keep encountering error "Error in par[1] : object of type 'closure' is not subsettable" when I try integrating Aa. This happens after calling Aa(par,m)-code provided. I am trying to ...
1
vote
1
answer
359
views
EM algorithm for a mixture of three normal distributions throws errors
I need to run an EM algorithm for a mixture of three normal distributions with unknows means and variances. My data points are a column with 500 rows. I am gonna take it as 'S'.
First I need to write ...
1
vote
0
answers
120
views
Do Expectation Maximization always have to converge after a finite number of iterations if the implementation is correct?
I'm trying to understand how this could be possible, I know EM algorithm has the property to increase the likelihood for each step. However, this does not imply convergence. My question is, if the ...
0
votes
0
answers
78
views
EM algorithm for clustering of Gaussian Mixture models
I have a multivariate dataset(4 dimensions) of 150 points. I should use EM algorithm to group the given set of points into clusters. The number of clusters to which these points fall into is also ...
0
votes
1
answer
299
views
EM algorithm with Beta Distribution in R
I would like to write R code to build the dirichlet mixture model.
The loglikelihood I used for the beta distribution is as below:
LL(α,β)=(α−1)nlnxi ̄+(β−1)nln(1−xi) ̄+nlnΓ(α+β)−nlnΓ(α...
0
votes
1
answer
101
views
How can one add/remove components from a sklearn gmm object for further warm start?
The sklearn.mixture object GaussianMixture provides the framework to fit a GMM to provided data, but how can one add/remove components from a sklearn gmm object for further warm start?
3
votes
1
answer
1k
views
Singular Covariance Matrix in Expectation Maximization
I was trying to code up an EM algorithm in python for clustering images of different types. My understanding of the EM algorithm is as follows:
Accordingly, I coded the same in python. Here's my code:...
1
vote
1
answer
1k
views
Pgmpy: expectation maximization for bayesian networks parameter learning with missing data
I'm trying to use the PGMPY package for python to learn the parameters of a bayesian network. If I understand expectation maximization correctly, it should be able to deal with missing values. I am ...
1
vote
1
answer
405
views
Expectation Maximization using a Poisson likelihood function
I am trying to apply the expectation-maximization algorithm to estimate missing count data but all the packages in R, such as missMethods, assume a multivariate Gaussian distribution. How would I ...
2
votes
0
answers
367
views
sketching a Gaussian Mixture plot
So basically I have these points and weights and covariance matrices and 1st mean=x1 and 2nd mean=x2.
What I am trying to do is to plot the clusters of that Gaussian distribution with ellipses ...
0
votes
1
answer
54
views
Defining the EM parameters
I have a list of observations where each data point is a pair of a time expression (e.g. night, morning) and an hour in a 12-hr clock (i.e. 1, 2, ..., 12): Y = {<e_i, h_i>}_i={1,...,N}. I would ...