I have an expression data. My hypothesis is that it is a combination of two models, negative binomial(NB) and Poisson. Currently, I'd like to do a simulation. Firstly I want to estimate the parameters of NB and Poisson using MLE in r. But I don't know how to create a negative log-likelihood function for two distribution. Anyone can help me?
Dear Liu,
from your description of data, it seems you could try a mixture model for your data. A mixture model assumes that your data don't arise from a single population, but rather from a linear combination of two or more subpopulations, each one with its own distribution and weight, under the constraint that all the weights must sum up to one. In your case, you could try to fit a Poisson mixture model and check whether two components are enough or more are needed. In R there are some packages that allow you to fit mixture models with ML methods or even with Bayesian methods. I'd suggest you to look at the package "Flexmix".
- Badges
- Science topic
- Similar topics
- Computer Science and Engineering
- Bioinformatics
- Bioinformatic Tools