Re: [R] GLM weights for the Poisson family

2014-02-06 Thread Milan Bouchet-Valat
I think you should have a look at svyglm() from the survey package. My two cents Le mercredi 05 février 2014 à 14:41 +1300, Rolf Turner a écrit : > You should direct your inquiry to R-help, not to me personally. I am > taking the liberty of cc-ing my reply back to the list. > > I really have

Re: [R] GLM weights for the Poisson family

2014-02-04 Thread Rolf Turner
You should direct your inquiry to R-help, not to me personally. I am taking the liberty of cc-ing my reply back to the list. I really haven't the time at the moment to think the issue through thoroughly, but off the top of my head: If you are going to use weighted log likelihoods then any

Re: [R] GLM weights for the Poisson family

2014-02-04 Thread Rolf Turner
On 04/02/14 20:12, IamRandom wrote: I am running a simple example of GLM. If I include weights when family="poisson" then the weights are calculated iteratively and $weights and $prior.weights return different values. The $prior.weights are what I supplied and $weights are the "posterior" wei

Re: [R] GLM weights for the Poisson family

2014-02-04 Thread David Winsemius
On Feb 3, 2014, at 11:12 PM, IamRandom wrote: > I am running a simple example of GLM. If I include weights when > family="poisson" then the weights are calculated iteratively and $weights and > $prior.weights return different values. The $prior.weights are what I > supplied and $weights are

[R] GLM weights for the Poisson family

2014-02-03 Thread IamRandom
I am running a simple example of GLM. If I include weights when family="poisson" then the weights are calculated iteratively and $weights and $prior.weights return different values. The $prior.weights are what I supplied and $weights are the "posterior" weights of the IWLS. If I include weig