Re: [R] Negative Binomial Regression - glm.nb

2013-02-28 Thread Ben Bolker
Martin Spindler Martin.Spindler at gmx.de writes: Dear all, I would like to ask, if there is a way to make the variance / dispersion parameter $\theta$ (referring to MASS, 4th edition, p. 206) in the function glm.nb dependent on the data, e.g. $1/ \theta = exp(x \beta)$ and to estimate

Re: [R] Negative Binomial Regression - glm.nb

2013-02-28 Thread Rolf Turner
On 02/28/2013 08:27 PM, Martin Spindler wrote: Dear all, I would like to ask, if there is a way to make the variance / dispersion parameter $\theta$ (referring to MASS, 4th edition, p. 206) in the function glm.nb dependent on the data, e.g. $1/ \theta = exp(x \beta)$ and to estimate the

[R] Negative Binomial Regression - glm.nb

2013-02-27 Thread Martin Spindler
Dear all, I would like to ask, if there is a way to make the variance / dispersion parameter $\theta$ (referring to MASS, 4th edition, p. 206) in the function glm.nb dependent on the data, e.g. $1/ \theta = exp(x \beta)$ and to estimate the parameter vector $\beta$ additionally. If this is

Re: [R] Negative Binomial Regression - glm.nb

2013-02-27 Thread Prof Brian Ripley
On 28/02/2013 07:27, Martin Spindler wrote: Dear all, I would like to ask, if there is a way to make the variance / dispersion parameter $\theta$ (referring to MASS, 4th edition, p. 206) in the function glm.nb dependent on the data, e.g. $1/ \theta = exp(x \beta)$ and to estimate the

[R] negative binomial regression with spatial weights matrix (not locations)

2011-07-19 Thread ChristyM
Hello, I have some data that exhibits a negative binomial distribution and also spatial structure. I would like to create a model that accounts for both. However, instead of locations, I have a distance matrix (cost matrix) describing the spatial relationships among the locations. I have tried

[R] negative binomial regression, unbalanced panel

2010-11-24 Thread hl
I am a student who is doing empirical work for his thesis and trying to switch to R. I am familiar with Stata, and at the moment I am trying to replicate some of my previous work. I have a large unbalanced panel data set, observations for different countries between 1970 and 2007. My

Re: [R] Negative Binomial Regression

2008-07-29 Thread Ben Bolker
jcarmichael jcarmichael314 at gmail.com writes: Hello. I am attempting to duplicate a negative binomial regression in R. SAS uses generalized estimating equations for model fitting in the GENMOD procedure. proc genmod data=mydata (where=(gender='F')); by agegroup; class id gender

Re: [R] Negative Binomial Regression

2008-07-29 Thread Prof Brian Ripley
On Tue, 29 Jul 2008, Ben Bolker wrote: jcarmichael jcarmichael314 at gmail.com writes: Hello. I am attempting to duplicate a negative binomial regression in R. SAS uses generalized estimating equations for model fitting in the GENMOD procedure. proc genmod data=mydata

Re: [R] Negative Binomial Regression

2008-07-29 Thread Ben Bolker
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Prof Brian Ripley wrote: | On Tue, 29 Jul 2008, Ben Bolker wrote: | | jcarmichael jcarmichael314 at gmail.com writes: | | | | Hello. | | I am attempting to duplicate a negative binomial regression in R. | SAS uses | generalized estimating equations

Re: [R] Negative Binomial Regression

2008-07-29 Thread Prof Brian Ripley
On Tue, 29 Jul 2008, Ben Bolker wrote: -BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Prof Brian Ripley wrote: | On Tue, 29 Jul 2008, Ben Bolker wrote: | | jcarmichael jcarmichael314 at gmail.com writes: | | Hello. | | I am attempting to duplicate a negative binomial regression in R. | SAS uses

Re: [R] Negative Binomial Regression

2008-07-29 Thread Ben Bolker
Prof Brian Ripley wrote: 'poisson' _family_, I presume? oops, yes. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide

[R] Negative Binomial Regression

2008-07-28 Thread jcarmichael
Hello. I am attempting to duplicate a negative binomial regression in R. SAS uses generalized estimating equations for model fitting in the GENMOD procedure. proc genmod data=mydata (where=(gender='F')); by agegroup; class id gender type; model count = var1 var2 var3 /dist=NB link=log