Re: [R] Great difference for piecewise linear function between R and SAS

2008-03-24 Thread Andrew Robinson
On Mon, Mar 24, 2008 at 01:09:35PM -0500, Paul Johnson wrote: > On Mon, Mar 24, 2008 at 8:11 AM, Duncan Murdoch <[EMAIL PROTECTED]> wrote: > > On 3/24/2008 9:03 AM, zhijie zhang wrote: > > > Dear Murdoch, > > > > > Thanks very > > > > I would use predict() instead. What you have there doesn't

Re: [R] Great difference for piecewise linear function between R and SAS

2008-03-24 Thread Paul Johnson
On Mon, Mar 24, 2008 at 8:11 AM, Duncan Murdoch <[EMAIL PROTECTED]> wrote: > On 3/24/2008 9:03 AM, zhijie zhang wrote: > > Dear Murdoch, > > > Thanks very > > I would use predict() instead. What you have there doesn't look as > though it uses the B-spline basis. > > The reference given in t

Re: [R] Great difference for piecewise linear function between R and SAS

2008-03-24 Thread Duncan Murdoch
On 3/24/2008 9:03 AM, zhijie zhang wrote: > Dear Murdoch, > Thanks very much for your rapid response. It helps me greatly. > If i want to write the model formula according to the estimets from > R, Is it correct for the below formula? I'm not very sure about it, and > i also hope that you c

Re: [R] Great difference for piecewise linear function between R and SAS

2008-03-24 Thread zhijie zhang
Dear Murdoch, Thanks very much for your rapid response. It helps me greatly. If i want to write the model formula according to the estimets from R, Is it correct for the below formula? I'm not very sure about it, and i also hope that you can recommend a good book on this area. I want to learn

Re: [R] Great difference for piecewise linear function between R and SAS

2008-03-24 Thread Duncan Murdoch
On 24/03/2008 7:06 AM, zhijie zhang wrote: > Dear Murdoch, >"Compare the predictions, not the coefficients.", this is the key point. >I have checked their predicted probability and their results are the > same. > What do you mean by "You're using different bases"? Could you please give >

Re: [R] Great difference for piecewise linear function between R and SAS

2008-03-24 Thread zhijie zhang
Dear Murdoch, "Compare the predictions, not the coefficients.", this is the key point. I have checked their predicted probability and their results are the same. What do you mean by "You're using different bases"? Could you please give me a little more info on it, so that i can go to find

Re: [R] Great difference for piecewise linear function between R and SAS

2008-03-24 Thread Duncan Murdoch
On 24/03/2008 5:23 AM, zhijie zhang wrote: > Dear Rusers, > I am now using R and SAS to fit the piecewise linear functions, and what > surprised me is that they have a great differrent result. See below. You're using different bases. Compare the predictions, not the coefficients. To see the

[R] Great difference for piecewise linear function between R and SAS

2008-03-24 Thread zhijie zhang
Dear Rusers, I am now using R and SAS to fit the piecewise linear functions, and what surprised me is that they have a great differrent result. See below. #R code--Knots for distance are 16.13 and 24, respectively, and Knots for y are -0.4357 and -0.3202 m.glm<-glm(mark~x+poly(elevation,2)+bs(di