Re: Multiple Regression w/ Polynomial-in-Y?

From: Paul Victor Birke (nonlinear_at_rogers.com)
Date: 09/25/04


Date: Sat, 25 Sep 2004 09:27:11 -0400
To: Frank Iannarilli <frankeye@cox.net>

Dear Frank

Firstly, I don't know what this is called. I was considering a model
such as C1 * ln(y) + C2 * exp(y) just the other day so it was with some
interest I saw your post. What did Karl Jung say about synchronicity??!!

What I would like to know is why you were movivated to do this.
You mentioned in your text the feeling or belief the RHS was nonlinear.
 From an Information-Theoretic perspective would we believe this
polynomial formulation would be somehow better.

Why did you suppress w0 the constant. Unless your formulation is a
probability summation I think you should not do that. Then w0 estimates
the error term I think

Anyways, well done and I hope we hear from some math gurus here

Paul Birke, P. Eng.

Frank Iannarilli wrote:
> Hi,
>
> Is what I'm tempted to call "Polynomial Root Regression" so obvious
> that nobody talks/writes about it?
>
> Here's the model, which is readily solved with multiple regression
> methods (also tried it successfully with shrinkage regression
> technique such as Partial Least Squares) - note well, the coefficients
> of the dependent-variable polynomial are unknown, and are estimated by
> the (standard) regression (in addition to the usual weights for the
> independent variables).
>
> Here, a quadratic-in-Y polynomial (extension to higher-order
> obvious):
>
> x.w = y + c*y^2 + error
>
> M indep. vars; N samples; (intercept term suppressed)
> Both w vector and c coefficient unknown/to be estimated
>
>
> |x(1,1), x(1,2)...x(1,M); -y(1)^2| |w(1)| |y(1)|
> |x(2,1), x(2,2)...x(2,M); -y(2)^2| |w(2)| |y(2)|
> |....... | |... | = |... |
> |x(N,1), x(N,2)...x(N,M); -y(N)^2| |w(M)| |y(N)|
> | c |
>
> Solve above with appropriate linear least-squares solver (e.g., OLS,
> or RR or PLS if multicollinearity). Yields w-vector and c coefficient
> for quadratic term for y.
>
> To predict a new y-value, just find the root y of the determined
> polynomial
>
> c*y^2 + y - x.w = 0
>
>
> I tried this successfully for a problem whose non-linearity was best
> modeled on the dependent-variable side, rather than forcing polynomial
> terms of independent variables (which I suspect would not be so good
> for my application). One drawback is needing to decide which root is
> the actual solution, but for many situations this is probably "easy".
>
> Thoughts? Does this have a name?
>
> Thanks!



Relevant Pages

  • Re: Multiple Regression w/ Polynomial-in-Y?
    ... > Both w vector and c coefficient unknown/to be estimated ... Yields w-vector and c coefficient ... > terms of independent variables (which I suspect would not be so good ... One drawback is needing to decide which root is ...
    (sci.math.num-analysis)
  • Re: Multiple Regression w/ Polynomial-in-Y?
    ... Yields w-vector and c coefficient ... > for quadratic term for y. ... > terms of independent variables (which I suspect would not be so good ... One drawback is needing to decide which root is ...
    (sci.stat.math)
  • Re: Multiple Regression w/ Polynomial-in-Y?
    ... Yields w-vector and c coefficient ... > for quadratic term for y. ... > terms of independent variables (which I suspect would not be so good ... One drawback is needing to decide which root is ...
    (sci.math.num-analysis)
  • Re: multicollinearity in regression
    ... I could use Analysis of Covariance but 2 of the independent variables ... I'm guessing that in the model with LOGSIZE, the LOGSIZE coefficient is ... multicollinearity: it may be that you can then see a sensible approach ... I always find it helpful to calculate the correlation coefficient ...
    (sci.stat.consult)
  • Re: multicollinearity in regression
    ... I could use Analysis of Covariance but 2 of the independent variables ... related to LOGSIZE variable which has a Variance Proportion of .99. ... I'm guessing that in the model with LOGSIZE, the LOGSIZE coefficient is ... I always find it helpful to calculate the correlation coefficient ...
    (sci.stat.consult)