# GMM_fitting_graph: Visualization of the result of spect_em_gmm In EMpeaksR: Conducting the Peak Fitting Based on the EM Algorithm

## Description

Visualization of the result of spect_em_gmm().

## Usage

 `1` ```show_gmm_curve(spect_em_gmm_res, x, y, mix_ratio_init, mu_init, sigma_init) ```

## Arguments

 `spect_em_gmm_res` data set obtained by spect_em_gmm() `x` measurement steps `y` intensity `mix_ratio_init` initial values of the mixture ratio of the components `mu_init` initial values of the mean of the components `sigma_init` initial values of the standard deviation of the components

## Details

Perform a visualization of fitting curve estimated by Gaussian mixture model.

## Value

Show the fitting curve and variation of the parameters.

## References

Matsumura, T., Nagamura, N., Akaho, S., Nagata, K., & Ando, Y. (2019). Spectrum adapted expectation-maximization algorithm for high-throughput peak shift analysis. Science and technology of advanced materials, 20(1), 733-745.

## Examples

 ``` 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32``` ```#generating the synthetic spectral data based on three component Gausian mixture model. x <- seq(0, 100, by = 0.5) true_mu <- c(35, 50, 65) true_sigma <- c(3, 3, 3) true_mix_ratio <- rep(1/3, 3) degree <- 4 y <- c(true_mix_ratio[1] * dnorm(x = x, mean = true_mu[1], sd = true_sigma[1])*10^degree + true_mix_ratio[2] * dnorm(x = x, mean = true_mu[2], sd = true_sigma[2])*10^degree + true_mix_ratio[3] * dnorm(x = x, mean = true_mu[3], sd = true_sigma[3])*10^degree) plot(y~x, main = "genrated synthetic spectral data") #Peak fitting by EMpeaksR #Initial values K <- 3 mix_ratio_init <- c(0.2, 0.4, 0.4) mu_init <- c(20, 40, 70) sigma_init <- c(2, 5, 4) #Coducting calculation SP_EM_G_res <- spect_em_gmm(x, y, mu = mu_init, sigma = sigma_init, mix_ratio = mix_ratio_init, conv.cri = 1e-8, maxit = 100000) #Plot fitting curve and trace plot of parameters show_gmm_curve(SP_EM_G_res, x, y, mix_ratio_init, mu_init, sigma_init) #Showing the result of spect_em_gmm() print(cbind(c(mu_init), c(sigma_init), c(mix_ratio_init))) print(cbind(SP_EM_G_res\$mu, SP_EM_G_res\$sigma, SP_EM_G_res\$mix_ratio)) print(cbind(true_mu, true_sigma, true_mix_ratio)) ```

EMpeaksR documentation built on April 20, 2021, 5:10 p.m.