function loglik=loglik(table,schaal)
    
e=size(table);
over=1:length(e)-1;
dims=[e(end) 1];
if ~isempty(over), marg_prob=franks_marginalize_pot(table,over,dims);
    else marg_prob=table;
end
loglik=sum(log(marg_prob)'+schaal);