similar to: Beta-coefficients for ZINB model

Dear all, I am modeling the incidence of recreational anglers along a stretch of coastline, and with a vary large proportion of zeros (>80%) have chosen to use a zero inflated negative binomial (ZINB) distribution. I am using the same variables for both parts of the model, can anyone help me with R code to compute overall marginal effects of each variable? My model is specified as follows:

Sorry bit of a Newbie question, and I promise I have searched the forum already, but I'm getting a bit desperate! I have over-dispersed, zero inflated data, with variance greater than the mean, suggesting Zero-Inflated Negative Binomial - which I attempted in R with the pscl package suggested on http://www.ats.ucla.edu/stat/R/dae/zinbreg.htm However my data is non-integer with some pesky

Dear all.. I have a respon variable y. Predictor variable are x1, x2, x3, x4, x5 (1) What is the syntax to get paramater estimation of ZINB Model by Newton Raphson (not BFGS) (2) What syntax to plot probability of observed & predicted of ZINB Thx. Regards Krist. [[alternative HTML version deleted]]

I want to generate 1,000 random samples of sample size=1,000 from ZINB. I know there is a rnegbin() to generate random samples from NB, and I know I can use the following process: do i=1 to 1000 n=0 do i=1 to 1000 if runi(1)>0.1 then x(i) = 0; else x(i)=rnegbin(); n=n+1; if n>1000 then stop; end; output; end; Anybody can help me out with the R code? Thanks very much ahead of time.

Hello, I am working with zero inflated models for a current project and I am getting wildly different results from R's zeroinfl(y ~ x, dist="negbin") command and Stata's zinb command. Does anyone know why this may be? I find it odd considering that zeroinfl(y ~ x, dist="poisson") gives identical to output to Stata's zip function. Thanks, --david [[alternative

Dear R User, Please, I am new to R. I want to overlay density plot for predictive interval pooled result in meta-analysis. http://addictedtor.free.fr/graphiques/graphcode.php?graph=114 Regards Frank Peter

Buenas tardes, Estoy estimando un modelo binomial negativo de ceros inflados (ZINB) utilizando el comando zeroinfl() del paquete pscl. Al ejecutarlo me da el siguiente aviso: Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred ¿Sabéis que significa y si puedo usar el modelo aún con ese aviso? ¿Los coeficientes son fiables? Muchas gracias, Miriam

Hi all, rmvnorm()can be used to generate the random numbers from a multivariate normal distribution with specified means and covariance matrix, but i want to specify the correlation matrix instead of covariance matrix for the multivariate normal distribution. Does anybody know how to generate the random numbers from a multivariate normal distribution with specified correlation matrix? What about

I get the following error when I try reading in a table. How are 1.1, 1.2, 1.3 duplicate row names? Thx. > table <- read.table('latestWithNumber.txt', header=T) Error in "row.names<-.data.frame"(`*tmp*`, value = c("1.1", "1.2", "1.3", : duplicate 'row.names' are not allowed Yongchuan

Hi, I am Olga Viedma. I am running a Zero-inflated negative binomial (ZINB) multi-level model using MCMCglmm package. I have a doubt. Can I use the "Liab" outputs as fitted data, instead of the predicted values from "predict"? The liab outputs fit very well with the observed data, whereas the predicted values are so bad. Thanks in advance, Olga Viedma D . Olga

Hi all, I have a problem to estimate Std. Error and t-value by ?polr? in library Mass. They result from the summary of a polr object. I can obtain them working in the R environment with the following statements: temp <- polr(formula = formula1, data = data1) coeff <- summary(temp), but when the above statements are enclosed in a function, summary reports the following error:

If I am buying a PC where the most compute intensive task will be running R and I do not have unlimited resources what trade-offs should I make? Specifically should I go for 1 - more memory, or 2 - faster processor, or 3 - something else? If it makes a difference I shall be running Windows on it and I am thinking about getting a portable which I understand makes upgrading more difficult. Extra

Muchas gracias por la respuesta. He mirado y los coeficientes no son altos pero sí tengo una gran cantidad de ceros en la variable dependiente (más del 90%). Sin embargo, al incluir otro tipo de variables independientes no me da ese aviso, dejando la misma variable dependiente. ¿Cómo podría utilizar stan/rstan de forma sencilla para diagnosticar el modelo? Muchas gracias El Lun, 9 de Abril de

Hi all, I have count data (number of flowering individuals plus total number of individuals) across 24 sites and 3 treatments (time since last burn). Following recommendations in the R Book, I used a glm with the model y~ burn, with y being two columns (flowering, not flowering) and burn the time (category) since burn. However, the residual deviance is roughly 10 times the number of degrees of

Hi, Anybody know what are the common Standard statistical methods used for the analysis of a dataset,and anybody know which of these methods give similar results Ram [[alternative HTML version deleted]]

I have been successfully been using R on my EEE PC for some while now. Although it uses Xandros I have been able to install various pieces of software by pretending that it was Debian etch. Johannes was kind enough to continue providing binaries on CRAN for a while but I think I am the only person still living in the past as CRAN etch is currently stuck at 2.11.0 and I cannot rely on him to

Hello! I have some data stored into 2 separate csv file. 1 file (called A.csv) (12 results named Group1, Group2, Group3, etc...) odds ratios, 2 file (called B.csv) 12 corresponded errors. How to import that data into R and make forest plot like I saw inside help file Rmeta and meta with included different font colors and names trough X and Y axis. I know for meta libb ... out <-

¿Quieres decir que para un nivel de una variable categorica todas las observaciones de la variable respuesta sean ceros? Gracias El Lun, 9 de Abril de 2018, 19:59, Carlos J. Gil Bellosta escribió: > ¿Podría ser que para algún nivel de alguna variable independiente > categórica solo hubiese ceros? En ese caso, casi seguro, aparecería ese > tipo de warning. > > El lun., 9 abr. 2018 a

In fact my first vignette full stop. I am intending to use Sweave. I have read the Sweave documentation and section 1.4 of the extensions manual and apart from (a) do not use split = TRUE (b) and include all the source components, there does not seem to be anything CRAN specific. 1 - Am I missing some other documentation? 2 - If I use a package from CTAN would that be considered within the

Dear Mr.Zeileis & all. (1)     Thx for your reply. Yes, I am talk about the function zeroinfl() from the package "pscl". I want to use Newton Raphson to get parameter             estimation ZINB, so I try this: ----------------------------------------------------------------------------------------------------------------------------------         > zinb <- zeroinfl(y