Serving the Quantitative Finance Community

 
User avatar
MikeCrowe
Topic Author
Posts: 0
Joined: January 16th, 2006, 8:20 am

Solve a noisy equation

November 3rd, 2006, 9:22 am

Suppose I have a noisy equation f(x) = a, and I want to solve this with as few samples from f(x) as possible. Can anyone suggest any good algorithms, or approaches/techniques for inventing a specific one.f(x) is a noisy approximation of an exact unknown function g(x). I know the following* g(x) is monotonically increasing and takes x [0,inf) and returns [0,1]* I know the likelihood distribution of g(x) given f(x)* equally I know the exact distribution of f(x) for any g(x)Any hints?
 
User avatar
zeta
Posts: 26
Joined: September 27th, 2005, 3:25 pm
Location: Houston, TX
Contact:

Solve a noisy equation

November 7th, 2006, 1:42 pm

How about kalman/recursive bayes on your series f(x) to estimate the coefficents in model g(x), which you solve analytically? Too cliche?
 
User avatar
zeta
Posts: 26
Joined: September 27th, 2005, 3:25 pm
Location: Houston, TX
Contact:

Solve a noisy equation

November 7th, 2006, 1:45 pm

or if your distributions are non-gauss, use EM algo to generate mixtures in terms of gauss, and run a series of kalman's?
 
User avatar
frontofficequant
Posts: 0
Joined: October 18th, 2006, 12:02 am

Solve a noisy equation

November 8th, 2006, 6:27 pm

There are several procedures for handling this. One of the most primative (yet most effective) is the "Robins-Monro" stochastic approximation algorithm.
 
User avatar
MikeCrowe
Topic Author
Posts: 0
Joined: January 16th, 2006, 8:20 am

Solve a noisy equation

November 9th, 2006, 8:39 am

Thanks zeta, but I was hoping not to model g(x) as this introduces extra errors, and I've not really got any justification for any given model.frontofficequant - that sounds interesting, could you point me in the direction of any good papers on this? What sort of algorithm is it? Direct search?The concern I've been having with most of the algorithms I've tried is that they cannot make use of the special information, such as monotonicity and knowing the conditional distributions, and so they seem very slow to converge.Additionally, what I really need is to be able to bound this from below, i.e. to find the largest x such that g(x) <= a, to some confidence, however the added twist is that i think f(x) generally converges on g(x) more slowly for smaller x.I'm hoping to do this in less than 100 steps.
 
User avatar
frontofficequant
Posts: 0
Joined: October 18th, 2006, 12:02 am

Solve a noisy equation

November 10th, 2006, 7:56 pm

The original paper is:Robbins, H. and Monro, S. (1951). A stochastic approximation method. Ann. Math. Stat. 22, 400-407. I'm sure there are more recent and refined version of this algorithm that apply to different situations. In fact, there are entire text books on stochastic approximation and it's convergence properties in several contexts.