similar to: RE : Dispersion parameter in Neg Bin GLM

Dear Lister, I am facing a strange problem fitting a GLM of the negative binomial family. Actually, I tried to estimate theta (the scale parameter) through glm.nb from MASS and could get convergence only relaxing the convergence tolerance to 1e-3. With warning messages: glm1<-glm.nb(nbcas~.,data=zonesdb4,control=glm.control(epsilon = 1e-3)) There were 25 warnings (use warnings() to see

Hali! Mi a legegyszer?bb m?dja annak, hogy egy pl. note = Notes.create(...) ut?n megkapjam az adott note ActiveRecord id-j?t (primary key)? Sajnos ez ugye 0, mivel frissen lett l?trehozva... ?tlet? G?bor

I have tried glm.nb in the MASS package, but many models (I have 250 models with different combinations of predictors for fish counts data) either fail to converge or even diverge. I'm attempting to use the negbin function in the AOD package, but am unsure what to use for the "random" term, which is supposed to provide a right hand formula for the overdispersion parameter.

Hi, can anyone help me fit betabinomial model to the following dataset where each iD is a cluster in itself , if i use package aod 's betabinom model it gives an estimate of zero to phi(the correlation coeficient ) and if i fix it to the anova type estimate obtained from icc( in package aod) then it says system is exactly singular. And when i try to fit my loglikelihood by

Hi, I have sort of a newbie question. I've seriously put a lot of effort into how to handle simple replicates in a repeated ANOVA design, but haven't had much luck. I really liked reading "Notes on the use of R for psychology experiments and questionnaires", by Jonathan Baron and Yuelin Li ( http://www.psych.upenn.edu/~baron/rpsych/rpsych.html ) but still didn't run across

Hi, I try to calculate the correlation between macroeconomic data from FRED vs Market Data However, since the timeseries are not in synch, the correlation fails. require(quantmod) USPBS =get(getSymbols("USPBS", src="FRED" )) USPBS = USPBS['1983-1-1::'] monDMANEMP = Cl(to.monthly(USPBS)) > length(monDMANEMP) [1] 312 > head(monDMANEMP) USPBS.Close

Hello, I am experiencing trouble with r-term 1.1.1 for windows 2000. When I start rterm, and I get a prompt, it seems to freeze after any assignment statement (so, if I type x_rnorm(1000), I don't get another prompt). Statements that don't involve assignment statements (like "plot(rnorm(1000))") seem to work, and I get a new prompt. The Rgui seems to work fine, and if I *can*

Dear all, I have a problem with an Rd file containing French accentuated characters. I have uploaded the file at http://filex.cirad.fr/get?k=cjW7lImMaNC6Ci2vX0H I have declared Encoding: latin1 in the package DESCRIPTION file and I have added \encoding{latin1} in the header of the Rd file. When I compile the package manual, I have LaTeX errors: ! Package inputenc Error: Unicode char \u8:?F

Hi all, I've done some anovas, found some significance and proceeded to do the post hoc tests to determine the source of the significance. I used the multcmp package as suggested by Andy Liaw (thank you). I am, however, unsure how to interpret the results of the call to simint. Could someone help me with this? I've attached the code and the results it produced. There are 3 groups

Hi there, I have relative abundance data for 13 mammal species that I collected at various sites that ranged in road density. I'm trying to determine the effect of road density on animal abundance across body sizes. For most species, I have data that was collected in one year but for a few species I have two years of complete data, and would like to use both. Since I have count data, I'm

Dear all, I have a puzzle regarding the estimation of Neg. Binomial event count model in R. I would greatly appreciate if anyone could shed some light on my puzzle. Using the glm.nb command, or the zelig command developed by Gary King et. al., I obtain the same point estimates in R as well as in Stata. However, if I write my own likelihood function to estimate a neg. binomial event count

Hi, I want to fit a GEE with a negative binomial distribution. I have uesd already a poisson glm and then neg binommial to deal with alot of dispersion. In my neg binomial residuals i have some patterns so i have implemented a GEE, but only with a poisson family as i couldnt with neg binomial. However the residual patterns in fact look worse here. When i try and put neg binomial family it

I am using glm.nb, a ~b*c ( b is categorical and c is continuous). when I run this model I get the warning message: Warning messages: 1: In theta.ml(Y, mu, sum(w), w, limit = control$maxit, trace = control$trace > : iteration limit reached 2: In theta.ml(Y, mu, sum(w), w, limit = control$maxit, trace = control$trace > : iteration limit reached What does this mean? -- Graduate

Fixes: https://bugs.freedesktop.org/show_bug.cgi?id=86618 Signed-off-by: Ilia Mirkin <imirkin at alum.mit.edu> --- src/gallium/drivers/nouveau/codegen/nv50_ir_emit_nv50.cpp | 2 ++ 1 file changed, 2 insertions(+) diff --git a/src/gallium/drivers/nouveau/codegen/nv50_ir_emit_nv50.cpp b/src/gallium/drivers/nouveau/codegen/nv50_ir_emit_nv50.cpp index 077eba8..3048f3d 100644 ---

Hi, Can someone tell me if it is possible to set the dispersion parameter constant when fitting a negative binomial glm in R? I've looked at the documentation and can't find the appropriate argument to pass. In STATA I can type: nbreg depvar [indepvar...], offset(offset) dispersion(constant). Thank you [[alternative HTML version deleted]]

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 This release provides an important security fix outlined in the release notes that follow. This is the latest stable release of Samba and the version that all production Samba servers should be running for all current bug-fixes. The source code can be downloaded from : http://download.samba.org/samba/ftp/ in the file samba-2.2.8a.tar.gz or

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 This release provides an important security fix outlined in the release notes that follow. This is the latest stable release of Samba and the version that all production Samba servers should be running for all current bug-fixes. The source code can be downloaded from : http://download.samba.org/samba/ftp/ in the file samba-2.2.8a.tar.gz or

Dear list, I'm currently trying to develop a model to assess clam yield potential in a lagoon. I'm using the zeroinfl function of the pscl package to fit a Zero-inflated negative binomial model, given the high occurrence of zero counts. I don't understand from the sentence in the pscl guide "Zero-inflated count models are a type of two-component mixture model, with a component

Hi Everyone: I'm going a little nuts here and am hoping someone might have some ideas to help out. Here is my problem: I am using the calendarHeatMap function (http://blog.revolutionanalytics.com/2009/11/charting-time-series-as-calendar-heat-maps-in-r.html) to plot some values of percentages above or below a watermark. In other words, I have a time series whose data can range arbitrarily

https://bugs.freedesktop.org/show_bug.cgi?id=86618 Bug ID: 86618 Summary: [NV96] neg modifiers not working in MIN and MAX operations Product: Mesa Version: git Hardware: Other OS: All Status: NEW Severity: normal Priority: medium Component: Drivers/DRI/nouveau