fit.nbinom {semisup} | R Documentation |
This function fits the semi-supervised negative binomial mixture model.
It is called by fit.wrap
.
fit.nbinom(y, z, phi, gamma, it.em, epsilon)
y |
observations:
numeric vector of length |
z |
class labels:
integer vector of length |
phi |
dispersion parameters:
numeric vector of length |
gamma |
offset:
numeric vector of length |
it.em |
(maximum) number of iterations in the |
epsilon |
convergence criterion for the |
This function returns the parameter estimates, the posterior probabilities, and the likelihood.
This is an internal
function.
The user functions are mixtura
and scrutor
.
# data simulation n <- 100 z <- rep(0:1,each=n/2) gamma <- runif(n=n,min=0,max=2) y <- rnbinom(n=n,mu=gamma*(5+2*z),size=1/0.05) z[(n/4):n] <- NA # model fitting fit.nbinom(y,z,phi=0.05,gamma=gamma, it.em=100,epsilon=1e-04)