getMixtureDensity: Get the mixture density

In RMixtCompUtilities: Utility Functions for 'MixtComp' Outputs

Get the mixture density

Description

Get the mixture density for each individual

Usage

getMixtureDensity(outMixtComp)

Arguments

outMixtComp

object of class MixtCompLearn or MixtComp obtained using mixtCompLearn or mixtCompPredict functions from RMixtComp package or rmcMultiRun from RMixtCompIO package.

Details

d(x_i) = \sum_k\pi_k P(x_i; \theta_k)

Value

a vector containing the mixture density for each individual.

Author(s)

Quentin Grimonprez

See Also

Other getter: getBIC(), getCompletedData(), getEmpiricTik(), getParam(), getPartition(), getType()

Examples

if (requireNamespace("RMixtCompIO", quietly = TRUE)) {
  dataLearn <- list(
    var1 = as.character(c(rnorm(50, -2, 0.8), rnorm(50, 2, 0.8))),
    var2 = as.character(c(rnorm(50, 2), rpois(50, 8)))
  )

  model <- list(
    var1 = list(type = "Gaussian", paramStr = ""),
    var2 = list(type = "Poisson", paramStr = "")
  )

  algo <- list(
    nClass = 2,
    nInd = 100,
    nbBurnInIter = 100,
    nbIter = 100,
    nbGibbsBurnInIter = 100,
    nbGibbsIter = 100,
    nInitPerClass = 3,
    nSemTry = 20,
    confidenceLevel = 0.95,
    ratioStableCriterion = 0.95,
    nStableCriterion = 10,
    mode = "learn"
  )

  resLearn <-RMixtCompIO::rmcMultiRun(algo, dataLearn, model, nRun = 3)

  d <- getMixtureDensity(resLearn)
}

RMixtCompUtilities documentation built on Sept. 22, 2023, 5:10 p.m.