[R] Poisson regression

2004-10-29 Thread Federico Gherardini
Dear all, First of all sorry if this is a dumb question... I'm trying to fit a linear model between the logarithm of two numerical variables (log(y) ~ log(x)). A log-log plot shows that the variance of log(y) is decreasing with the mean of log(x), in other words the points are quite dispersed

Re: [R] poisson regression with robust error variance ('eyestudy')

2004-06-09 Thread Lutz Ph. Breitling
okay, so now the bootcov works fine. aren't the lower bootstrap variances just what Karla is talking about when she writes on the website describing the eyestudy that i was trying to redo in the first place: Using a Poisson model without robust error variances will result in a confidence

Re: [R] poisson regression with robust error variance ('eyestudy')

2004-06-05 Thread Lutz Ph. Breitling
Thank you very much for your comments! however, i still do not get it right. robcov() accepts fit objects like lrm or ols objects as arguments, but obviously not the glmD objects (or at least not as simple as that). below some code to demonstrate. what am i still doing wrong? thx for your efforts-

Re: [R] poisson regression with robust error variance ('eyestudy')

2004-06-05 Thread Frank E Harrell Jr
Lutz Ph. Breitling wrote: Thank you very much for your comments! however, i still do not get it right. robcov() accepts fit objects like lrm or ols objects as arguments, but obviously not the glmD objects (or at least not as simple as that). below some code to demonstrate. what am i still doing

[R] poisson regression with robust error variance ('eyestudy')

2004-06-02 Thread Lutz Ph. Breitling
Dear all, i am trying to redo the 'eyestudy' analysis presented on the site http://www.ats.ucla.edu/stat/stata/faq/relative_risk.htm with R (1.9.0), with special interest in the section on relative risk estimation by poisson regression with robust error variance. so i guess rlm is the function

Re: [R] poisson regression with robust error variance ('eyestudy')

2004-06-02 Thread Thomas Lumley
On Wed, 2 Jun 2004, Lutz Ph. Breitling wrote: Dear all, i am trying to redo the 'eyestudy' analysis presented on the site http://www.ats.ucla.edu/stat/stata/faq/relative_risk.htm with R (1.9.0), with special interest in the section on relative risk estimation by poisson regression with

Re: [R] poisson regression with robust error variance ('eyestudy')

2004-06-02 Thread Frank E Harrell Jr
Thomas Lumley wrote: On Wed, 2 Jun 2004, Lutz Ph. Breitling wrote: Dear all, i am trying to redo the 'eyestudy' analysis presented on the site http://www.ats.ucla.edu/stat/stata/faq/relative_risk.htm with R (1.9.0), with special interest in the section on relative risk estimation by poisson

Re: [R] Poisson Regression

2003-03-28 Thread Martyn Plummer
On Fri, 2003-03-28 at 12:49, Meles MELES wrote: Hi, somebody asked me to do a Poisson regression on cancer incidence over years to see wether if there is a descending or an ascending tendancy. I tried with R, but it complains that the data are not integer but floating poing data. So,