[R] poisson and peritonitis
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