function Pjj = getcondP1(P,m,n,jj,ii,y)
% Pjj(x) = P{ X[jj]=x | X[ii]=y }
if ~isempty(intersect(ii,jj))
  error('ii and jj must not intersect');
end

Pjj = getcondP(P,m,n,jj,ii,y);
if sum(Pjj) < .9
    if length(ii)>=1
        ii = ii(2:end);
        y = y(2:end);
        Pjj = getcondP1(P,m,n,jj,ii,y);
    else
        error('deficient probability')
    end
end
return