Comments (6)
There seem to be experimental functions in broom.mixed, but these don't seem to work?
Moreover, I can't seem to find any other way to extract the residual degrees of freedom from the model (normally this would be given by glance()), is there any other way to extract this?
from glmmadaptive.
from glmmadaptive.
Thank you for the quick response. Is there any way to extract the residual df from a MixMod object?
from glmmadaptive.
from glmmadaptive.
Yeah, I ran into that too, maybe a temporary workaround can be added to allow functionality that depends on residual degrees of freedom?
For example from the glmmTMB package?
(Don't know if this is an acceptable solution, but this would allow the package to be used with for example MICE).
##' @importFrom stats df.residual
##' @method df.residual glmmTMB
##' @export
## TODO: not clear whether the residual df should be based
## on p=length(beta) or p=length(c(theta,beta)) ... but
## this is just to allow things like aods3::gof to work ...
## Taken from LME4, including the todo
##
df.residual.glmmTMB <- function(object, ...) {
nobs(object)-length(object$fit$par)
}
from glmmadaptive.
Perhaps it would allow it to work with MICE, but the relevant questions are where these degrees of freedom will be used and whether giving these degrees of freedom may lead to incorrect results.
from glmmadaptive.
Related Issues (20)
- emmeans functions in two-part models: only gives output for "main" model HOT 2
- Possibility to evaluate marginal log-likelihood HOT 2
- Error in chol.default(X[[i]], ...) HOT 1
- Random intercept only for zero-inflated portion, not count portion HOT 1
- Documentation contains dead links HOT 1
- May I know how to get the effect size of a two-part mixed effects model? HOT 4
- emmeans produces strange (wrong) estimates from beta regression after mixed_model HOT 1
- Adding "sandwich" argument to mixed_model HOT 1
- Issue with data prep for backward formulation HOT 4
- A problem with the censored.normal() family? HOT 2
- emmeans support doesn't return the link HOT 4
- Predict function for hurdle.lognormal response
- hurdle.lognormal() needs to trigger error message when zi_fixed is omitted HOT 1
- Confusion about the values returned by vcov(MixMod) HOT 4
- More than one random effect? HOT 2
- predict(type = 'subject_specific') not working in model with not positive definite Hessian matrix HOT 3
- "Error in eigen(Y, symmetric = TRUE): infinite or missing values in 'x'" HOT 1
- A binomial model with log link function doesn't display summary HOT 1
- post_vars Matrix Names HOT 1
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from glmmadaptive.