R/lambda_m_step.R

In MixedPoisson: Mixed Poisson Models

Defines functions lambda_m_step

Documented in lambda_m_step

lambda_m_step <-function(variable, X, offset){
g=glm(variable~-1+X, offset=offset,  family=poisson(log))
lambda.new=exp(X%*%g$coefficients) 

outlist = list(lambda=lambda.new, beta=g$coefficients, glm=g)
    class(outlist) = "lambda_m_step"
    return(outlist)
}

Try the MixedPoisson package in your browser