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