glm.cmp
documentation.predict.cmp
that failed to handle new data with factors in it. autoplot
as an alias to gg_plot
. Credit to Emi Tanaka (@emitanaka) for this suggestion.broom
tidier
methods support. Specifically added method for tidy()
, glance()
and augment()
. summary()
was rewritten in order to support tidier
methods.vcov()
, confint()
, influence()
, hatvalues()
, rstandard()
, cooks.distance()
.glm.cmp()
no longer uses getnu()
. rcomp()
a bit by precalculating all the dcomp()
values. Credit to Guilherme Parreira (@guilhermeparreira) for the issue request. offset
term cannot be incorporated properly in the mean model. Credit to Sean Hardison (@seanhardison1) for finding this bug.model.matrix()
to extract model matrix from a fitted object. glm.cmp()
to allow varying dispersion. You can now link the dispersion parameters to some covariates via a log-link.fit_glm_cmp_const_nu()
and
fit_glm_cmp_vary_nu()
. print()
, summary()
are updated to support the updated glm.cmp()
. sitophilus
dataset to demonstrate the updated glm.cmp()
function. data(sitophilus)
M.sit <- glm.cmp(formula = ninsect ~ extract, formula_nu = ~extract,
data = sitophilus)
summary(M.sit)
gg_plot()
, gg_histcompPIT()
and gg_qqcompPIT()
are added to provide the ggplots version of the diagnostic plots. NEWS.md
file to track changes to the package. comp_mu_loglik_log_nu_only()
to facilitate the optimisation. Z()
, the normalizing constant function, approximates its true value via (a fixed) truncation. This means the approximation would fail if the mean is large.
The followings are implemented as a fix:
logZ()
is created, based on a similar function in the cmpreg
package of Ribeiro Jr, Zeviani & Demétrio (2019), and will supersede Z()
due to its superior numerical stability.glm.cmp()
, dcomp()
, pcomp()
, qcomp()
, rcomp()
and functions that calculate expected values are updated to take advantage of these changes.
sum(0:500*dcomp(0:500,100,1.2))
sum(0:500*dcomp(0:500,150,0.8))
qcomp(0.6, 150, 1.2)
fish
dataset as a proof of concept that glm.cmp()
can handle some larger count data. data(fish)
M.fish <- glm.cmp(species~ 1+log(area), data=fish)
max(M.fish$fitted.values)
comp_lambdas()
now has the ability to scale up \& down lambdaub
so that the correct $\lambda$s can be found even if they are outside the preset boundary. model.matrix()
now retrieves the design matrix of the model properly. glm.cmp()
gains a few standard glm arguments: start
, contrasts
, na.action
, subset
.Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.