Log in

Experimental Gaussian kernel function for genetic optimized neural net

I worked recentrly on transfer functions in order to use the native genetic optimizer under Metatrader 4 platform. 

We are still at the experimental stage. 

1. The sample percetrone code with gaussian transfer function

A percetron code would look like this.

double perceptron() 
{
double w1 = MathExp(-x1*x1/2);// This is the Gaussian transfer function
double w2 = MathExp(-x2*x2/2);// Here we can specify the number if the hidden neurones
double w3 = MathExp(-x3*x3/2);
double w4 = MathExp(-x4*x4/2);

double p1 = iCustom(NULL,Symbol(),"PFE",x5,true,5,0,0);// This is the input
double p2 = iCustom(NULL,Symbol(),"PFE",x6,true,5,0,0);
double p3 = iCustom(NULL,Symbol(),"PFE",x7,true,5,0,0);
double p4 = iCustom(NULL,Symbol(),"PFE",x8,true,5,0,0);

return(w1*w2*w3*w4* p1 + w1*w2*w3*w4* p2 + w1*w2*w3*w4*p3 + w1*w2*w3*w4* p4);
}

// From this line depends the whole architecture of the percetrone. We  multiply the input by the value

//coming from the transfer function, in our example it is between 0 and 1.

So this is a sample code. Here we do not devide the output because we are interested if the value is above or below 0. 

It is also important that the inputs are normalized and scaled properly. The technical indicator PFE is very appropriate because it is properly scaled and normalized varying frome -1 to +1.

You can see a practical example of the scaling and normalization in the neural indicator code by Jaguar. 

Available for registered users.

 

The use of genetic optimizer limits the number of hidden neurones we can practically use. In the Neuroshell Trader as a commercial software a simmilar application limits itself in the neural net architecture only to two hidden neurones (). And that is considered as a good neural net implementation practice. (Neural indicators for Neuroshell are available as additional addon, you can refer to the help files).

 

 2. The transfer function 


As it was previously stated we can choose different transfer functions in this part of the code.

The gaussian transfer function is:

double w1 = MathExp(-x1*x1/2);

x1 means in practice that we are looking for the best avalues of x1 using the genetic optimizer.

the value of x would be from -3 tp +3. In practice I can suggest a step size of 0.01.

 

image
От 08 ноември 2011

 

3. Towards a more sophisticated kernel function.

Yesterday I watched a very impressive video about the possibilities of Metatrader 5 for visualisation of the 3 D visualisation space.

So then came into my mind the idea, as we are limited with the genetic optimizer by the number of the hidden neurones we can use (of course you can use many hidden neurones but you will otpimize forever).

So why not to use just one, or two hidden neurones but we will have two parameters in the transfer function, not just one.

So the code will look like this:

double w1 = MathExp(-x1*x1-y1*y1); 

Here we optimize two parameters x1 and x2.

On the first shot you can see what you will get if you do not use a genetic optimizer.


 

image
От 23 септември 2011

Here on the second shot is what you will practically get using a genetic optimizer searching for the best values of x and y.

For the variable x the start from -3, the step is 0.01 and the stop is 3

For the variable y the start from -3, the step is 0.01 and the stop is 3

So the idea is that we just one hidden neurone using a genetic optimization with two variables we could achieve 3 dimensional  hypersphere with just one radial unit.

image
От 23 септември 2011

Comments

  • forexnes 4211 days ago

    The theme is very interesting! This link may be helpful :)

    http://www.khanacademy.org/math/statistics/v/normal-distribution-excel-exercise

  • JohnLast 4211 days ago

    Thanks for the link.

    I am trying to do some other things today but I do not achieve more profitability.

    So as it is gaussian transfer function, I do not need to have it centered arround 0. Basically when I have two neurones:

    -one gaussian function will be cetered at +0.5 for example

    double w2 =  MathExp((-(x2-0.5)*(x2-0.5))/2);

     

    -the other will be centered at -0.5

    double w3 =  MathExp((-(x3+0.5)*(x3+0.5))/2);

     

    For the 3d gaussian transfer function it will be

    3d centered at +1

    double w1 =  MathExp(-(x1-0.5)*(x1-0.5)-y1*y1);

    3d centered at -1

    double w2 =  MathExp(-(x2+0.5)*(x2+0.5)-y2*y2);

     

    Accordingly the search space for the genetic has to be modified.

    -2.5 to 3.5 for the center 0.5

    -3.5 ro 2.5 for the center -0.5

  • JohnLast 4202 days ago

    I see also the Gauss kernel plottet like that.

     

    #property copyright "Copyright 2011, MetaQuotes Software Corp."
    #property link "http://www.mql5.com"
    #property version "1.00"
    //--- input parameters
    input double x=-3.0;
    input double y=-3.0;
    input double gamma=1;
    //+------------------------------------------------------------------+
    //| Tester function |
    //+------------------------------------------------------------------+
    double OnTester()
    {
    //---
    double res=MathExp(-gamma*(x*x+y*y -2*x*y ));
    //---
    return(res);

     

    image
    От Поле за пускане