Re: [R] Consistency of Logistic Regression

2010-11-12 Thread Benjamin Godlove
n the model? it could be > just a parameterization difference. > > albyn > > On Thu, Nov 11, 2010 at 12:41:03PM -0500, Benjamin Godlove wrote: > > Dear R developers, > > > > I have noticed a discrepancy between the coefficients returned by R's > glm() > &g

[R] Consistency of Logistic Regression

2010-11-11 Thread Benjamin Godlove
Dear R developers, I have noticed a discrepancy between the coefficients returned by R's glm() for logistic regression and SAS's PROC LOGISTIC. I am using dist = binomial and link = logit for both R and SAS. I believe R uses IRLS whereas SAS uses Fisher's scoring, but the difference is something

[R] feedback/question on function update()

2010-09-18 Thread Benjamin Godlove
Hi, First let me say I am a big fan of R and appreciate all your time and effort. The update() function does not seem to work in a for loop. Consider the following: mdat <- matrix(c(1,2,3, 11,23,13, 12,4,8), nrow = 3, ncol=3, byrow=TRUE) reg <- lm(mdat[7:9]~1) for(i in 1:2) { reg <- update(reg,