SSest: E-Step For Variance Estimate

In JoaoPedro2536/univ.em:

Description

E-Step For Variance Estimate

Usage

ssest(theta, bl, bu, muupdate, freq)

Arguments

theta	Initial theta value, vector containing mean and sd.
bl	Lower bound of the intervals, values in vector, starting from -inf.
bu	Upper bound of the intervals, values in vector, ending with +inf.
muupdate	Is the updated estimates of mu.
freq	Frequency over the intervals, values in vector.

Value

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

Examples

 library(em.univ)
 library(stats)


 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))
 
 ssnew <- em.univ::ssest(theta=c(67,2),
              bl=simdataaaa$simul_data[,1,1],
              bu=simdataaaa$simul_data[,2,1],
              muupdate=mest(theta=c(67,2),
                            bl=simdataaaa$simul_data[,1,1],
                            bu=simdataaaa$simul_data[,2,1],
                            freq=simdataaaa$simul_data[,3,1]),
              freq=simdataaaa$simul_data[,3,1])
 
 ssnew

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