similar to: minor quirk in rnbinom

Just for the records (was on R-help). Uwe Ligges -------- Original Message -------- Subject: Re: [R] minor quirk in rnbinom Date: Fri, 14 Jun 2002 18:47:59 +0200 From: Uwe Ligges <ligges@statistik.uni-dortmund.de> To: Mike Lonergan <mel@mcs.st-and.ac.uk> CC: r-help <r-help@stat.math.ethz.ch> References: <NEBBJLLBCLEDFLCCBKDAEEEKCGAA.mel@mcs.st-andrews.ac.uk> Mike

I am trying to reset the default arguments in the rnbinom function with the following example code: params <- c("size" = 1, "mu" = 1) formals(rnbinom)[names(params)] <- params rnbinom(n = 10) It returns the following: Error in rnbinom(n = 10) : argument "prob" is missing, with no default If I set the defaults with this code: params <- c("size"

The documentation states : An alternative parametrization (often used in ecology) is by the mean ?mu?, and ?size?, the dispersion parameter. However, this fails : > rnbinom(10, mu = 100, size = 0) [1] NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN Warning message: In rnbinom(10, mu = 100, size = 0) : NAs produced For dispersion set to 0, it should work like drawing from a Poisson distribution.

wanfairos wrote: > > Dear members; > > I have one problem to display my simulated data. Basically, I've generated > a set of rnorm (std.norm.mat)and rnbinom (gen.data) data to test on the > performance of FE Poisson model estimates > > ... Long self-contained example shorted > > my.data<-rep(9999,sim.size) > > for(i in 1: sim.size){ > ... >

Dear R-helpers, let us assume, that I have the following dataset: a <- rnbinom(200, 1, 0.5) b <- (1:200) c <- (30:229) d <- rep(c("q", "r", "s", "t"), rep(50,4)) data_frame <- data.frame(a,b,c,d) In a first step I run a glm.nb (full code is given at the end of this mail) and want to predict my response variable a. In a second step, I would

(Sent this to r-help back in October, got no comments, forgot to re-submit it as a bug report.) There are a few inconsistencies, at least, in some of the functions that generate random deviates from particular distributions (I think they're bugs because they're inconvenient, but maybe someone can make an argument for the current behavior). If people think these are really bugs I can

Dear list, I ran into a strange behaviour of the pnbinom function - or maybe I just made a stupid mistake. First thing is that pnbinom seems to be very slow. The other - more interesting one - is that I get two different curves when I plot the estimated density and the density given by pnbinom. Shouldn't it be the same? This is only the case, I think, if I use the parameter size = 1. I

Hi R-Users! I fitted a negative binomial distribution to my count data using the function glm.nb() and obtained the calculated parameters theta (dispersion) and mu. I would like to simulate values from this negative binomial distribution. Looking at the function rnbinom() I was looking at the relationship between the two possible parametrizations of the negative binomial and found that for this

Hi, The code below is giving me this error message: Error in while (err > eps) { : missing value where TRUE/FALSE needed In addition: Warning messages: 1: In dnbinom(x, size, prob, log) : NaNs produced 2: In dnbinom(x, size, prob, log) : NaNs produced I know from the help files that for dnbinom "Invalid size or prob will result in return value NaN, with a warning", but I am not able

How can I add an overall plot title to these four plots? I would like to have something that says, "Distribution Comparisons": par(mfrow = c(2, 2)) # Plot 1 plot(rnorm(10),type="l",col="red") title(main = list(paste("Normal"), ????????????????????????? col="black", cex = 1.0)) # Plot 2 plot(rpois(10,

Hello I want to use the nlminb function but I have the objective function like characters. I can summarize the problem using the first example in the nlminb documentation. x <- rnbinom(100, mu = 10, size = 10) hdev <- function(par) -sum(dnbinom(x, mu = par[1], size = par[2], log = TRUE)) nlminb(c(9, 12), objective=hdev) With the last instructions we obtain appropriate results. If I have

Hi, I wonder if predict.gam is supposed to work with family=negbin() definition? It seems to me that the values returned by type="response" are far off the observed values. Here is an example output from the negbin examples: > set.seed(3) > n<-400 > dat<-gamSim(1,n=n) > g<-exp(dat$f/5) > dat$y<-rnbinom(g,size=3,mu=g) >

Hi, dear R experts , I am new. I met this problem when I am trying to learn how to use the nlminb() function. I tried the example which the document provides ( as the following code ) and R gives no response . I don't know whether it is running or not and it takes a very long time but still output nothing so I just close the session window. my question is is there any method

Hello, Sample code: library("vcd") dummy <- rnbinom(200, size=1.5, prob=0.8) gf <- goodfit(dummy, type="nbinomial", method="MinChisq") plot(gf) I would like to: 1. add a lgened stating the bars show the actual counts and the red dots - the fit. 2. show the goodness-of-fit values calculated somewhere on an empty white space ob the plot. But... the legend

Can anyone explain to me why the AIC values are so different when using glm.nb and glm with a negative.binomial family, from the MASS library? I'm using R 1.8.1 with Mac 0S 10.3.4. >library(MASS) > dfr <- data.frame(c=rnbinom(100,size=2,mu=rep(c(10,20,100,1000),rep(25,4))), + f=factor(rep(seq(1,4),rep(25,4)))) > AIC(nb1 <- glm.nb(c~f, data=dfr)) [1] 1047 >

Dear list members, while I appreciate the possibility to deal with overdispersion for count data either by specifying the family argument to be quasipoisson() or negative.binomial(), it estimates just one overdispersion parameter for the entire data set. In my applications I often would like the estimate for overdispersion to depend on the covariates in the same manner as the mean. For example,

Dear Experts, Probably trivial, but I am struggling to get what I want: I need to know how the number of required trials to get a certain number of successes. By example: How many trials do I need to have 98% probability of 50 successes, when the a priory probability is 0.1 per trial. The Negative binomial function may do the job (not sure): NegBinomial {stats} The Negative Binomial

Dear R users I would like to simulate underdispersed Poisson and binomial distributions somehow. I know you can do this for overdispersed counterparts - using rnbinom() for Poisson and rbetabinom() for binomial. Could anyone share functions to do this? Or please share some tips for modifying existing functions to achieve this. Thank you very much for your help and time Shinichi

On 1/20/20 12:33 PM, Martin Maechler wrote: > > It's really something that should be discussed (possibly not > here, .. but then I've started it here ...). > > The NEWS for R 3.0.0 contain (in NEW FEATURES) : > > * Functions rbinom(), rgeom(), rhyper(), rpois(), rnbinom(), > rsignrank() and rwilcox() now return integer (not double) > vectors.

Ugh, sounds like competing priorities. * maintain type consistency * minimize storage (= current version, since 3.0.0) * maximize utility for large lambda (= proposed change) * keep user interface, and code, simple (e.g., it would be easy enough to add a switch that provided user control of int vs double return value) * backward compatibility On 2020-01-20 12:33 p.m., Martin Maechler