 function [co] = KEJR1_HR_initialization(mult,post_init)
%function [co] = KEJR1_HR_initialization(mult)
%function [co] = KEJR1_HR_initialization(mult,post_init)
%Initialization of the exponentiated JensenRenyi kernel1 estimator defined according to the relation:
%K_EJR1(f_1,f_2) = exp[u x H_R((y^1+y^2)/2)], where H_R is the Renyi entropy, (y^1+y^2)/2 is the mixture of y^1~f_1 and y^2~f_2 with 1/21/2 weights, u>0.
%
%Note:
% 1)The estimator is treated as a cost object (co).
% 2)We use the naming convention 'K<name>_initialization' to ease embedding new estimators for kernels on distributions.
% 3)This is a meta method: the Renyi entropy estimator can arbitrary.
%
%INPUT:
% mult: is a multiplicative constant relevant (needed) in the estimation; '=1' means yes, '=0' no.
% post_init: {field_name1,field_value1,field_name2,field_value2,...}; cell array containing the names and the values of the cost object fields that are to be used
% (instead of their default values). For further details, see 'post_initialization.m'.
%OUTPUT:
% co: cost object (structure).
%Copyright (C) 2013 Zoltan Szabo ("http://www.gatsby.ucl.ac.uk/~szabo/", "zoltan (dot) szabo (at) gatsby (dot) ucl (dot) ac (dot) uk")
%
%This file is part of the ITE (Information Theoretical Estimators) Matlab/Octave toolbox.
%
%ITE 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.
%
%This software is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of
%MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
%
%You should have received a copy of the GNU General Public License along with ITE. If not, see <http://www.gnu.org/licenses/>.
%mandatory fields:
co.name = 'EJR1_HR';
co.mult = mult;
%other fields:
co.u = 1; %assumption: u>0
co.alpha = 0.95; %assumption: 0<=alpha<1
co.member_name = 'Renyi_kNN_k'; %you can change it to any Renyi entropy estimator
%post initialization (put it _before_ initialization of the members in case of a meta estimator):
if nargin==2 %there are given (name,value) cost object fields
co = post_initialization(co,post_init);
end
%initialization of the member(s):
co.member_co = H_initialization(co.member_name,mult,{'alpha',co.alpha}); %{'alpha',co.alpha}: the 'alpha' field of the member is also set here
