construct_scpt.m

来自「The BNL toolbox is a set of Matlab funct」· M 代码 · 共 21 行

function smallCPTs=construct_sCPT(lin_pred,equiv_class,gausskwadnodes)
%smallCPTs=construct_sCPT constructs conditional probability tables 
%small: data not yet incorporated -> same CPTs for all persons

%does not always work (permute dimensions should be checked)

smallCPTs=cell(size(lin_pred));
for i=1:length(lin_pred)
    s=size(lin_pred{i});
    f=find(equiv_class==i,1);
    nodenr=find(gausskwadnodes.nodenrs==f,1);
    if ~isempty(f)
    [x,w]=herzo(gausskwadnodes.nodesizes(nodenr);%no rounding erros for guaaskwad nodes
    %by taking logit and transforming back to probs
    probs=w';
    else
    probs=multinom_logistic(lin_pred{i});
    end
    if isvector(probs) probs=probs(:);, end;
    smallCPTs{i}=probs;
end