R/rqres.R

In lotze/COMPoissonReg: Conway-Maxwell Poisson (COM-Poisson) Regression

rqres.zicmp = function(y, lambda, nu, p, control = NULL)
{
	F = function(y) {
		pzicmp(y, lambda, nu, p, control = control)
	}
	rqres(y, F)
}

rqres.cmp = function(y, lambda, nu, control = NULL)
{
	F = function(y) {
		pcmp(y, lambda, nu, control = control)
	}
	rqres(y, F)
}

rqres = function(y, F, eps = 1e-6)
{
	n = length(y)
	FL = F(y - eps)
	FU = F(y)
	u = runif(n, min = FL, max = FU)
	qres = qnorm(u)
	return(qres)
}