RegulatoryModulesStringKernel.h

Go to the documentation of this file.

/*
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 3 of the License, or
 * (at your option) any later version.
 *
 * Written (W) 2009 Sebastian J. Schultheiss and Soeren Sonnenburg
 * Copyright (C) 2009 Max-Planck-Society
 */

#ifndef _REGULATORYMODULESSTRINGKERNEL_H___
#define _REGULATORYMODULESSTRINGKERNEL_H___

#include "lib/common.h"
#include "kernel/StringKernel.h"
#include "features/SimpleFeatures.h"

namespace shogun
{
class CRegulatoryModulesStringKernel: public CStringKernel<char>
{
    public:
        CRegulatoryModulesStringKernel(void);

        CRegulatoryModulesStringKernel(int32_t size, float64_t width, int32_t degree, int32_t shift, int32_t window);

        CRegulatoryModulesStringKernel(CStringFeatures<char>* lstr, CStringFeatures<char>* rstr, 
            CSimpleFeatures<uint16_t>* lpos, CSimpleFeatures<uint16_t>* rpos, 
            float64_t width, int32_t degree, int32_t shift, int32_t window, int32_t size=10);

        virtual ~CRegulatoryModulesStringKernel();

        virtual bool init(CFeatures* l, CFeatures* r);

        virtual EKernelType get_kernel_type() { return K_REGULATORYMODULES; }

        inline virtual const char* get_name() const { return "RegulatoryModulesStringKernel"; }
        
        void set_motif_positions(
            CSimpleFeatures<uint16_t>* positions_lhs, CSimpleFeatures<uint16_t>* positions_rhs);
        
    protected:
        virtual float64_t compute(int32_t idx_a, int32_t idx_b);
        
        float64_t compute_wds(char* avec, char* bvec, int32_t len);

        
        void set_wd_weights();

    protected:
        float64_t width;

        int32_t degree; 
        int32_t shift;
        
        int32_t window;

        CSimpleFeatures<uint16_t>* motif_positions_lhs;

        CSimpleFeatures<uint16_t>* motif_positions_rhs;

        float64_t* position_weights;

        float64_t* weights;
};
}
#endif /* _REGULATORYMODULESSTRINGKERNEL_H__ */