Re: [R] newbie problem using Design.rcs

2008-12-23 Thread Frank E Harrell Jr
you are doing this inside { } you will need to do print(Function(f)) Note this is a simplified version of the rcs with 2 redundant terms to avoid writing in terms of differences in cubes. Frank Sincerely, sp --- On Tue, 12/23/08, Frank E Harrell Jr wrote: From: Frank E Harrell Jr

Re: [R] newbie problem using Design.rcs

2008-12-23 Thread sp
/23/08, Frank E Harrell Jr wrote: > From: Frank E Harrell Jr > Subject: Re: [R] newbie problem using Design.rcs > To: to_rent_2...@yahoo.com > Cc: "David Winsemius" , r-help@r-project.org > Date: Tuesday, December 23, 2008, 4:57 PM > sp wrote: > > 2. I didn&#

Re: [R] newbie problem using Design.rcs

2008-12-23 Thread Frank E Harrell Jr
f this is NOT a good way to test fitting, could you pls tell me a better way? Respectfully, sp --- On Tue, 12/23/08, Frank E Harrell Jr wrote: From: Frank E Harrell Jr Subject: Re: [R] newbie problem using Design.rcs To: "David Winsemius" Cc: to_rent_2...@yahoo.com, r-help@r-project.org

Re: [R] newbie problem using Design.rcs

2008-12-23 Thread sp
t the regression curve differs from the actual points by as high as 10^5 with 3 knots and roughly -10^5 with 4 knots as I make y=x^2, y=x^3 If this is NOT a good way to test fitting, could you pls tell me a better way? Respectfully, sp --- On Tue, 12/23/08, Frank E Harrell Jr wrote: > Fr

Re: [R] newbie problem using Design.rcs

2008-12-23 Thread Frank E Harrell Jr
In addition to David's excellent response, I'll add that your problems seem to be statistical and not programming ones. I recommend that you spend a significant amount of time with a good regression text or course before using the software. Also, with Design you can find out the algebraic for

Re: [R] newbie problem using Design.rcs

2008-12-23 Thread David Winsemius
On Dec 22, 2008, at 11:38 PM, sp wrote: Hi, I read data from a file. I'm trying to understand how to use Design.rcs by using simple test data first. I use 1000 integer values (1,...,1000) for x (the predictor) with some noise (x+.02*x) and I set the response variable y=x. Then, I try rcs

[R] newbie problem using Design.rcs

2008-12-22 Thread sp
Hi, I read data from a file. I'm trying to understand how to use Design.rcs by using simple test data first. I use 1000 integer values (1,...,1000) for x (the predictor) with some noise (x+.02*x) and I set the response variable y=x. Then, I try rcs and ols as follows: m = ( sqrt(y1) ~ ( rcs(x1