[Lemon-user] user defined network optimum function for gene regulatory network optimization

[tomasz jankowski]

Hi All,

I'm looking some c++ library which help me find best regulatory gene network.
lets asume I have collected genotypes of three genes and length of mice tail in group of 1000 mices.
To find regulatory gene network which will describe  in most accurate way relations between genes in context of length of tail i need to have posibility to define my own function which calculate accuracy. Based on these function optimization algorith will be rearenging network topology until reach the best one.

I need to have posibility to asign values and operatos */+- to nodes and directed edges. (graph can star with many nodes but always end with only one) these values and operators will be used to calculate end value. node walue and oparator will depend on genotype. forexample if AA node has value 3 andoperator *. if BB value 4 and operator -.  values and operator of edge depend on genotypes of nodes between they are. for  example in three genes network if 1AA--->2AA--->3AA first edge has value 2 
and operator + but in 1AA--->2BB-->3AA first edge has value 0.5 and operator -.

after graph topology will be defined and values and operators for nodes and edges will be asigned based on genotypes. the final walue can be calculated . just by using values and operators from starting nodes to the end   node. 

after that will be calculated correlation in group of 1000 mices between length of mice tail and value obtained from regulatory network. 

the higher correlation the better network.

Can you please tell me if lemon library is what i'm looking for? will be possible to implement optimization of such a regulatory network with lemon? if yes, can you please give me some advices for best strategy for that task?

Thank you,
tomasz