-
Notifications
You must be signed in to change notification settings - Fork 61
/
Copy pathprior_unif.m
67 lines (56 loc) · 1.34 KB
/
prior_unif.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
function p = prior_unif(varargin)
%PRIOR_UNIF Uniform prior structure
%
% Description
% P = PRIOR_UNIF creates uniform prior structure.
%
% See also
% PRIOR_*
%
% Copyright (c) 2009 Jarno Vanhatalo
% Copyright (c) 2010 Aki Vehtari
% This software is distributed under the GNU General Public
% License (version 3 or later); please refer to the file
% License.txt, included with the software, for details.
ip=inputParser;
ip.FunctionName = 'PRIOR_UNIFORM';
ip.addOptional('p', [], @isstruct);
ip.parse(varargin{:});
p=ip.Results.p;
if isempty(p)
init=true;
p.type = 'Uniform';
else
if ~isfield(p,'type') && ~isequal(p.type,'Uniform')
error('First argument does not seem to be a valid prior structure')
end
init=false;
end
if init
% set functions
p.fh.pak = @prior_unif_pak;
p.fh.unpak = @prior_unif_unpak;
p.fh.lp = @prior_unif_lp;
p.fh.lpg = @prior_unif_lpg;
p.fh.recappend = @prior_unif_recappend;
end
end
function [w, s, h] = prior_unif_pak(p, w)
w=[];
s={};
h=[];
end
function [p, w] = prior_unif_unpak(p, w)
w = w;
p = p;
end
function lp = prior_unif_lp(x, p)
lp = 0;
end
function lpg = prior_unif_lpg(x, p)
lpg = zeros(size(x));
end
function rec = prior_unif_recappend(rec, ri, p)
% The parameters are not sampled in any case.
rec = rec;
end