Why do we use quasi-ML methods when estimating a quasi-Poisson regression? To make clear what I want to ask, I want to begin with the negative binomial distribution. The first two moments are $E(y)=\mu$ with variance $Var(y)=\mu+\mu^2/k$. When we solve the scorefunction of the log-likelihood with an iterationalgorithm, we get an estimator for $\widehat{\mu}$ (because of the estimator of $\widehat{\beta}$) and an estimator for the shape parameter $\widehat{k}$. So there are two parameters to be estimated.
When I look at the quasi-Poisson approach we have mean $E(y)=\mu$ and variance $Var(y)=\phi\cdot\mu$. This approach is not solvable since we do not know the probability function of $y$. Thus we use quasi-ML methods.
My question is: Why don't we create a distribution where we have the two moments of the quasi-Poisson regression and use a classic ML-estimation? It should be solvable with classic ML-methods since an estimation with negbin distribution (where we have also two estimators for $\mu$ and $k$ compared to $\mu$ and $\phi$).
I'm curious about your answers!
 A: You can't just "create a distribution" with a given form for the mean-variance relationship. In this case and others where quasi-likelihood is used, there is no probability distribution with the properties you desire. 
That's why it's called quasi-likelihood: 

"a quasi-likelihood function is not the log-likelihood corresponding to any actual probability distribution".

A: There are two ways here. The first way is to model the scaled Poisson data by the mixture of Normal distributions with parameters α1N(μ1,ϕ⋅μ1)+α2N(μ2,ϕ⋅μ2)+... You can easily write ML estimation for α, μ and ϕ. This approach works for relatively large mean values μ>10. The μ is small, you can use "scaled Poisson" distribution.
A: This is a great question. The issue with creating a distribution is that two moments are not sufficient to define a distribution. In fact, an infinite number of moments is required to perfectly define a single distribution. There is short discussion on this topic here on mathoverflow.
Simply put, there are an infinite number of distributions which have the desired first two moments. Therefore, which of the infinite number are appropriate to define for ML? There is no way to answer this question. Hence why Quasi-Likelihood methods are used instead.
