statdivlab / corncob

Count Regression for Correlated Observations with the Beta-binomial
102 stars 22 forks source link

Option for non-count estimator (estimating equations interpretation) #153

Closed adw96 closed 10 months ago

adw96 commented 10 months ago

Currently, corncob is set up to only take integer $W_ij$'s, which makes sense as it outputs MLEs from a beta-binomial model. However, we could consider a generalization for non-integer W_ij's, using the beta-binomial model to motivate estimating equations, which would still give us CANE estimates. Useful for modeling coverages! (and needed for a future project for which corncob is a dependency)

TODO: generalized corncob to estimate under non-integer W_ij's.

Must output a warning, because then people can throw anything in, without any guarantee that what they are doing is sensible, and previously this was often caught by the non-integer error.

Another option is to have a optional flag that defaults to FALSE allow_noninteger... maybe I like this idea better?

adw96 commented 10 months ago

Implementation of this currently on testing 1c3dd4a. Some notes

Still TODO

No need to review yet, @svteichman -- I'll open a PR when it's ready.

bryandmartin commented 10 months ago

This is a great new feature!

svteichman commented 10 months ago

@adw96 I did some reading and I think a good practice could be adding the authors of VGAM as authors in the DESCRIPTION file with the tag "cph" for copyright holder and then add a comment for the author saying that they are the creator of the VGAM package. Additionally, I think we should add a comment for the dbetabinom_cts_mod function stating that it is a modified version of VGAM::betabinom.ab. Here's an example of attributing package authorship like this in the shiny package.

svteichman commented 10 months ago

Or if copyright holder is not quite right in this situation (I'm not totally sure what is appropriate given the GPL-3 license for VGAM), they could instead by listed as a contributor with "ctb"

svteichman commented 10 months ago

Closing because this is implemented in merged pull request #155