Clara,
The most likely problem is that theta==0 for some of the
calculations. Or mu is NA, I suppose.
-Peter Ehlers
milton ruser wrote:> Dear Clara,
>
> Thanks for the reply. I am forwarding your message to the list, ok.
> When I wrote was a way of get further information to help the helpers.
>
> happy holidays,
>
> milton
>
> ---------- Forwarded message ----------
> From: Clara Br?ck <clara_brueck at web.de>
> Date: 2009/12/30
> Subject: Re: [R] Negbin Error Warnings
> To: milton ruser <milton.ruser at gmail.com>
>
>
>
>
> Dear Milton,
>
> Thanks for your email. I ran the NBR-model and then tried to compute the
> first differences of some
> expected values:
> m3.x1 <-setx(m3, advisory=1)
> m3.x2 <-setx(m3, advisory=0)
> m3.s1 <-sim (m3,x=m3.x1, x1=m3.x2)
>
> However, I then get 50 error warnings, and there all like this:
> In rpois(k, (mu * rgamma(k, theta))/theta) : NAs produced
>
>
> Bests
> Clara
>
>> -----Urspr?ngliche Nachricht-----
>> Von: "milton ruser" <milton.ruser at gmail.com>
>> Gesendet: 30.12.09 15:06:04
>> An: Clara Br?ck <clara_brueck at web.de>
>> CC: r-help at r-project.org
>> Betreff: Re: [R] Negbin Error Warnings
>
> Dear Clara,
>> could you, please, share the summary() of your data.
>>
>> bests
>>
>> milton
>>
>> 2009/12/30 Clara Br?ck <clara_brueck at web.de>
> > Hi,
>> I ran a negative binomial regression (NBR) using the Zelig-package
>> and the negbin model.
>> When I then try to use the simumlation approach using the setx () and
>> sim() functions to calculate expected values
>> and first difference for different levels of one of my independent
>> variables, I get 50 errors warnings, telling me that the calculation
>> rpois produced NAs. However, the data I use doesn't have any NAs.
>> Does this mean that the NBR-model is the wrong model for this data?
>>
>> Thanks in advance.
>>
>> ___________________________________________________________
>> Preisknaller: WEB.DE <http://web.de/> DSL Flatrate f?r nur 16,99
> Euro/mtl.!
>> http://produkte.web.de/go/02/
>>
>> ______________________________________________
>> R-help at 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-<http://www.r-project.org/posting->
>> guide.html
>> and provide commented, minimal, self-contained, reproducible code.
>>
>>
>>
>
>
> ___________________________________________________________
> Preisknaller: WEB.DE <http://web.de/> DSL Flatrate f?r nur 16,99
Euro/mtl.!
> http://produkte.web.de/go/02/
>
> [[alternative HTML version deleted]]
>
>
>
> ------------------------------------------------------------------------
>
> ______________________________________________
> R-help at 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 commented, minimal, self-contained, reproducible code.
--
Peter Ehlers
University of Calgary
403.202.3921