[R] poisson and peritonitis

Troels Ring tring at mail1.stofanet.dk
Mon Jan 22 22:47:56 CET 2001

Dear friends. I have a very simple question. In an old paper (Vonesh,Stat 
Med 1990;9:263-271), the rate of peritonitis in 81 peritoneal dialysis 
patients is estimated. The estimated intercept and an effect for race (0,1) 
in a generalized linear model with poisson distribution is given. The 
patients are observed a different length of time but time is not in the 
model, so I cannot see how is accounted for that. Anyhow, the observed 
number of zero incidences (35) is compared to the fitted zero incidences 
(said to be 33.54).
I would guess the way to do that was to let
mm <- exp(0.00149+race*0.7888) (his intercept and effect for race)
sum(dpois(0,mm)) but this doesn't yield but 25.64225

I know the times, but how could I get to the results of Vonesh without 
letting it have an effect ?

In hope of forgiveness if I have totally misunderstood what is going on


r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch

More information about the R-help mailing list