mu_zmcem_mult: Estimate Mu for Sample

In JoaoPedro2536/univ.em:

Description

This is the function for calculation mean vector from the generated samples.

Usage

mu_zmcem_mult(data, generated_sample)

Arguments

data	Data in the form of multivariated grouped data.
generated_sample	Array, which is the simulated (generated) data.

Value

Return the vector of means.

Examples

library(MASS)
library(em.univ)
library(tmvtnorm)

simulateddata = em.univ::mult_simul(mm = c(68,68),
                          ss = base::matrix(c(3,2,2,6),2,2) ,
                          n_data_sets = 1,
                          breaks_x = c(-Inf,64,65,66,67,68,69,70,71,72,Inf),
                          breaks_y = c(-Inf,64.2,65.2,66.2,
                                       67.2,68.2,69.2,70.2,
                                       71.2,72.2,Inf),
                          lower_x = base::rep(c(-Inf,64,65,66,67,
                                          68,69,70,71,72),10),
                          lower_y = c(base::rep(-Inf,10),
                                      base::rep(64.2,10),
                                      base::rep(65.2,10),
                                      base::rep(66.2,10),
                                      base::rep(67.2,10),
                                      base::rep(68.2,10),
                                      base::rep(69.2,10),
                                      base::rep(70.2,10),
                                      base::rep(71.2,10),
                                      base::rep(72.2,10)),
                          upper_x = base::rep(c(64,65,66,67,68,69,70,71,72,Inf),10),
                          upper_y = c(base::rep(64.2,10),
                                      base::rep(65.2,10),
                                      base::rep(66.2,10),
                                      base::rep(67.2,10),
                                      base::rep(68.2,10),
                                      base::rep(69.2,10),
                                      base::rep(70.2,10),
                                      base::rep(71.2,10),
                                      base::rep(72,2,10),
                                      base::rep(Inf,10))
)

mu <- c(67,67)
sigma <- matrix(c(3.1, 2.16, 2.16, 6.05), 2, 2)

generated_sample<- em.univ::zmcem_mult(data=simulateddata[,,1],
                              mu=mu,
                              sigma=sigma)

out<- em.univ::mu_zmcem_mult(data=simulateddata[,,1],
                     generated_sample=generated_sample)


out

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