Daniel@0: function CPD = set_params(CPD, varargin)
Daniel@0: % SET_PARAMS Set the parameters (fields) for a softmax_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: % (Let ns(i) be the size of node i, X = ns(X), Y = ns(Y), Q1=ns(dps(1)), Q2=ns(dps(2)), ...
Daniel@0: %   where dps are the discrete parents; if there are no discrete parents, we set Q1=1.)
Daniel@0: %
Daniel@0: % weights - (W(:,j,a,b,...) - W(:,j',a,b,...)) is ppn to dec. boundary
Daniel@0: %           between j,j' given Q1=a,Q2=b,... [ randn(X,Y,Q1,Q2,...) ]
Daniel@0: % offset  - (offset(j,a,b,...) - offset(j',a,b,...)) is the offset to dec. boundary
Daniel@0: %           between j,j' given Q1=a,Q2=b,... [ randn(Y,Q1,Q2,...) ]
Daniel@0: % clamped     - 'yes' means don't adjust params during learning ['no']
Daniel@0: % max_iter    - the maximum number of steps to take [10]
Daniel@0: % verbose     - 'yes' means print the LL at each step of IRLS ['no']
Daniel@0: % wthresh     - convergence threshold for weights [1e-2]
Daniel@0: % llthresh    - convergence threshold for log likelihood [1e-2]
Daniel@0: % approx_hess - 'yes' means approximate the Hessian for speed ['no']
Daniel@0: %
Daniel@0: % e.g., CPD = set_params(CPD,'offset', zeros(ns(i),1));
Daniel@0: 
Daniel@0: args = varargin;
Daniel@0: nargs = length(args);
Daniel@0: glimsz = prod(CPD.sizes(CPD.dpndx));
Daniel@0: for i=1:2:nargs
Daniel@0:   switch args{i},
Daniel@0:    case 'discrete',     str='nothing to do';   
Daniel@0:    case 'clamped',      CPD = set_clamped(CPD, strcmp(args{i+1}, 'yes'));
Daniel@0:    case 'max_iter',     CPD.max_iter = args{i+1};
Daniel@0:    case 'verbose',      CPD.verbose = strcmp(args{i+1}, 'yes');
Daniel@0:    case 'max_iter',     CPD.max_iter = args{i+1};
Daniel@0:    case 'wthresh',      CPD.wthresh = args{i+1};
Daniel@0:    case 'llthresh',     CPD.llthresh = args{i+1};
Daniel@0:    case 'approx_hess',  CPD.approx_hess = strcmp(args{i+1}, 'yes');
Daniel@0:    case 'weights',      for q=1:glimsz, CPD.glim{q}.w1 = args{i+1}(:,:,q); end; 
Daniel@0:    case 'offset',
Daniel@0:     if glimsz == 1
Daniel@0:       CPD.glim{1}.b1 = args{i+1};
Daniel@0:     else
Daniel@0:       for q=1:glimsz, CPD.glim{q}.b1 = args{i+1}(:,q); end; 
Daniel@0:     end
Daniel@0:    otherwise,  
Daniel@0:     error(['invalid argument name ' args{i}]);       
Daniel@0:   end
Daniel@0: end