nm_expected: Calculate the expected probabilities

In mcomas/normalmultinomial: Normal-Multinomial Distribution

Description

Calculate the probabilities conditioned on parameters mu and sigma more likely to generate dataset X

Usage

nm_expected(X, mu, sigma, nsim = 1000, parallel.cluster = NULL)

Arguments

X	count data set
mu	paramater mu
sigma	parameter sigma
nsim	number of simulation used to calculated the expected values
parallel.cluster	instance of a Socket Cluster (see function makeCluster from package parallel)

Value

the expected probabilities

Examples

X = rnormalmultinomial(100, 100, rep(0,4), diag(4))
fit = nm_fit(X, verbose = T)
P = nm_expected(X, fit$mu, fit$sigma)
head(P)

mcomas/normalmultinomial documentation built on May 22, 2019, 3:15 p.m.