[R] solving non-linear system of equations

From: HAKAN DEMIRTAS <demirtas_at_uic.edu>
Date: Tue 15 Aug 2006 - 01:30:02 EST

Didn't get any useful response to the following question. Trying again.

I can't seem to get computationally stable estimates for the following system:

Y=a+bX+cX^2+dX^3, where X~N(0,1). (Y is expressed as a linear combination of the first three powers of a standard normal variable.) Assuming that E(Y)=0 and Var(Y)=1, one can obtain the following equations after tedious algebraic calculations:

  1. b^2+6bd+2c^2+15d^2=1
  2. 2c(b^2+24bd+105d^2+2)=E(Y^3)
  3. 24[bd+c^2(1+b^2+28bd)+d^2(12+48bd+141c^2+225d^2)]=E(Y^4)-3

Obviously, a=-c. Suppose that distributional form of Y is given so we know E(Y^3) and E(Y^4). In other words, we have access to the third and fourth raw moments. How do we solve for these four coefficients? I reduced the number of unknowns/equations to two, and subsequently used a grid approach. It works well when I am close to the center of the support, but fails at the tails. Any ideas?

Hakan Demirtas

R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Tue Aug 15 01:34:55 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Tue 15 Aug 2006 - 04:22:18 EST.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.