wolffd@0: function CPD = set_fields(CPD, varargin) wolffd@0: % SET_PARAMS Set the parameters (fields) for a tabular_CPD object wolffd@0: % CPD = set_params(CPD, name/value pairs) wolffd@0: % wolffd@0: % The following optional arguments can be specified in the form of name/value pairs: wolffd@0: % wolffd@0: % CPT, prior, clamped, counts wolffd@0: % wolffd@0: % e.g., CPD = set_params(CPD, 'CPT', 'rnd') wolffd@0: wolffd@0: args = varargin; wolffd@0: nargs = length(args); wolffd@0: for i=1:2:nargs wolffd@0: switch args{i}, wolffd@0: case 'CPT', wolffd@0: if ischar(args{i+1}) wolffd@0: switch args{i+1} wolffd@0: case 'unif', CPD.CPT = mk_stochastic(myones(CPD.sizes)); wolffd@0: case 'rnd', CPD.CPT = mk_stochastic(myrand(CPD.sizes)); wolffd@0: otherwise, error(['invalid type ' args{i+1}]); wolffd@0: end wolffd@0: elseif isscalarBNT(args{i+1}) wolffd@0: p = args{i+1}; wolffd@0: k = CPD.sizes(end); wolffd@0: % Bug fix by Hervé Boutrouille 10/1/01 wolffd@0: CPD.CPT = myreshape(sample_dirichlet(p*ones(1,k), prod(CPD.sizes(1:end-1)), CPD.sizes)); wolffd@0: %CPD.CPT = myreshape(sample_dirichlet(p*ones(1,k), prod(CPD.sizes(1:end-1))), CPD.sizes); wolffd@0: else wolffd@0: CPD.CPT = myreshape(args{i+1}, CPD.sizes); wolffd@0: end wolffd@0: wolffd@0: case 'prior', wolffd@0: if ischar(args{i+1}) & strcmp(args{i+1}, 'unif') wolffd@0: CPD.prior = myones(CPD.sizes); wolffd@0: elseif isscalarBNT(args{i+1}) wolffd@0: CPD.prior = args{i+1} * normalise(myones(CPD.sizes)); wolffd@0: else wolffd@0: CPD.prior = myreshape(args{i+1}, CPD.sizes); wolffd@0: end wolffd@0: wolffd@0: %case 'clamped', CPD.clamped = strcmp(args{i+1}, 'yes'); wolffd@0: %case 'clamped', CPD = set_clamped(CPD, strcmp(args{i+1}, 'yes')); wolffd@0: case 'clamped', CPD = set_clamped(CPD, args{i+1}); wolffd@0: wolffd@0: case 'counts', CPD.counts = args{i+1}; wolffd@0: wolffd@0: otherwise, wolffd@0: %error(['invalid argument name ' args{i}]); wolffd@0: end wolffd@0: end wolffd@0: wolffd@0: