28 questions
- Bountied 0
- Unanswered
- Frequent
- Score
- Trending
- Week
- Month
- Unanswered (my tags)
1
vote
2
answers
102
views
easystats parameters package doesn't respect variable labels for pscl models
Consider this code:
options(parameters_labels = TRUE)
mydf <- pscl::bioChemists
labelled::var_label(mydf) <- list(art="MyCount", fem="MyGender",
...
0
votes
1
answer
25
views
pscl model with no zero-inflation
I have 3 outcome variables that I want to model with a negative binomial regression and 1 count variable I want to model with a zero-inflated negative binomial.
If I use MASS::glm.nb for the first ...
0
votes
1
answer
671
views
Interpretation of hurdle models with estimated marginal means
My goal is to interpret the coefficients of a hurdle model through estimated marginal means. I prefer to interpret probabilities (back-transformed from the logit scale), rather than log-odds (model ...
3
votes
1
answer
1k
views
`zeroinfl` from the `pscl` and `countreg` packages give very different results. Why?
I am experimenting with the mdvis dataset from the COUNT package of R for a teaching purpose. I fitted a zero-inflated negative-binomial model using the zeroinfl function from the pscl and countreg ...
2
votes
1
answer
317
views
Difference in results between glmmTMB and zeroinf (in pscl)
I applied glmmTMB and zeroinfl in pscl to the same dataset. I obtained the identical coefficients for the conditional part, but the coefficients for the binary part are somewhat different. Any idea on ...
0
votes
0
answers
281
views
Hurdle model: In sqrt(diag(object$vcov)) : NaNs produced with no data for some levels
Good afternoon, I am trying to run a hurdle model with two dependent variables: probability of disease persistence (0 or 1; the 'hurdle') and disease prevalence rate given there was disease ...
3
votes
1
answer
892
views
Is there a way to plot hurdle model results pscl package or a way to plot the count - zero truncated negbin part of hurdle model in r?
Is there a way to plot hurdle model results in R? I was able to plot the zero part (binomial with logit link) of the hurdle model (below) but I can't figure out how to plot the count part of the model ...
1
vote
0
answers
124
views
Back transforming pscl hurdle model output in R
Here's the setup of my hurdle model from the pscl package:
m1 <- hurdle(p_eaten ~ state + species + year, data = herb, dist = "negbin",
zero.dist = "binomial")
...
0
votes
1
answer
802
views
How to handle zero-inflated (semi-)continuous data in R?
I would like to model / fit Value on explanatory variables Type and Material (Value ~ Material + Type). Having a look at the sample test data provided here, one could see that Material X has all zero ...
0
votes
1
answer
177
views
Extracting z values from hurdle model output, 'pscl' package in R
I am using the hurdle() function from the pscl package to run negative binomial hurdle models in R.
I would like to extract only the list of z values from the model output and assign those values to a ...
2
votes
0
answers
162
views
how to generate confident intervals for predictions of hurdle models
I'm using 'pscl' package to build hurdle models, in order to model the abundance of birds as a function of environmental variables.
I am interested in predicting the probability of occurrence and ...
1
vote
0
answers
396
views
PSCL Returning NAs for zeroinfl negbin model
I am trying to run a Zero-Inflated Negative Binomial Count Model on some data containing the number of campaign visits by a politician by county. (Log Liklihood tests indicate Negative Binomial is ...
0
votes
0
answers
209
views
glmmadmb no longer supported in lsmeans. Alternative packages?
My data is zero inflated so I'm running a zero-inflated model using glmmamdb:
Model3z <- glmmadmb(Count3 ~ Light3 + (1|Site3), zeroInflation = T, family= "poisson", data = dframe3)
However, when ...
0
votes
0
answers
2k
views
How to run zero inflated poisson in R if zeroinfl function cannot be found?
I am using R version 3.6.0 and I am following this guide in order to run a zero-inflated poisson regression on cholera data: https://stats.idre.ucla.edu/r/dae/zip/
I have installed the pscl package ...
1
vote
1
answer
300
views
bias-corrected percentile bootstrap
I am trying to run a bias-corrected percentile bootstrap for 3 different effect sizes and ran across some code from a predecessor. While this code functions, the run time takes FAR too long. It takes ...