sigmaMCEM: E-Step estimate the parameter of sigma.

In JoaoPedro2536/univ.em:

Description

E-Step estimate the parameter of sigma.

Usage

sigmamcem(data, simzz, mupd)

Arguments

data	Contingency table, matrix format of three columns, first column lower limits, second column upper limits, and third column observed frequencies. The current state of the parameters, vector containing mean and sd.
simzz	The matrix, which is the samples simulated over the intervals, of the zmcem function.
mupd	Is the updated estimates of mu.

Value

Return sigma which are the estimates of variance in E-step.

Examples

 simdataaaa <- em.univ::univ_simul(ncol_matrix=1,
                         n=50,
                         nclass = 10,
                         mean = 68,
                         sd = 1.80,
                         fr_breaks=c(62,64,66,68,70,72,74,76,78))
 
 zm_matriz <- em.univ::zmcem(theta=c(67,2),data = simdataaaa$simul_data[,,1])
 
 mu_estimate <- em.univ::mumcem(data = simdataaaa$simul_data[,,1], simz = zm_matriz)
 
 sigma_estimate <- em.univ::sigmamcem(data = simdataaaa$simul_data[,,1],
                            simz = zm_matriz,
                            mupd = mu_estimate)
 sigma_estimate

JoaoPedro2536/univ.em documentation built on Dec. 18, 2021, 1:38 a.m.