How many em iterations are needed for convergence

Assignment Help Simulation in MATLAB
Reference no: EM131003294

Problem 1 [EM]

Write a 2-dimensional RNG for a Gaussian mixture model (GMM) pdf with 2 sub-populations. Use any function/sub-routine available in your language of choice.

Implement the expectation maximization (EM) algorithm for estimating the pdf parameters of 2-D GMMs from samples.

Compare the quality and speed your GMM-EM estimation on 300 samples of different GMM distributions (e.g. spherical vs ellipsoidal covariance, close vs well-separated subpopulations).

Problem 2 [Testing Faith]

Download the "old faithful" data set from blackboard. This has samples of a 2-D random variable: the first dimension is duration of the geyser eruption, the second is waiting time for the next eruption. Apply your GMM-EM algorithm to fit the data to a GMM pdf.

How many EM iterations are needed for convergence? Plot a contour plot of your final GMM pdf. Overlay the contour plot with a scatterplot of the data set. How would you use the GMM pdf estimates to cluster the data?

Problem 3 [Noise in GMM-EM]

Modify your GMM-EM routine by sampling and injecting Gaussian noise into the old faithful data at each iteration. Scale the noise to a fraction of the standard deviation in each dimension. And let the noise standard deviation decay at each iteration (e.g. inversely proportional to the square of the iteration counter).

Compare the average convergence time of the GMM-EM with and without noise. Plot the average convergence time for different initial noise standard deviations.

Turn in:

- A summary of your experiments including any relevant plots
- brief discussions of the results
- a print out of your code.

The following link gives a detailed tutorial about EM algorithms with close emphasis on EM for GMMs in section 3.2: https://mayagupta.org/publications/EMbookGuptaChen2010.pdf

You are not required to go through the whole paper in detail. But I recommend it.

Section 2.2 of this next paper also discusses the update equations for the GMM-EM algorithm: https://sipi.usc.edui-kosko/Noisy-Clustering-Neural-Networks.pdf

You can find more details about clustering with GMM-EM in the second paper linked above.

Reference no: EM131003294

Questions Cloud

Pick a health-related organization and study its website : There are many sites on the Internet that have information on healthcare and health information technology. Pick a health-related organization and study its website. Prepare a 1-2 paper describing the information and resources found on the site. E..
What the benefits of having a flexible budget : What the benefits of having a flexible budget? What are some of the downsides of having a flexible budget? Do you think the budget busters are accurate? Why or why not
Calculate the monthly specific returns : Calculate the monthly specific returns (i.e. the residuals of the regressions) for CBA, WES and BHP for the 84-month in-sample period -  Calculate the historical variance of the monthly specific returns for CBA, WES and BHP for the 84 month in-sampl..
Why is stock-based compensation added to net income : How is the undelivered portion of Microsoft's sales of Windows XP Professional recorded initially?
How many em iterations are needed for convergence : How many EM iterations are needed for convergence? Plot a contour plot of your final GMM pdf. Overlay the contour plot with a scatterplot of the data set.
Spring 2002 exam : Ten years ago, an advertising agency took a random sample of 20 personal computer owners who used brands such as Dell, Gateway, Hewlett-Packard, etc. At the time, the agency recorded a satisfaction score, with a maximum of 100 possible points.
Do you feel that the tools and technologies are appropriate : Then, review the tools and technologies posted by your classmates. Imagine that you were going to be tasked to assist them with their project. Based on their recommendations, do you feel that the tools and technologies are appropriate? Would you r..
Describe the steps involved in the innovation planning : Analyze the internal and external factors related to designing innovation strategies. Describe how they relate to one another and to the overall innovation process.
Find the function f for the given differential equation : Find the function f(x,y) whose differential is df = (x+y)-1 dx + (x+y)-1 dy and which has the value f(1,1)=0. Do this by performing a line integral on a rectangular path from (1,1) to (x1,y1) where x1>0 and y1>0.

Reviews

Write a Review

Simulation in MATLAB Questions & Answers

  What is the span of the columns

(a) Generate a random 3x5 matrix ?? with entry integers uniformly between -100 and 100. Hint: Search for help on the functions rand and randi. (b) Put the matrix ?? in reduced echelon form using the command rref. (c) Are the columns of A linearly ..

  Calculate the work done by integration

Takes an input argument and calculate the work done by integration usingn segment trapezoidal rule

  Write and implement a matlab program to apply 3 x 3 mask

Write and implement a matlab program to apply 3 x 3 mask

  Use matlab to plot the function

Plot the original periodic square wave on the same graph. Comment on the difference between the original periodic square wave and its truncated Fourier series presentation.

  Calculate the mid-point of the hermite curve

Calculate the mid-point of the Hermite curve that fits the points and the tangent vectors

  Voice transmission between bluetooth devices

Need someone to simulate a Bluetooth module in matlab showing that voice transmission between Bluetooth devices

  Solve the linear system

Compute the rank of S manually and use the MATLAB function to verify the result - Solve the linear system of equations Ax = b using MATLAB.

  Develop and simulate in matlab model

Develop and simulate in matlab model of trash to steamp lant in Baltimore

  Write down the expression for the laplace transform

Write down the expression for the Laplace Transform Y(s) of the solution. Output using a comment in your script. Obtain the solution y(t) using inverse Laplace transform in Symbolic Toolbox. Output using a comment in your script

  Queuing simulator

Queuing Simulator: Consider a communications router that can route (service) EXACTLY 800 messages per second (mps). Assume that messages arrive at an AVERAGE rate of 850 messages per second. Assume that the arriving messages follow an exponential dis..

  Generator and degradations - acquisition

Plot the autocorrelation function and show that other than peak there are only three values in the autocorrelation. What are the values of the autocorrelation function?

  Create a simulink model for the scara robot depicted

Create a Simulink Model for the SCARA robot depicted - The motor specification is appended to this brief - create a subsystem to model the motors taking into account that the inertia of Arms 1 and 2 combined is not a constant but depends upon the..

Free Assignment Quote

Assured A++ Grade

Get guaranteed satisfaction & time on delivery in every assignment order you paid with us! We ensure premium quality solution document along with free turntin report!

All rights reserved! Copyrights ©2019-2020 ExpertsMind IT Educational Pvt Ltd